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ABSTRACT 


This  report  describes  technical  findings  in  four  main  subject 
areas:  semiconductors  and  semiconductor  devices,  insulating  and  semi¬ 
conducting  glasses,  measurement  techniques,  and  magnetic  materials. 

For  silicon  doped  with  gold  and  phosphorus,  a  utudy  of  the  photo¬ 
magnetoelectric  effect  yields  a  basic  theoretical  relationship  between 
the  excess  hole  and  excess  electron  concentrations  and  an  experimental 
determination  of  the  effective  lifetime  for  temperatures  between  21 
and  84°K.  A  pragmatic  approach  is  proposed  and  empirically  tested  con¬ 
cerning  the  dependence  that  the  density  of  quantum  states  in  a  semicon¬ 
ductor  shows  on  the  impurity  concentration.  For  a  variety  of  ceramic 
materials,  the  electric-field  breakdown  shown  is  related  to  the  various 
environmental  and  material  variables.  Based  on  available  data  con¬ 
cerning  material  structure,  a  study  is  made  to  characterize  the  relation¬ 
ship  borne  in  semiconducting  glasses  between  heterogeneity  in  structure 
and  resultant  electrical  properties.  Errors  in  x-ray  spectrochemical 
analysis  of  powdered  materials  due  to  pressing  variations  are  quanti¬ 
tatively  examined.  The  several  means  for  analysis  of  the  defect  struc¬ 
ture  of  solid  surfaces  are  discussed  together  with  the  relation  that 
this  defect  structure  bears  to  the  surface  electronic  and  chemical 
properties.  Small-angle  x-ray  data  on  the  low-temperature  aging  of 
aluminum-zinc  alloys  is  shown  to  provide  a  method  for  determining  the 
particle  size  distribution  as  a  function  of  aging  time.  A  new  approach 
is  proposed  and  evaluated  for  characterizing  the  trapped-flux  inter¬ 
action  between  magnetic  films  and  nearby  electrical  conductors. 
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SUMMARY 


This  report  for  the  fifth  semiannual  period  of  contract  support 
describes  technical  findings  in  four  main  subject  areas:  semiconductors 
and  semiconductor  devices,  insulating  and  semiconducting  glasses,  measure¬ 
ment  techniques,  and  magnetic  materials. 

SEMICONDUCTORS  AND  SEMICONDUCTOR  DEVICES: 

For  silicon  doped  with  gold  and  phosphorus  the  photomagnetoelectric 
(PME)  effect  was  studied  theoretically  and  observed  experimentally  for 
temperatures  between  21  and  84#K.  A  generalized  diffusion  equation  is 
formulated  and  solved  with  the  impurity  centers  included  in  the  charge 
balance  equation  and  for  arbitrary  intensity  and  steady  Illumination. 

When  the  density  of  gold  greatly  exceeds  the  density  of  the  free  carriers, 
the  charged  impurities  maintain  charge  neutrality;  this  yields  a  range  of 
Injection  where  the  excess  hole  concentration  can  be  expressed  as  the  sum 
of  two  terms:  the  first  proportional  to  the  excess  electron  concentration 
and  the  second  to  the  square  of  that  concentration.  If  the  quadratic  term 
dominates,  the  PME  current  Is  proportional  to  the  4/3  power  of  the  photo¬ 
conductance.  Measurements  show  this  power-law  dependence,  from  which  one 
can  deduce  that  the  effective  lifetime  ranges  between  36  nsec  to  63  ysec. 

The  density  of  quantum  states  in  a  semiconductor  depends  on  the  impurity 
concentration;  a  study  is  reported  concerning  a  pragmatic  approach  to  the 
problem  of  describing  this  dependence.  The  approach  combines  heretofore 
separate  theories  concerning  the  various  contributors  to  this  dependence, 
utilizing  an  approximation  that  enables  application  of  the  superposed  theories 
over  the  range  of  impurity  concentration  of  technological  interest.  To 
demonstrate  the  practical  utility  of  the  approach,  the  prediction  yielded 
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by  its  application  to  p-type  gallium  arsenide  as  compared  to  an  experi¬ 
mentally  determined  density  of  states.  The  agreement  is  satisfactory. 

INSULATING  AND  SEMICONDUCTING  GLASSES: 

The  electric  field  a  ceramic  material  can  withstand  may  determine  its 
applicability  in  various  uses  in  the  electronic  technology:  as  an  insulator, 
as  a  protective  surface  in  various  devices,  as  an  Isolator  in  integrated  cir¬ 
cuits,  and  as  a  material  for  encapsulation.  The  several  theories  of  break¬ 
down  mechanisms  in  ceramics  are  reviewed  and  the  influence  of  environmental 
and  materials  variables  on  the  breakdown  strength  is  discussed.  The  discussion 
includes  consideration  of  the  following  variables:  thickness,  temperature, 
ambient  atmosphere,  electrode  shape  and  composition,  surface  finish,  field 
frequency  and  waveform,  porosity,  crystalline  anisotropy,  amorphous  structure, 
and  composition.  The  relative  importance  of  these  factors  receives  attention 
for  a  variety  of  ceramic  materials. 

Heterogeneities  present  in  semiconducting  glasses  play  a  key  role  in 
determining  the  electrical  properties  of  the  glass.  Based  on  available  data 
concerning  material  structure,  yielded  by  small-angle  scattering  x-ray  studies, 
and  concerning  the  dc  and  ac  conductivity  of  VjOj-KPO^  semiconducting  glass, 
a  theoretical  study  ia  made  to  characterize  the  relationship  borne  between 
structure  and  electrical  properties.  Various  alternative  models  of  the  con¬ 
duction  mechanisms  are  proposed  and  evaluated. 

Errors  in  x-ray  spectrochemical  analysis  of  powdered  materials  due  to 
sample  pressing  variations  are  quantitatively  examined.  It  is  shown 
that  the  x-ray  intensity  of  a  pressed  sample  is  a  function  of  the  sample  density. 
MgO,  Ti02,  Cr20j  and  Th02  samples  all  exhibited  three  major  regions  of  x-ray 
response  in  the  Intensity  vs  density  studies.  Analytical  errors  of  as  much 
as  SOX  can  result  if  measurements  are  made  without  regard  to  pressed  density 
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variation.  The  behavior  is  explained  in  terms  of  density  gradients  pro¬ 
duced  during  pression.  Appreciable  thickness  dependence  of  x-ray  results  in 
MgCO^  samples  is  also  discussed. 

MEASUREMENT  TECHNIQUES: 

The  development  of  ultra-high  vacuum  capabilities  and  sophisticated 
observational  methods  has  enabled  the  analysis  of  the  defect  structure  of 
solid  surfaces  and  the  relation  of  this  defect  structure  to  many  physical 
properties.  In  this  study  various  means  for  characterizing  the  defect  struc¬ 
ture,  such  as  optical,  x-ray  topographs,  and  transmission  electron-microscopy 
are  described  with  major  emphasis  being  placed  on  the  observation  and  inter¬ 
pretation  of  surface  defects  by  the  field-ion  microscopic  technique.  The 
influence  of  these  defects  on  the  surface  electronic  and  chemical  properties 
receives  attention. 

Small-angle  x-ray  data  on  the  low-temperature  aging  of  three  aluminum- 
zinc  alloys  affords  a  method  for  determining  the  particle  size  distribution 
as  a  function  of  aging  time.  The  study  yields  a  series  of  particle  size 
distribution  curves  and  particle  diameter  growth  paths,  and  the  evolution  of 
the  particle  size  distribution  can  be  followed.  The  observed  growth  can  be 
described  as  a  typical  particle  coarsening  process  with  competitive  growth. 
The  zonal  state  immediately  after  quench  is  characterized  by  a  more  mono- 
disperse  zone  size  distribution  than  at  later  aging  times. 

MAGNETIC  MATERIALS: 

The  trapped-flux  interaction  between  magnetic  films  and  nearby  elec¬ 
trical  conductors  is  a  well-known  phenomenon  experimentally  and  has  consti¬ 
tuted  the  basic  mechanism  of  NDRO  memories.  Previous  theoretical  analyses 
of  this  process,  however,  have  proved  so  formidable  computationally  and  con¬ 
ceptually  as  to  limit  their  usefulness.  A  new  approach  is  suggested  that. 


despite  its  simplicity,  leads  to  results  that  agree  well  with  those  that  are 
predicted  by  the  more  complicated  methods.  The  approach  consists  basically 
in  appealing  to  the  analogy  of  an  RC  transmission  line.  Subsequent  use  of 
this  analog  permits  geometrical  factors  to  be  represented  as  electric  cir¬ 
cuit  components  yielding  thereby  conceptual  and  computational  advantages. 
Taking  into  account  the  geometry  of  the  film  insulators  and  conductors,  the 
approach  provides  a  basis  for  estimating  the  magnitude  and  time  duration  of 
switching  transients.  The  model  developed  explains  the  experimental  obser¬ 
vation  that  eddy-current  torques  play  an  insignificant  role  in  film  flux 
reversal  for  devices  with  large  stripline  separation. 
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I.  Introduction 


The  original  general  objective  of  this  research  program  was  to  estab¬ 
lish  at  the  University  of  Florida  a  "Center  of  Competence  in  Solid-State 
Materials  and  Devices."  From  the  efforts  expended  in  developing  this 
center  of  competence  have  evolved  technical  findings:  technical  findings 
concerning  such  materials  as  glass  ceramics,  semiconducting  glasses,  mag¬ 
netic  films,  degenerate  materials,  and  degenerate  semiconductors,  concerning 
devices  made  from  these  materials;  concerning  measurement  techniques;  and 

concerning  methods  of  fabrication.  The  first  of  the  findings  in  these  various 

1-4 

areas  are  described  in  four  previous  scientific  reports.  To  Scientific 
Report  No.  1  the  reader  is  referred  for  a  more  detailed  statement  of  the 
research  objectives  than  given  here  and  for  a  discussion  of  the  means  to 
be  used  in  achieving  these  objectives.  The  present  report  sets  forth  major 
findings  of  the  fifth  semiannual  period  of  contract  support.  In  the  pre¬ 
sentation  to  follow,  Section  II  describes  the  results  of  research  concerning 
semiconductors  and  semiconductor  devices.  Section  III  reports  findings  con¬ 
cerned  with  insulating  and  semiconducting  glasses.  Section  IV  describes 
advances  made  in  measurement  techniques,  and  Section  V  deals  with  research 
concerning  magnetic  materials. 
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It*  Sam lconductors  and  S— iconductor  Devices  (A.  J.  Broderscn,  E.  R.  Chenettc, 

R.  W.  Gould,  L.  L.  Bench,  J.  J.  Bren,  S.  S.  LI,  F.  A.  Lind holm,  C.  T.  Sah, 

A.  van  der  Ziel) 

A.  LOU  TEMPERATURE  PBOTOMAGNETO ELECTRIC  PROPERTIES  OF  GOLD-DOPED  n-type 
SILICON  (J.  Agraz-G.  and  S.  S,  Li) 

1.  Introduction 

The  purpose  of  this  study  was  to  investigate  the  effects  of  heavy  gold 
doping  upon  the  transport  of  excess  carriers  in  n-type  silicon  by  means  of 
the  photomngnctoelcctric  (PME)  effect.  This  effect  was  observed  between  21* 
and  84*K  in  gold-  and  phosphorus-doped  silicon  and  the  small  Ball  angle  theory1 
has  been  extended  to  account  for  the  effects  of  the  impurities. 

For  constant  recombination  lifetime,  no  trapping  and  constant  ambipolar 
diffusivity,  the  PME  short  circuit  current  is  proportional  to  the  photocon¬ 
ductivity1.  The  PME  effect  has  been  studied  in  silicon  under  these  conditions 
for  small^  and  large  signal^’*.  The  trapping  by  the  gold  centers  has  signifi¬ 
cant  effects  on  the  diffusion  of  excess  carriers  when  the  equilibrium  carrier 
density  is  smaller  than  the  gold  density.  This  case  occurs  in  silicon  over- 
compensated  by  gold  at  room  temperature  and  in  undercompensated  silicon  at 
very  low  temperatures.  Our  experiments  at  low  temperature  show  a  PME  short 
circuit  current  proportional  to  the  4/3  power  of  the  photoconductance.  This  power 
law  is  also  obtained  from  our  theoretical  analysis. 

2.  Theory 

The  experimental  situation  is  schematically  shown  in  Fig.  1.  The  theory 
for  the  PME  effect1  involves  the  solution  of  the  continuity  equation  for  the 
carriers  injected  at  the  illuminated  surface.  We  have  obtained  this  solution 
under  the  following  assumptions: 

a.  The  magnetic  flux  density  is  very  small. 

b.  The  carriers  diffuse  in  one  direction  only. 

c.  Photoinjection  and  surface  recombination  can  be  represented  as 
surface  boundary  conditions. 

d.  Charge  neutrality  is  preserved  and  results  in  a  relationship  between 
the  excess  densities  of  elections  and  holes. 
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e.  The  recombination  rate  is  a  function  of  injection  only. 

f.  The  injected  carriers  do  not  reach  the  dark  surface. 

2.1  Charge  Neutrality 

For  charge  neutrality,  the  injection  dependent  densities  must  satisfy 

the  equation:  + 

p  +  N,  +  Nn  -  n  +  N.  +  n.  (1) 

r  Au  D  Au  D 


A  detailed  solution  of  this  equation  is  reported  elsewhere  .  Here  we 
only  consider  the  results  of  direct  interest  to  our  experimental  work  on  the 
PME  current  and  photoconductivity  measurements. 

Shockley-Read 's^  (S-R  model)  and  Sah-Shockley's2(S-S  model)  statistics 
are  used  to  relate  the  densities  of  the  electrons  trapped  in  the  impurity 
centers  to  the  densities  of  free  carriers.  A  relationship  between  the  excess 
carrier  densities  £p  and  An  can  be  deduced  from  Eq.  (1)  by  making  use  of 
S-R  and  S-S  statistics,  which  takes  the  simple  form  of  power  laws  within 
certain  ranges  of  injection^.  Of  particular  interest  here,  are  the  low 
temperature  results.  The  low  and  moderate  injection  solution  of  Eq.  (1) 
takes  the  form 


where 


and 


This  solution  is  valid  for 


Ap 

"  ri  An 

+  r2  An2 

(2) 

ri 

-  h>- 

NAul/(Vl/2> 

(3) 

r2 

n 

*-j 

*-■ 

a 

o 

f0r  ND  >  NAu 

(A) 

r2 

■  VWl/2  for  "d  ‘  V 

(5) 

for 

An  <<  (I 

r  -1 

(6) 

0  < 

2  JYl/2Y-l/2* 

and  can  be  separated  into  c  small  injection  linear  range  (i.e.:  Ap 


TjAn) 
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and  an  intermediate  quadratic  range  (i.e.:  Ap  *  r^An*-). 

This  analysis  also  shows  that  the  electron  recombination  lifetime 
changes  slowly  with  injection? 

2.2  Diffusion  Equation 

The  PME  current  is  related  to  the  diffusion  current,  and  the  photocon¬ 
ductivity  is  proportional  to  the  average  injected  carrier  density.  The  con¬ 
tinuity  equation  for  electrons  is  given  by 


where 


+  qR 


0 


qnp  E  +  qD 
^  n  y  s  n  dy 

qpu  E  -  qD 
P  y  p  dy 


(7) 

(8) 

(9) 


The  Dember  field  (E^)  in  Eqs. 
current  in  the  y-direction  is 
can  be  reduced  to: 


(8)  and  (9)  is  independent  of  y 

zero.  With  these  considerations 

.  ^  dn 

J  •  qD  7 

ny  dy 


and  the  total 
,  Eqns.  (8)  and  (9) 

(10) 


where 


D-D 


p  +  n 


& 

dn 


n 


p  +  bn 


,  and  p 


p  +  Ap,  n  -  n  +  An 
o  o 


(ID 


is  a  generalized  diffusivity  function.  The  continuity  equation  in  Eq.(7) 
then  takes  the  form: 

57  +  *■  0  <12> 

The  solution  of  Eq.(12)  gives  an  implicit  expression  for  the  carrier 

profile:  An 

y  -  Y  -  /  ( -D""d-An - nr-)  (13) 

AUy  (  2/DRdAn  C1) 
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The  constants  of  integration  are  determined  by  the  boundary  conditions  at 
the  illuminated  and  dark  surfaces.  If  the  injected  carriers  do  not  reach  the 
back  surface,  the  constant  of  integration  is  zero.  This  can  be  shown  from 
the  boundary  condition  at  this  surface: 


ny!y 


+  S .  An 
d  Y 


0 


(14) 


Usually  the  computation  of  Eq.  (13)  is  only  an  intermediate  step  in 
the  calculation  of  a  measurable  quantity.  Therefore,  a  complete  evaluation 
of  this  solution  is  not  necessary. 


2.3  The  PME  Effect  and  Photoconductance 

The  PME  short  circuit  current,  in  the  limit  of  small  magnetic  flux 
density  is:*  y 

*PME  "  ~e  '  Jny  d* 

where  0  -  |0j  +  |0p|  ,  0R-  -u^B,  0p  -  ypHB  . 


(15) 


For  the  thick  sample,  using  Eq.(10)  we  can  write 


PME 


An0 

0  /  q  D  dAn 
0 


(16) 


Similarly  the  photoconductance  can  be  written  as: 

Y 

AG  -  qy  /  (An  +  r^)  dy  (17) 

n  o  D 

°r  An. 

AG  -  qy  /  (An  +  rfy*  - -  (18) 

"0  b  4n  m 

[2  /  DR  dAn)1' ^ 

0 

Eqs . (16)  and  (18)  relate  IpME  and  G  parametrically  through  An^.  In  the  low 
injection  linear  range  Ap  ■  T^An,  Solving  Eqs. (10),  (11),  (16)  and  (17)  one 
finds 


6 


e(D./t  ) 
1  n 


1/2 


PME  p  (1  *  F./b) 
n  x 


A  G 


°q  DiQ0 

Si+  VTn 


(19) 


where 


(n  +  An)r. 

0  ■  D 

1  n  bn  +  T,  An 


(20) 


Eq.(19)  is  consistent  with  the  result  obtained  by  Amith 
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Eq.(19)  predicts  a  linear  relationship  between  and  AG.  In  this 
range,  the  photoconductance,  PME  current  and  the  photon  flux  density  are  pro¬ 
portional  to  the  first  power  of  the  surface  injection  (and  therefore  to  each 
other).  The  effective  diffusivity  (D^)  and  the  lifetime  were  found  to  be  con¬ 
stant  in  thie  range.  The  linear  relationship  between  IpME  and  AC  has  been 

4 

observed  recently  by  Li  in  n-type  silicon  at  low  temperatures. 

In  the  intermediate  injection  range,  where  the  steady  state  trapping  of 
photo-injected  carriers  by  gold  centers  becomes  significant,  the  relation 
between  Ap  and  An  is  quadratic,  i.e., 


Ap  -  r2  An 


(21) 


By  substituting  Eq.(21)  into  (11)  and  realizing  the  fact  that  in  this 
region  Ap  >>  pQ,  An  >>  nQ  and  Ap  <  An,  the  diffusivity  function  D  in  Eq.(ll) 


reduced  to  the  form: 


3Dnr7 

D  -  — r —  An 


(22) 


Substituting  Eq.(22)  into  Eq.(16)  yields  the  PME  short  circuit  current: 

3D  r,q  7 

I _ -  e(  "-■■)  (An  ;  (23) 


PME 


An 


And  by  substituting  Eq.(22)  and  R  ■  —  into  Eq.(18)  one  finds 

n 


AG  »  |i  q  j2r,D_t_/b  (An  ) 
n  w  l  n  n  o 


3/2 


(24) 


From  Kq8.(23)  and  (24),  the  'relationship  between  I_._  and  AG  is  given 

rrlE 
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r 


in  the  folJuwiiu;  fm ut: 


1’ME 


,  r,D  1/3  Ar  4/3 

*  l'“  V 


(25) 


wliicli  predicts  a  4/3  power  law  between  Ip^  and  AG.  This  power  law  relation- 

sliip  was  observed  in  gold-doped  n-type  silicon  at  low  temperatures  as  will  be 

reported  in  the  next  seel  ion. 

3.  Experimental  Results  and  Discussion 

We  report  here  the  experimental  results  on  three  gold-doped  silicon 

samples.  The  samples  were  prepared  by  diffusion  of  gold  into  phosphorus-doped 

silicon  bars.  The  impurity  concentrations,  listed  in  Table  1,  were  estimated 

from  Hall  data  obtained  between  20®K  and  300°K. 

By  comparing  the  Hall  data  of  the  different  samples  and  the  gold  diffusion 

data,  we  have  arrived  at  values  for  the  phosphorus  and  gold  densities  as 

listed  in  Table  I.  The  activation  energy  listed  is  deviated  from  the  accepted 

value  of  0.044  eV  for  phosphorus  at  low  donor  concentrations.  However,  it  is 

9 

consistent  with  measurements  by  Long  and  Myers  of  0.043  eV  for  a  donor  con- 

35  -3  10 

centration  of  4.5  x  10  cm  and  by  Swartz  who  obtained  0.032  eV  for 
1 7  -3 

1.8  x  10  cm  of  donor  concentration.  The  activation  energy  of  the  phosphorus 

is  estimated  for  samples  1,2  and  3  from  the  slope  of  the  low  temperature  plot 
3/2 

of  R  T  versus  1/T,  which  yields  a  value  of  0.033  eV. 


To  check  the  adequacy  of  our  model,  we 


investigate  the  temperature 


dependence  of  the  1’MK  short  circuit  current,  Ip^t  and  compare  with  the  theo¬ 


retical  expression.  From  F,q.(25)  we  find  that 

I  «  f 
l'ME  2 


(26) 


This  is  the  factor  that  contributes  most  to  the  dependence  of  I  on  tem- 

PME 


perature  because 


Au 


NDYl/2nlD 


(27) 
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Table  1.  Impurity  concentrations  from  Hall  effect ,  conductivity  and  gold  diffusion  data 


where 


(28) 


n1D  =  Nc  exp[-(Ec-ED)/kT] 

Since  n^  decreases  exponentially  with  temperature  with  the  activation 
energy  of  the  phosphorus  level,  it  is  found  from  Eq.(26)  and  (27)  that  Ip^E 
should  increase  with  decreasing  temperature  with  one  third  of  the  activation 
energy  of  phosphorus.  The  observed  Ip^E  for  AG  =  0.1  mho/cm  in  sample  3  is 
plotted  in  Fig.  5.  The  slope  indicated  is  0.001  eV  which  is  about  one-third 
of  the  activation  energy  we  observed  for  phosphorus.  This  is  in  good  agree¬ 
ment  with  the  results  obtained  from  the  Hall  measurements.  The  experimental 

procedures  for  the  l’ME  and  photoconductance  measurements  were  described  in 
3  4 

detail  previously  *  . 

In  Sample  1,  gold  was  diffused  at  1200°C  and  annealed  at  800°C.  This 
resulted  in  a  light  compensation  by  gold.  This  sample  shows  normal  band  con¬ 
ductivity  in  the  liquid  hydrogen  temperatures.  The  mobility  increases  mono- 
tonically  with  decreasing  temperature  in  consistency  with  the  theory  for 
small  compensation.  The  PME  short  circuit  current  at  24.7°K  is  proportional 

to  the  4/3  power  of  the  photoconductance  which  is  consistent  with  our  theory 

2 

for  the  range  of  injection  where  Ap  <*  An  .  This  result  is  shown  in  Fig.  2. 

Sample  2  was  prepared  by  diffusing  gold  at  1200°C  and  then  quenching  to 
room  temperature.  Tito  compensation  of  phosphorus  by  gold  for  this  case  was 
about  30  percent.  Impurity  conduction  by  hopping  was  observed  below  22®K, 
mixed  impurity  and  bind  conduction  between  22°  and  32°K  and  normal  band  con¬ 
duction  above  32°K.  The  PME  current  follows  the  4/3  power  law  at  32°K,  but 
shows  anomalous  behavior  below  29.3°K  where  the  impurity  conduction  becomes 
important.  This  anomalous  behavior  is  attributed  to  the  interaction  of  in¬ 
jected  carriers  with  impurity  conduction. 
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Weber/m 


Fig.  3  PME  short  circuit  current  versus  photocon¬ 
ductance  for  sample  2. 

N_  =  1.8  x  1017cm-3,  N4  -  8  x  1016cm"3. 

D  Au 
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In  sample  3,  gold  was  diffused  at  1200°C  and  quenched  to  room  tempera¬ 
ture.  The  gold  overcompensation  observed  in  this  sample  establishes  the  gold 

16  -3 

concentration  by  this  diffusion  method  ut  about  5  x  10  cm  .  The  4/3  pover 
law  was  observed  between  20.8*  and  84*K  in  good  agreement  with  our  model.  (Sec 
Eq.<23).) 

A  nonlinear  PME  effect  in  gold-doped  silicon  has  been  observed  In  the 
injection  range  characterized  by  a  4/3  power  law  for  the  PME  current  depen¬ 
dence  on  the  photoconductancc.  This  power  law  dependence  is  related  to  the 
steady  state  trapping  of  the  photo-injected  carriers  by  the  gold  centers.  From 
Eq. (23)  and  the  results  shown  in  Figs.  2,  3,  and  4,  one  can  deduce  the  electron 

effective  lifetime  J  v,  ._  t  .  The  corresponding  values  are  listed  In  Table  2. 

“  1/2  n 

Since  no  data  are  available  in  literature  for  the  ratio  of  electron  and 

hole  capture  rate  at  the  gold  acceptor  centers  at  low  temperatures ,  it  is 

impossible  to  calculate  the  electron  lifetime  t  from  the  PME  and  PC  measure- 

n 

ments  for  the  present  case.  However,  it  is  possible  to  determine  from 
photoconductivity  decay  experiment  at  low  temperature  and  then  calculate  the 
*1/2  from  the  PME  and  PC  measurements.  The  increase  in  temperature  of 
for  sample  3  is  consistent  with  the  fact  that  Yj^  lncrea8es  with  temperature 

*t*  X 

because  C^  increases  faster  with  temperature  than  Cr. 

The  PME  open  circuit  voltage  observed  in  the  gold-doped  silicon  is  found 

several  orders  of  magnitude  greater  than  that  observed  in  samples  without 
2  3  4  , 

gold  '  '  .  The  formulation  for  the  steady  state  diffusion  equation  shown  here 
is  general  with  respect  to  trapping  and  recombination  processes. 


IS 


r 


Tabic*  2.  The  Jh72  product . 


T(#K) 

Kn  Tn  (U8CC) 

Sample  1 

24.7 

63 

Sample  2 

30.1 

0.71 

Sample  3 

84.0 

0.036 

30.4 

0.075 

26.0 

0.098 

21.6 

0.362 

20.8 

0.448 

.  

,  . .  .  ...  ...  T  ri  — 
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LIST  or  SYMBOLS 


p,  n 
Ap,  An 
n 


D 


nau’  nd 
Yl/2’  Y-l/2 


R 

S, 


Ally,  Ally 


n 


PME 

AC-' 


n 


II) 


Total  density  of  holes  and  electrons,  respectively 
Excess  density  of  holes  and  electrons,  respectively 
Equilibrium  electron  density 

Density  of  electrons  in  the  shallow  donor  level 

Density  of  positively  and  negatively  charged  gold  centers, 
respectively 

Total  density  of  gold  and  shallow  donor  impurities,  respectively 

B  Ratio  of  the  hole  and  electron  capture  rate  for  the  gold 
acceptor  and  donor  levels,  respectively 

Recombination  rate 

Dark  surface  recombination  velocity 

Thickness  of  the  sample 

Excess  electron  density  at  the  illuminated  and  dark  surface, 
respectively. 

Electron  recombination  lifetime 
Electron-hole  mobility  ratio 
PME  current  per  unit  sample  width 

Photoconductaucc  of  unit  width  -  unit  length  sample 
Nc  cxp[-(Ec-K0)/kT) 

Illuminated  surface  recombination  velocity 
Photon  flux 
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B.  IMPURITY  CONCENTRATION  DEPENDENCE  OF  THE  DENSITY  OF  STATES  IN  SEMICONDUCTORS 
(D.  Dale  Kleppinger  and  F.  A.  Lindholm) 

I.  Introduction 

(1  2  3) 

The  conventional  engineering  theory  of  semiconductors  ’  *  has  re¬ 
ceived  wide  use  in  describing  the  behavior  of  semiconductor  devices.  This 
theory  is  based  on  various  approximations  that  make  design  and  analysis 
manageable;  among  these  approximations  is  the  neglect  of  the  dependence  of 
the  density  of  quantum  states  on  impurity  concentration.  In  predicting  the 
behavior  of  many  semiconductor  devices ,  this  neglect  Introduces  no  serious 
error;  but  for  others,  in  which  highly-doped  regions  play  a  key  role,  it  may 
cause  the  theory  to  fail  utterly  in  portraying  certain  aspects  of  device 
performance. 

That  the  density  of  quantum  states  does  indeed  depend  on  impurity  con¬ 
centration  is  implied  by  the  following  experimental  evidence.  In  their  study 

(A) 

of  silicon,  Pearson  and  Bardeen  found  that  the  impurity  ionization  energy 
decreased  with  increasing  impurity  concentration.  Debye  and  Conwell^  found 
the  same  behavior  in  germanium.  Conwell ^  described  measurements  of  car¬ 
rier  concentration  and  mobility  versus  temperature  and  impurity  concentration 
in  silicon  and  germanium.  To  explain  the  data,  she  used  a  theory  that  sup¬ 
posed  the  existence  of  two  species  of  mobile  carriers:  the  conventional  free 
carriers  and  carriers  in  states  associated  with  the  impurities.  Swartz^ 
studied  resistivity  in  silicon  as  a  function  of  temperature  and  impurity  con¬ 
centration.  He  found  that,  depending  on  temperature  and  impurity  concentra¬ 
tion,  three  modes  of  conduction  can  exist:  normal  band  conduction  in  the 
states  of  the  host  lattice  band;  non-band  conduction  (hopping)  in  states 
associated  with  the  impurities;  and  Impurity  band  conduction,  again  in  states 
associated  with  the  impurities.  Finally,  measurement  of  the  optical  absorp- 
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tion  coefficient  for  various  impurity  concentrations  gives  direct  evidence 

that  the  band  gap,  and  hence  the  density  of  quantum  states,  depends  on  im- 

(8) 

purity  concentration'  . 

This  dependence  also  manifests  itself  in  device  characteristics.  By 
regarding  the  density  of  quantum  states  to  depend  on  impurity  concentration, 

(9) 

Kane  has  accounted  for  the  non-zero  valley  current  shown  by  tunnel  diodes. 
Kauffman  and  Bergh^^  and  Buhanan^)  have  explained  the  observed  anomalous 
temperature  sensitivity  of  transistor  gain,  h__,  as  being  due  to  a  smaller 
band  gap  in  the  heavily  doped  emitter  than  exists  in  the  more  lightly  doped 
base. 

The  present  paper  develops  an  approach  that  enables  quantitative  account 
to  be  taken  of  the  relationship  between  the  impurity  concentration  and  the 
density  of  quantum  states.  The  approach  involves  superposing  heretofore 
separate  theories  concerning  the  various  contributors  to  this  dependence. 
Combined  with  a  key  approximation,  which  is  justified,  about  the  extent  of 
their  validity,  the  superposed  theories  afford  a  description  for  the  whole 
range  of  impurity  concentration  that  is  of  technological  interest. 

Section  II  below  sketches  some  aspects  of  the  conventional  theory.  In 
tills  section  and  subsequently,  we  treat  as  an  example  n-type  material,  unless 
otherwise  indicated;  the  same  treatment  holds  also,  however,  for  p-type  mate¬ 
rial.  Section  III  focuses  on  the  density  of  quantum  states  arising  from  the 

impurity  atoms;  and  Section  IV  focuses  on  the  density  of  states  in 
t he  host  lattice  band.  Section  V  sets  forth  and  justifies  the  approach  by 
which  we  combine  these  contributors  to  the  density  of  states  to  yield  a  com¬ 
prehensive  description  of  the  dependence  on  impurity  concentration.  To 
demonstrate  the  practical  utility  of  the  approach,  a  prediction  based  on  its 
use  for  p-type  gallium  arsenide  is  compared  in  Section  VI  to  an  experimental 
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observation.  The  approach  holds  also  for  material  other  than  gallium  arsenide, 
and  Section  VII  indicates  some  results  to  be  included  in  future  papers. 

II.  Convent ijnal  Theory 

The  conventional  engineering  theory  of  semiconductors  makes  three  assump¬ 
tions  about  the  density  of  quantum  states:  The  density  of  states  due  to  the 
host  lattice  conduction  band  is  proportional  to  the  square  root  of  energy; 
the  density  of  states  due  to  the  impurity  atoms  is  a  delta  function  positioned 
in  the  energy  gap  and  separated  from  the  edge  of  the  conduction  band  by  the 
ionization  energy  of  the  Impurity;  and,  there  is  no  impurity  concentration 
dependence  for  either  of  these  densities  of  states.  Let  us  briefly  examine 
these  assumptions. 

The  assumption  that  the  conduction  band  density  of  states  is  proportional 
to  the  square  root  of  energy  is  valid  only  near  the  bottom  of  the  conduction 
band  and  only  for  intrinsic  material  or  for  material  with  small  or  moderate 
impurity  concentration.  (The  appendix  contains  a  discussion  of  such  relative 
terms  as  "small"  and  "moderate.")  At  large  impurity  concentrations  the  elec¬ 
trons  donated  to  the  crystal  become  so  numerous  that  their  interaction  with 
the  atomic  cores  alters  the  conduction  band.  The  alteration  takes  the  form 
of  a  deviation  from  the  square  root  dependence.  This  concentration  dependent 
effect  is  sometimes  called  the  formation  of  "band  tails." 

For  the  contribution  to  the  density  of  states  due  to  the  impurity  atoms, 
the  interactions  between  the  individual  impurity  atoms  increase  as  the  con¬ 
centration  increases,  i.e.,  as  the  average  inter-impurity  distance  decreases. 

The  result  of  these  interactions  is  like  that  occurring  when  many  atoms  are 

(12) 

brought  closer  together  in  a  conceptual  development  of  crystal  properties.  .  The 
energy  levels  of  the  individual  atoms  split  and  form  a  quasi-continuous  band  of 
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energies,  the  "Impurity  band 


II 


The  approach  to  be  described  provides  a  means  for  simultaneous  inclusion 
of  the  concentration  dependence  of  the  band  tails  and  the  Impurity  band. 

III.  Impurity  States 

For  small  impurity  concentrations,  implying  negligible  interaction 
between  impurity  atoms,  the  energy  dependence  of  the  impurity  states  is  ade¬ 
quately  described  by  a  delta  function.  As  the  impurity  concentration  in¬ 
creases,  interaction  begins  to  become  important  and  the  (spatially  averaged) 
density  of  impurity  states  begins  to  spread.  Most  of  the  states  still  reside 

at  the  non-degenerate  level,  and  the  spread  occurs  around  this  energy.  This 

(13) 

approach  is  based  on  the  calculations  of  James  and  Ginzberg  which  dealt 
with  a  one-dimensional,  random  lattice.  They  found  that  the  density  of  states 
is  nearly  Gaussian  in  form.* 

For  a  three-dimensional  crystal,  the  logical  extension  is  to  assume  that 

the  density  of  impurity  states  is  still  described  by  a  Gaussian  whose  mean 

is  given  by  the  non-degenerate  ionization  energy  and  whose  standard  deviation 

depends  on  impurity  concentration.  This  function  will  be  essentially  a  delta 

function  at  small  concentrations  and  will  spread  around  the  non-degenerate 

(14) 

ionization  energy  as  the  impurity  concentration  increases.  Morgan  has 
developed  a  Gaussian  function  for  describing  the  impurity  states  that  we  will 
use  together  with  our  assumption  which  places  the  mean  of  the  distribution 
at  the  ionization  energy. 

Morgan's  equation  for  the  density  of  impurity  states  is 

*  Since  the  impurity  atoms  are  assumed  to  be  randomly  distributed,  the 
impurity  band  resulting  from  their  interactions  will  vary  in  width  in  energy 
throughout  the  material.  Spatial  averaging  of  this  varying  width,  necessary 
for  comparison  to  experiments,  results  in  an  impurity  band  that  does  not  have 
sharp  edges  but  "tails  off." 
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(1) 


p(E)  -  (2N) (2ir)  2  o'1  exp[-^-]  , 

2o 

where  ^111 

o  -  (^y-)2  A2N  2  exp(-  (6.2419xlOU)  (2) 

£ 

and  A  is  the  screening  length.  Equation  (1)  has  been  normalized  so  that 

the  doubly  Infinite  energy  Integral  of  p(E)  gives  the  total  number  of  impurity 

states,  2N  (including  spin  degeneracy).  The  equations  have  been  converted 

from  the  cgs-esu  system  of  units  used  by  Morgan  to  a  modified  cgs-esu 

system  of  units  in  Which  energy  is  measured  in  electronvolts. 

Using  the  screened  Coulomb  potential  screening  length  for  degenerate 

material,  1.  *  i  A 

A  -  [2(f)6  (S-)2  ^N6]"1  ,  (3) 

we  obtain  for  Eq.(l) 

1 

p(E)  -  3. 651x10®  (  e)  4  (~)  4  N,5833ex 

The  units  of  p(E)  are  reciprocal  electronvolts  per  cubic  centimeter;  1i  is 
in  erg-sec.  Equation  (A)  has  been  derived  assuming  that  tq  «  X  or,  essen¬ 
tially,  rQ  ■  0.  The  quantity  rQ  is  a  radius  which  Morgan  uses  to  define  a 
sphere  near  r  ■  0  where  no  impurity  sites  are  located.  Its  value  is  usually 
approximately  one-half  of  the  lattice  constant.  Although  Morgan  points  out 
that  can  be  used  to  include  skewness  in  the  distribution  function,  we  will 
consider  only  a  pure  Gaussian,  leaving  for  future  study  the  question  of  the 
effects  of  skewness. 

IV.  Conduction  Band  States 

As  we  noted  previously,  near  the  bottom  of  the  conduction  band  the  den¬ 
sity  of  quantum  states  in  the  conduction  band  is  proportional  to  the  square 
root  of  energy  —  but  only  if  the  material  is,  at  most,  lightly  doped.  At  large  Impu- 


-(E-y)2(e)2(g-)21 

9.551x10’18N,8334J 
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rity  concentration,  conduction  electrons  become  so  numerous  that  they  alter 
the  density  of  states  through  their  interactions.  The  main  effect  resulting 
from  the  interactions  is  the  screening  of  the  field  contributed  by  the  atomic 
cores,  thereby  altering  the  periodic  potential  distribution  of  the  system 
and  hence  the  allowed  energies.  Bonch-Bruyevich has  developed  a  theory 
which  considers  this  screening  and  describes  the  associated  density  of  states 
in  heavily  doped  semiconductors. 

Bonch-Bruyevich  deals  with  semiconductors  with  very  large  impurity  con¬ 
centration,  i.e.,  semiconductors  in  which  the  impurity  band  and  the  conduction 
band  have  overlapped.  The  equation  for  the  density  of  states  in  such  a  mate¬ 
rial  is  *  2^ 

p(E)  “  ^$”7 3  +  I  X(x)  +  I3(x)  +  I  3(x)]  (5) 

16tt  an  ~r  -7*  7 *  -r 


where 


(6) 


a 


7_  3  l  5_ 

12  2  x2  ..12 

?r _ e  n  N 

1  1_  3  1 

,2  ,12  4  4 

2  3  cm 


(7) 


and  the  I's  are  Bessel  functions  with  imaginary  argument.  By  making  appro¬ 
priate  mathematical  approximations,  Bonch-Bruyevich  reduces  this  equation 
to  two  simpler  equations  which  apply  for  particular  ranges  of  energy.  These 
equations  are  as  follows. 

For  energies  near  the  Fermi  level  the  equation  becomes 

2 

p(E)  -  p(l  -  S-r)  (8) 

0  4E2 

where  pq  is  the  usual  parabolic  density  of  states  function, 
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3  1 

*2  2  2  i3  -1 

Pq(E)  =  [  (2m  ;  EZ]  [ 2ttZ  A  ]  . 


Equation  (8)  shows  that  near  the  Fermi  level  the  density  of  states  differs 
only  slightly  from  that  in  lightly  doped  material.  It  should  be  recalled 
that  in  heavily  doped  semiconductor  material  the  Fermi  level  is  expected  to 
be  in  the  allowed  band.  Note  that  at  higher  energies  the  density  of  states 
approaches  the  lightly  doped  form. 

The  second  simpler  equation  applies  for  energies  near  the  conduction 
* 

band  edge.  The  density  of  states  in  this  energy  range  will  be  altered  to 
a  greater  degree  as  shown  by  Eq.  (10)^^. 

P (E)  -  [(2m*)2  a2]  [2tt2  t>3  T Cr) ] _1  [  1  +  “1  .  (10) 

*r(f) 


This  equation  is  valid  for  energies  greater  than  -E^ ,  where 

1_ 

_  r  *  4, ro  2  *2.-1  3.10 

E,  =  [m  e  ] [2e  n  ]  (Na  ) 

1  o 


and 


a 

o 


*  2 
m  e 


(U) 

(12) 


is  the  first  Bohr  radius  in  the  crystal.  Equation  (10)  gives  a  linear 
dependence  on  energy  for  the  density  of  states  near  the  band  edge,  rather 
than  the  conventional  square-root  dependence. 

As  has  been  pointed  out  above,  Bonch-Bruyevich's  equations  apply  to 


*  Since  we  will  be  discussing  heavily  doped  material  which  has  its  host  lat¬ 
tice  band  distorted  from  the  intrinsic  or  lightly  doped  case,  it  is  necessary 
to  establish  a  convention.  The  term  "band  edge"  will  refer  to  the  ene  :gy  at 
which  the  density  of  states  becomes  zero  in  non-degenerate  or  lightly  doped 
material.  When  it  becomes  necessary  to  discuss  a  band  edge  in  connection  with 
the  heavily  doped  material,  the  modifier  "effective"  will  be  used  to  avoid 
confusion.  Similarly,  the  energy  reference  in  this  discussion  will  be  the 
band  edge,  i.e.,  energies  in  the  (lightly  doped)  allowed  band  will  be  posi¬ 
tive  and  energies  in  the  (lightly  doped)  forbidden  band  will  be  negative. 

Both  of  these  conventions  will  be  followed  for  the  conduction  band  and  the 
valence  band. 
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heavily  doped  semiconductors.  However,  an  inspection  of  Eqs.  (8),  (7),  and  (11) 
shows  that,  as  the  impurity  concentration  decreases,  the  density  of  states  func¬ 
tion  approaches  the  intrinsic  or  lightly  doped  expression.  In  view  of  this  beha¬ 
vior,  we  assume  that  Bonch-Bruyevich's  equations  are  valid  at  impurity  concen¬ 
trations  below  those  at  which  the  impurity  band  has  overlapped  the  conduction  band. 

V.  The  Approach 

We  seek  a  description  of  the  density  of  quantum  states  applicable  for  all 
impurity  concentrations.  A  full  and  rigorous  quantum  mechanical  treatment  of  this 
problem  having  yet  to  be  done,  we  make  an  approximation  that  employs  knowledge 
available  to  date.  We  assume  the  dependence  of  the  density  of  states  to  be  given 
by  the  superposition  of  the  theories  due  to  Morgan  and  to  Bonch-Bruyevich.  Notice 
that  this  superposition  purports  to  hold  for  all  impurity  concentrations.  Upon 
this  point  rests  the  limit  of  validity  of  the  approximation,  for  the  individual 
theories  themselves  hold,  strictly  speaking,  only  for  restricted  ranges  of  concen¬ 
tration. 

To  explore  the  limits  of  validity,  consider  the  asymptotic  behavior.  For 
small  impurity  concentrations,  which  were  of  no  interest  to  Bonch-Bruyevich,  ob¬ 
serve  nonetheless  that  his  equations  reduce  to  the  conventional  density  of  states 
appropriate  to  a  lightly  doped  semiconductor.  So  also  does  Morgan's  equation, 
his  Gaussian  function  becoming  narrow  and  approximating  a  delta  function.  For 
very  large  concentrations,  Bonch-Bruyevich's  theory  is  valid  The  superposi¬ 

tion  of  Morgan's  and  Bonch-Bruyevich's  equations,  however,  fails  to  strictly 
account  for  the  states  in  the  impurity  band  having  been  removed  from  the  host  lat¬ 
tice  band.  This  deficiency  will  introduce  negligible  error  in  the  calculation  of 
the  total  density  of  states  provided  the  impurity  concentration  remains  small  com¬ 
pared  to  the  density  of  host  lattice  atoms.  The  same  restriction  on  the  impurity 
concentration  is  also  required  by  Morgan's  theory.  For  many  impurities  and  semi¬ 
conductors,  the  solubility  limit  of  impurities  will  ensure  that  this  restriction 


26 


is  not  violated.  Thus  the  superposition  at  large  densities  introduces  practically 
no  error  in  the  total  density  of  states  calculated.  Nonetheless  it  could  con¬ 
ceivably  introduce  sizeable  error  for  states  near  the  band  edge  because  it  is  here 
that  the  superposition  is  influenced  comparably  by  its  two  components.  As  will  be 
shown  in  the  following  section,  however,  the  superposition  predicts  a  density  of 
states  in  good  accord  with  an  empirically  determined  density  of  states  for  an 
impurity  concentration  exceeding  the  effective  density  of  states  (N^  in  the 
example).  For  silicon,  for  example,  this  is  nearly  the  maximum  Impurity  density 
seen  in  the  portions  of  a  device  that  dominate  its  performance. 

From  a  pragmatic  viewpoint,  then,  the  superposition  is  an  approximation  that 
one  may  expect  to  serve  adequately  over  the  entire  range  of  impurity  concentra¬ 
tion  that  is  of  technological  interest. 

VI.  The  Density  of  StAtes  in  p-type  Gallium  Arsenide 

Empirical  justification  for  superposing  the  theories  of  Morgan  and  Bonch- 
Bruyevich  follows  from  examination  of  the  recent  experimental  data  of  Mahan  and 
Conley Using  a  tunneling  theory,  Mahan  and  Conley  inferred  the  density  of 
states  in  p-type  gallium  arsenide  from  measurements  of  the  current-voltage  charac¬ 
teristics  in  gallium  arsenide-gold  Schottky  junctions.  They  report  measurements 

18  -3  18  -3 

made  with  samples  containing  zinc  impurity  densities  of  5.4x10  cm  and  9.9x10  cm  . 
A.  Qualitative  Discussion 

The  density  of  states  that  Mahan  and  Conley  infer  for  the  more  Impure  sample 
is  given  in  Fig.  1,  which  also  shows  our  resolution  of  their  data  into  probable  com¬ 
ponents  based  on  the  discussion  in  foregoing  sections.  The  left-hand  portion  of 
the  Gaussian-like  curve  was  obtained  by  simply  reflecting  the  right-hand  portion 
about  the  peak.  The  next  step  is  to  subtract  the  Gaussian-like  curve  from  the 
total  curve  to  obtain  the  third  or  remainder  curve.  Notice  that  this  resolution  of 
the  experimental  data  is  the  reverse  of  our  proposed  approach  of  superposing  the 

the  results  of  Morgan  and  Bonch-Bruyevich.  We  can  identify  the  Gaussian-like  curve 


27 


(oxeos  XaeuiTqae)  sa;e;s  J°  Atsuoq 


28 


Figure  1.  The  resolution  of  the  experimentally  determined  density  of 
states  in  p-type  GaAs.  (After  G.  D.  Mahan  and  J.  W.  Conley,  "The  Density 
of  States  In  Metal -Semi conductor  Tunneling,"  Applied  Physics  Letters,  11: 
29  (1967).) 


with  the  density  of  impurity  states  and  the  remainder  curve  as  the  density  of 
states  due  to  the  host  lattice  band  (valence  band  in  this  experimental  case). 
Notice  that  the  host  lattice  states  appear  as  though  they  can  be  described  by  a 
linear  energy  dependence  near  the  zero  of  energy  (the  band  edge)  and  by  a  para¬ 
bolic  function  beyond  the  band  edge.  This  observation  lends  support  to  our  pro¬ 
posed  use  of  Bonch-Bruyevich's  results. 

Mahan  and  Conley's  data  supports  our  hypothesis  that  the  non-degenerate  ion* 
izatlon  energy  remains  the  mid-point  about  which  the  impurity  band  develops  as 
the  impurity  concentration  Increases.  The  peak  of  the  Gaussian-like  curve 
occurs  at  -0.04584  ev  for  both  of  their  impurity  concentrations.  The  Gaussian- 
like  curve  of  Mahan  and  Conley  also  has  a  wider  spread  in  energy  at  the  higher 
impurity  concentration  than  at  the  lower. 

The  resolution  of  Mahan  and  Conley's  data  as  we  have  performed  it  is  further 

(18) 

supported  by  an  observation  of  Casey  and  Panishv  .  They  state  that  zinc  in 

18  -3 

gallium  arsenide  is  fully  ionized  for  concentrations  greater  than  5x10  cm  . 

Thus  we  can  conclude  that  the  ionization  energy  has  decreased  to  zero  as  a  result 

of  the  overlapping  of  the  impurity  and  host  lattice  bands.  Both  of  Mahan  and 

18  —3 

Conley's  samples  are  more  heavily  doped  than  5x10  cm  ;  therefore,  it  is  con¬ 
sistent  that  their  density  of  states  data  should  show  an  overlapping  of  the  two 
bands. 

The  assumption  that  the  mid-point  of  the  Gaussian-like  curve,  -0.04584  ev, 

is  the  lightly  doped  Impurity  ionization  energy  is  lent  further  validity  by 
(19) 

Huang  ,  who  places  the  ionization  energy  at  -0.0363  ev.  Since  his  material 
16  -3 

was  doped  at  5x10  cm  ,  we  expect,  in  view  of  Casey  and  Panish's  comment,  that 
the  ionization  energy  would  have  decreased  from  its  lightly  doped  value.  More¬ 
over,  Huang  states  that  hydrogenic  acceptor  impurities  are  expected  to  have 
ionization  energies  in  the  range  of  0.030  ev  to  0.050  ev  in  gallium  arsenide. 

Thus  the  support  for  our  assumption  concerning  the  mid-point. 
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B.  Quantitative  Discussion 

We  now  compare  the  quantitative  predictions  of  the  approach  to  Mahan 
and  Conley's  experimental  data.  Before  proceeding,  however,  let  us  rewrite 
into  convenient  forms  the  equations  that  we  will  be  using.  In  the  following 
equations,  energies  are  in  electronvolts,  concentrations  are  in  reciprocal 
cubic  centimeters,  and  the  density  of  states  is  in  reciprocal  electronvolts 
per  cubic  centimeter. 

The  density  of  impurity  states  is  given  by  rewriting  Eq.  (A)  as 


1  *1 

p(E)  -  3. 651x10®  e4  (—)4  N*5833  exp 

m 


-(E-Ejj)2  e2(m*/m)2 


9.551xl0-18  N*8334 


(13) 


In  Eq.  (13)  we  have  replaced  y  with  Ep  by  using  our  assumption  that  the  im¬ 
purity  band  will  develop  with  the  conventional  ionization  energy  as  its  mid¬ 
point.  The  density  of  states  near  E  ■  0  ev  is  given  by  Eq.  (10)  which  becomes 

P(E)  -  1.309X1017  <„*/„) 11/8  c3/S  N5/M  ♦ 


6. 267xl025  e3/8  (m*/m)13/8  N _5/24  E, 


(1A) 


when  the  various  universal  constants  are  used  to  evaluate  some  of  the  factors. 
Similarly,  Eq.  (8)  for  the  density  of  states  near  the  Fermi  level  becomes 


p (E)  -  6.8125xl021  (m*/m)3/2  E1/2  - 

A.065xl03  (m*/m)  e"3/2  N5/6  E-3/2  .  (15) 


The  first  term  of  Eq.  (15)  is  the  parabolic  density  of  states  function  used 
in  the  conventional  theory. 

There  are  sr:r.e  comments  that  need  to  be  made  about  the  ranges  of  validity 
of  Eqs .  (14)  and  (15).  For  the  impurity  concentrations  that  we  will  consider. 
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Eq.  (14)1*  valid,  according  to  Bonch-Bruyevich' a  lower  limit  -E^,  for  all 
energlea  at  which  p(E)  la  positive  and  therefore  meaningful.  Bonch-Bruyevich, 
however,  does  not  give  an  upper  limit  of  validity  for  Eq.  (14),  nor  does  he 
give  a  lower  limit  for  Eq.  (IS).  Because  of  this  we  have  adopted  the  fol¬ 
lowing  scheme.  Eqs.  (14)  and  (IS)  have  two  Intersections  at  energies  above 
* 

the  band  edge.  In  the  range  of  energy  between  these  intersections,  we 
apply  linear  weighting  functions  to  the  two  equations  and  use  the  sum  of  the 
weighted  equations  to  calculate  the  density  of  states.  The  weighting  func¬ 
tions  were  chosen  to  force  Eq.  (14)  to  be  the  only  contributor  to  the  sum  at 
the  lower  energy  Intersection  and  Eq.  (IS)  to  be  the  only  contributor  at 
the  higher  energy  Intersection. 

We  now  compare  the  predlctlona  of  these  equations  to  the  experimental 
data.  Flrat,  conalder  the  Gauss lan-1 Ike  curve  of  Fig.  1.  The  experimental 
data  can  be  analyzed  to  find  the  parameters  of  a  Gaussian  function  that  will 
then  be  an  approximate  mathcsuitlcal  description  of  the  data.  The  mean  of 
the  Gauaaian  la  0.04584  ev  as  discussed  in  Section  V.  The  standard  deviation, 
o,  can  be  calculated  from  the  data  by  taking  the  ratio  of  two  equations  of 
the  form  of  Eq.  (1),  one  evaluated  at  the  mean  and  one  evaluated  at  the  energy 
that  gives  the  function  the  value  of  one-half  Its  peak.  This  yields 

o  -  0.02799  ev  (16) 

for  Malian  and  Conley's  data  as  resolved  In  Fig.  1.  Morgan's  theoretical 
equation,  Eq.  (13)  yields 

o  -  0.03035  ev,  (17) 

when  the  material  parameters  of  GaAs  are  used  (m*  ■  0.7132  m^^\ 

(21) 

e  ■  12  ),  The  agreement  Is  satisfactory.  By  using  the  normalization 

*  The  other  interaectlons  do  not  occur  at  meaningful  energies. 
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property,  we  remove  the  arbitrariness  of  Mahan  and  Conley's  scale  of  ordi¬ 
nates.  The  peak  of  the  Gausslan-like  curve  is  then 

p(p)  ■  1.376x10^  ev”*  cm”"*  (18) 

Morgan's  equation  predicts  a  peak  value  of 

p(p)  ■  1.401x10^0  ev  *  cm  \  (19) 

which  agrees  satisfactorily  with  Eq.  (18). 

Figure  2  shows  that  Morgan's  equation  satisfactorily  matches  our  reso¬ 
lution  of  Mahan  and  Conley's  data.  The  dashed  Gaussian  curve  is  the  theoret¬ 
ical  prediction  of  Morgan's  equation,  while  the  solid  Gaussian-llke  curve  Is 
the  resolved  data. 

Figure  2  also  shows  a  (solid)  straight  line  which  has  been  visually 
fitted  to  the  resolved  curve  in  the  energy  range  near  E  »  0  ev.  The  equation 
for  this  straight  line  is 

p(E)  -  3.633xl020  +  9.911xl021  E,  (20) 

after  converting,  as  described  above,  from  Mahan  and  Conley's  arbitrary 
scale.  Bonch-Bruyevich's  equation,  Eq.  (14)  yields 

p(E)  -  2.936xl020  +  1.013xl022  E.  (21) 

The  agreement  is  again  satisfactory  and  within  the  errors  of  experiment 
and  graphical  analysis.  The  dashed  straight  line  in  Fig.  2  is  described  by 
Eq.  (21).  The  downward  shift  of  the  theoretical  line  with  respect  to  the 
fitted  line  is  approximately  19Z  at  E  ■  0  ev. 

VII.  Discussion 

Superposition  of  the  expressions  of  Morgan  and  Bonch-Bruyevich  yields 
a  description  of  the  density  of  states  in  a  semiconductor  as  a  function  of 


32 


f 


Che  impurity  concentration.  The  description  is  an  approximation  that  one 
may  expect  to  show  the  widest  discrepancy  with  reality  at  asymptotic  values 
of  concentration.  But  for  small  impurity  concentration,  the  superposition 
reduces  appropriately;  and,  for  large  concentrations,  the  limit  of  Impurity 
solubility  will  tend  to  protect  the  validity  of  the  approximation.  The  approach 
yielded  predictions  in  good  agreement  with  Mahan  and  Conley's  data  for  p-type 
gallium  arsenide. 

From  the  resultant  density  of  states  one  can  calculate  the  dependence 
on  impurity  concentration  of  various  parameters  that  interest  both  the  phys¬ 
icist  and  the  device  engineer:  the  Fermi  level,  transition  energies,  etc. 

Though  in  this  paper  we  have  explicitly  applied  the  approach  only  to  p-type 
gallium  arsenide,  it  applies  also  to  other  semiconductors.  Subsequent  papers 
will  deal  with  the  application  to  silicon  and  with  the  problem  of  Including 
the  complex  dependence  on  impurity  concentration  in  a  tractable  theory  of  pn 
junctions. 
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Appendix 


The  purpose  of  this  appendix  is  to  give  some  meaning  to  the  relative 
terms  "small,"  "large,"  etc.,  in  reference  to  impurity  concentration. 

A  small  impurity  concentration  is  one  for  which  the  impurity  atoms  have 
a  negligible  interaction  with  each  other.  Then  the  impurity  states  may  be 
treated  as  spatially  localized  states  in  which  an  electron  will  have  a  negli¬ 
gible  probability  of  tunneling  to  another  impurity  state.  At  small  concen¬ 
trations  the  energy  spectrum  of  the  impurity  states  is  adequately  described 
by  a  delta  function.  Electrical  conduction  takes  place  in  the  host  lattice 
bands . 

A  moderate  impurity  concentration  is  one  for  which  the  impurity  atoms 
have  a  sufficient  degree  of  interaction  among  themselves  that  there  is  a 
•mall  but  finite  probability  of  transfer  between  impurity  states.  At  temper¬ 
atures  near  0*K,  electrical  conduction  by  "hopping"  between  im¬ 
purity  states  will  be  detectable  in  this  concentration  range.  The  conductivity 
at  these  low  temperatures  will  be  enhanced  by  compensating  impurities  which 
would  have  the  effect  of  creating  empty  states  to  which  hopping  can  occur. 

At  the  upper  end  of  this  moderate  concentration  range  there  is  a  further 
possible  result  of  the  interactions  arising  because  of  a  difference  between 
the  usual  tight-binding  treatment  and  the  actual  circumstances  that  the  im¬ 
purities  experience.  (This  difference  is  in  addition  to  the  difference  in 
the  dielectric  medium  in  the  two  cases.)  The  usual  assumption  in  the  tight- 
binding  treatment  is  that  the  atoms  are  brought  together  and  formed  into  a 
perfect,  periodic  lattice.  This  is  not  the  case,  however,  for  impurities  in 
a  semiconductor.  The  impurities  assume  random,  substitutional  positions  in 
the  host  lattice,  and,  therefore,  any  given  atom  interacts  to  different  degrees 
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with  each  of  its  neighbor  impurity  atoms.  Hence,  along  spatial  paths  of 
locally  large  impurity  concentration,  there  will  be  splitting  of  the  impurity 
energy  levels  and  the  formation  of  threadlike  energy  bands  through  the  crystal. 

When  the  impurity  atoms  interact  sufficiently  with  each  other  that  the 
impurity  states  are  no  longer  localized  and  an  impurity  band  exists  quite 
generally  throughout  the  crystal,  the  impurity  concentration  is  described  as 
large.  In  this  concentration  range  the  conduction  band  is  noticeably  dis¬ 
torted  by  the  formation  of  band  tails.  Electrical  conduction  at  low  tempera¬ 
tures  occurs  in  the  Impurity  band  and  is  degraded  by  compensating  impurities 
which  act,  in  this  case,  to  remove  mobile  carriers  from  the  impurity  band 
instead  of  providing  more  available  states  as  compensation  does  in  moderately 
doped  material. 

A  very  large  impurity  concentration  is  characterized  by  the  merging  of 
the  conduction  band  states  and  the  impurity  band  states.  The  impurity  ion¬ 
ization  energy  has  become  zero.  Electrical  conduction  in  this  case  has  prop¬ 
erties  which  are  nearly  temperature  insensitive  and  similar  to  the  type  of 
behavior  found  in  metals. 


*  If  the  impurities  did  in  fact  occur  in  a  perfectly  periodic  pattern  in 
the  lattice,  they  would  not  be  able  to  cause  the  scattering  events  which 
manifest  themselves  as  the  impurity-determined  components  of  mobility:  the 
ionized  impurity  and  neutral  impurity  mobilities (15), 
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Symbol  List 


The  number  in  parentheses  following  the  definition  indicates  the  equa¬ 
tion  in  which  the  symbol  is  first  found.  Symbols  bearing  the  subscript  "o" 
designate  quantities  evaluated  in  thermal  equilibrium  unless  otherwise  speci¬ 
fied  in  the  following  list.  The  cgs-esu  system  is  used  with  the  exception 
that  energy  is  expressed  in  electronvolts . 

a  defined  in  Eq.  (7) 

aQ  the  Bohr  radius  in  a  crystal  (12) 

E  general  energy  (1) 

E^  defined  in  Eq.  (11) 

ED  the  energy  of  the  donor  level  in  the  forbidden  gap  (13) 

e  the  magnitude  of  the  electronic  charge  (2) 

h  Planck’s  constant 

t  h/  (2 tt )  (3) 

I  Bessel  function  with  imaginary  argument  (5) 

m  free  electron  rest  mass  (4) 

m*  the  density  of  states  effective  mass  (3) 

N  impurity  concentration  (1) 

rQ  a  parameter  in  Morgan's  theory  (See  text  below  Eq.  (4)) 

x  a  normalized  energy  defined  in  Eq.  (6) 

e  the  relative  dielectric  constant  (2) 

y  the  mean  of  the  Gaussian  function  (1) 

p  the  density  of  quantum  states  (1) 

pQ  the  parabolic  density  of  quantum  states  used  in  conventional  theory  (8) 

o  the  standard  deviation  of  the  Gaussian  function  (1) 

A.  screening  length  (2) 
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III.  Insulating  and  Semiconducting  Glasses  (E.  R.  Chenette,  R.  W.  Gould, 

L.  L.  Hench,  J.  J.  Hren) 

A.  DIELECTRIC  BREAKDOWN  OF  CERAMICS  (G.  C.  Walther  and  L.  L.  Hench) 
Introduction 

Ceramics  are  widely  used  in  the  electrical  industry  for  insulation, 
capacitors,  and  encapsulation.  In  each  of  these  applications  the  ceramic 
material  is  exposed  to  a  voltage  gradient  and  it  must  withstand  the  gradient 
for  the  operating  life  of  the  system.  Failure  occurs  when  an  electrical 
short  develops  across  the  material  and  such  a  failure  is  called  dielectric 
breakdown.  The  voltage  gradient,  expresses  as  volts/cm,  sufficient  to  pro¬ 
duce  the  short  is  termed  the  breakdown  strength  of  the  material. 

Breakdown  strengths  of  ceramics  vary  widely  due  to  many  facotrs.  Some 
of  the  most  important  include:  thickness,  temperature,  ambient  atmosphere, 
electrode  shape  and  composition,  surface  finish,  field  frequency  and  wave 
form,  porosity,  crystalline  anisotropy,  amorphous  structure,  and  composition. 

A  detailed  characterization  of  the  importance  of  these  factors  on  breakdown 
strengths  of  a  variety  of  ceramics  has  not  been  accomplished.  A  general 
comparison  of  several  variables  can  be  made  from  the  data  accumulated  in 
Table  I. 

Strengths  as  large  as  several  million  volts/cm  are  reported  for  micron 
thick  thin  films.  Because  they  are  so  thin  the  voltage  insulation  of  such 
films  is  small,  however.  Thus,  bulk  ceramics  of  many  centimeters  in  thickness 
are  required  to  insulate  large  electrical  power  voltages.  Table  I  shows  that 
the  breakdown  strength  unfortunately  decreases  drastically  to  levels  of  only 
several  thousand  volts/cm  for  bulk  ceramics. 

The  marked  changes  in  strength  occurring  with  changing  thickness  are  due 
to  changes  in  breakdown  mechanisms.  The  temperature  effects  on  breakdown 
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shown  in  Table  1  are  due  to  the  strong  influence  of  thermal  energy  on  break¬ 
down  mechanisms.  It  is  the  objective  of  this  paper  to  review  the  various 
theories  of  breakdown  mectianisms  and  discuss  the  influence  of  experimental 
and  materials  variables  on  the  magnitude  of  breakdown  strengths  obtained  for 


ceramics. 


Definitions 


When  the  temperature  (i.e.,  thermal  energy)  of  the 
lattice  or  its  electrons  reaches  a  value  during  the  appli¬ 
cation  of  an  electric  field  such  that  the  conductivity 
increases  rapidly  and  results  in  permanent  damage  to  the 
material  dielectric  breakdown  occurs.  There  are  three  basic 
types  of  breakdown  called  intrinsic,  thermal  and  avalanche 
breakdown.1  Historically,  there  has  also  been  three  pseudo¬ 
types  of  breakdown  described.  The  psuedo- types  are  termed 

2 

discharge,  electrochemical  and  mechanical  breakdown. 

Psuedo- types  can  be  considered  to  be  produced  by  one  or 
more  of  the  three  basic  mechanisms.  Each  of  the  types  of 
failure  will  be  briefly  defined  and  then  the  fundamental 
breakdown  mechanisms  will  be  discussed  in  detail. 

Dielectric  discharge  is  associated  with  a  gaseous  break¬ 
down  in  the  pores  or  at  the  surface  of  a  solid  material. 
Electrochemical  breakdown  is  a  result  of  a  gradual  deteri¬ 
oration  of  insulating  properties  through  chemical  reactions 
until  breakdown  occurs  by  one  of  the  basic  mechanisms. 
Mechanical  breakdown  is  due  to  cracks,  defects  and  other 
stress  raisers  distorting  the  applied  field  and  thus  pre¬ 
cipitating  failure. 

Experimentally,  intrinsic  breakdown  is  found  to  be 
primarily  field  dependent  in  that  the  applied  field  deter¬ 
mines  when  the  electron  temperature  reaches  the  critical 
level  for  breakdown.  Observations  that  intrinsic  breakdown 
occurs  at  or  below  room  temperature  and  occurs  in  very  short 
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time  intervals,  approximately  a  microsecond  or  less,  is 
strong  evidence  that  it  is  electronic  in  nature.  The  name 
"intrinsic"  is  used  because  breakdown  by  this  mechanism  is 
independent  of  the  sample  or  electrode  geometry  used  (pro¬ 
vided  no  field  distortion  occurs)  or  of  the  wave  form 
applied.  Hence  the  value  of  applied  field  required  to  cause 
intrinsic  breakdown  at  a  given  temperature  is  a  property 
solely  dependent  on  the  material. 

Breakdown  observed  between  room  temperature  and  approxi¬ 
mately  300°C  is  described  as  thermal  breakdown,  because 
reaching  the  critical  thermal  energy  for  breakdown  is  pri¬ 
marily  influenced  by  the  ambient  temperature,  and  not  the 
electric  field.  Also,  thermal  breakdown  is  dependent  on 
the  rate  of  application  of  the  field.  Slow  increases  in 
field  cause  breakdown  in  milliseconds  to  minutes,  with  the 
value  of  breakdown  being  influenced  by  sample  geometry. 

For  faster  field  pulses,  breakdown  is  independent  of  geometry 
and  breakdown  strength  increases  with  shorter  pulse  times. 

Avalanche  breakdown  is  related  to  intrinsic  breakdown 
in  that  it  occurs  at  relatively  low  temperatures  and  short 
times.  However,  thermal  properties  of  the  material  are 
used  to  describe  the  breakdown  behavior,  so  it  is  properly 
a  combination  of  thermal  and  intrinsic  mechanisms.  Thin 
samples,  such  as  dielectric  films,  undergo  avalanche 
destruction.  There  is  a  statistical  variation  in  breakdown 
strength  with  short  pulses  for  this  mechanism  of  breakdown. 


Pre-breakdown  noise  is  found  with  slower  pulses  which  indi¬ 
cates  a  sequential  type  of  lattice  destruction  is  occurring. 


Breakdown  Mechanisms 


The  basic  theoretical  approach  used  to  describe  dielectric 
breakdown  in  solids  is  to  develop  an  energy  balance  equation 
of  the  form 


A(T0,F,a)  -  B(T0,a)  (1) 

where  A(T0>F,a)  ■  energy  gained  by  the  material  from  the 

applied  field 

B(To,a)  *  energy  dissipated  by  the  material 

Tq  ■  lattice  temperature 

F  *  applied  field 

a  *  energy  distribution  parameter,  which 

depends  on  the  model  proposed 
Thus  A  *  B  is  the  limiting  condition  for  breakdown. 

Intrinsic  Breakdown 

Theories  of  intrinsic  breakdown  can  be  classified  as 
those  dealing  primarily  with  the  electron- lattice  energy 
transfer  and  those  considering  changes  in  the  electron  energy 
distribution  in  the  material. 

In  electron- lattice  interaction  models  the  behavior 
of  the  material  is  approximated  by  considering  a  single 
electron  of  average  E  (i.e.,  a  =  E) .  Consequently,  there 
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is  a  low  electron  density,  a  small  probability  of  electron- 
electron  interaction,  and  only  electron- lattice  energy 
transfer  is  allowed.  The  mechanisms  contributing  to  this 
transfer  are: 

(1)  Lattice  vibrations  in  a  dipolar  field 

(2)  Electron  shell  distortion  accompanying  the  dipolar 
field- lattice  vibrations 

(3)  Short  range  electron  shell  distortion  in  a  non¬ 
polar  field 

The  mechanism  chosen  to  describe  the  electron-lattice  energy 
transfer  depends  on  the  model  of  the  material  used.  The 
problem  is  to  determine  the  value  of  E  required  for  the 
energy  balance  and  then  to  calculate  F  ,  the  critical  field 
strength,  that  causes  the  critical  temperature  for  break¬ 
down,  T  ,  to  be  reached, 
c 

One  model,  called  the  Frohlich  high  energy  criteria,3 
assumes  that  multiplying  the  number  of  conduction  electrons 
to  a  very  large  value  will  destroy  the  lattice.  Thus, 

E  =  I,  the  ionization  energy  or  gap  between  the  valence 
and  conduction  band.  Breakdown  occurs  when  the  ionization 
rate  exceeds  the  recombination  rate  causing  the  conduction 
electron  density  to  increase  irreversibly.  Figure  1  shows 
the  energy  balance  equation  schematically.  The  value  of  E 
at  the  intersection  of  the  A  and  B  curves  indicates  the 
energy  an  electron  must  have  tc  be  accelerated  by  the 
appropriate  field,  F.  The  high  energy  electron  resulting 
from  a  recombination  collision  will  have  E  <  I  so  that  only 
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those  fields  F  >  F  can  cause  the  electron  to  be  accelerated. 
—  c 

resulting  in  another  ionization  and  thus  increasing  the 

density  of  conduction  electrons. 

The  above  model  has  been  questioned1  on  the  following 

grounds:  (1)  An  ionization-recombination  steady  state  may 

not  be  realizable  in  the  short  time  before  breakdown,  and 

(2)  a  steady  state  might  not  occur  physically:  the  energy 

balance  might  be  maintained  by  some  non- steady  state  process 

such  as  radiative  recombination  (i.e.,  light  emitted  upon 

recombination,  as  observed,  but  not  published,  by  Kelly).1 

Another  model  is  the  low  energy  criteria  due  to  von 
4  S 

Hippel  and  Callen.  They  assumed  the  limiting  condition 

for  breakdown  to  occur  when  B(E',TQ)  is  a  maximum,  with  E' 

being  the  maximum  average  electron  energy  required  for  steady 

state  conditions.  This  means  that  F'  is  the  field  needed  to 

c 

accelerate  all  the  conduction  electrons  against  the  lattice 
influence  and  is  also  shown  on  Figure  1.  This  model  faces 
criticism  similar  to  the  above  in  addition  to  the  improb¬ 
ability  of  achieving  the  extremely  high  field,  F^.,  necessary 
to  accelerate  all  the  conduction  electrons  simultaneously. 

Stratton**  attempted  to  extend  the  above  ideas  to  non¬ 
polar  materials  but  agreement  between  theory  and  experiment 
is  poor.1  It  should  be  mentioned  that  in  general  the  actual 
breakdown  field,  Ffi,  must  be  greater  than  F£,  the  latter 
being  only  the  value  necessary  to  cause  an  irreversible 
imbalance  in  the  energy  transfer  equation. 
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Energy  E 

FIGURE  1 

Schematic  representation  of  the  energy  balance  relation 
for  the  Frohlich  High  Energy  Criteria,  A(F  ,T  E)  =  B(E,T  ), 
and  the  Von  Hippel-Callen  Low  Energy  Criteria?  ° 

A(F*  »T  »E)  =  B(E\T  ). 


Intrinsic  theories  based  on  changes  in  the  electron 
energy  distribution  of  the  solid  have  also  been  advanced. 
These  theories  are  based  on  five  contributing  factors: 

(1)  Acceleration  of  electrons  due  to  the  applied  field 

(2)  Collisions  between  conduction  electrons 

(3)  Collisions  between  conduction  electrons  and  the 
lattice 

(4)  Ionization  or  recombination  of  electrons  to  and 
from  the  valence  band  or  traps 

(5)  Diffusion  due  to  a  field  gradient 

This  large  number  of  theoretical  variables  can  be  reduced 
by  ignoring  ionization-recombination  processes.  This  is 
possible  because  it  can  be  assumed  that  most  electrons  in 
dielectrics  occupy  energy  states  intermediate  between  those 
low  enough  to  allow  recombination  or  high  enough  to  permit 
ionization.  The  diffusion  mechanism  can  also  be  dismissed 
by  not  considering  the  influence  of  the  field  source. 

The  energy  distribution  in  the  solid  can  be  classified 
under  assumptions  of  either  a  low  or  high  density  of  con¬ 
duction  electrons.  Assuming  a  low  electron  density  permits 
ignoring  electron-electron  collisions.  For  low  energy 
electrons  a  Boltzmann  distribution  is  obtained.  However, 

for  high  energy  electrons,  Frohlich  showed  that  a  valid  dis- 

7 

tribution  function  could  not  be  described.  The  distribu¬ 
tion  approaches  infinity  with  increasing  energy  unless 
ionization-recombination  properly  makes  the  theory  an 
avalanche  type. 
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For  high  electron  density  theories  it  is  assumed  that 

electron-electron  collisions  are  so  important  that  they 

determine  the  electron  energy  distribution.  Electron- 

lattice  interactions  determine  the  critical  value  of  electron 

temperature  Tc,  usually  greater  than  the  lattice  temperature 

7  -  x 

T  .  Critical  conduction  electron  densities  of  n„  •  10  cm 
o  c 

„  14  .3 

and  nc  -  10  cm  have  been  calculated  for  polar  and  non¬ 
polar  materials,  respectively.  Thus  for  these  models, 
a  »  Tc  and  the  breakdown  strength  is  calculated  by  deter¬ 
mining  the  field,  F^,  necessary  to  make  T  increase  to 
infinity. 

An  electron  density  distribution  model  for  a  pure  crys- 

o 

talline  solid  has  been  proposed  by  Frohlich  and  Paranjape. 
With  few  or  no  defects  there  would  be  few  traps.  Conse¬ 
quently,  in  the  energy  dissipation  function,  B,  ionizing 
collisions  were  also  assumed  to  be  of  minor  significance 

and  therefore  is  one  point  of  criticism.  Another  less  than 

17  -3 

satisfactory  aspect  of  the  model  is  that  for  n  ■  10  cm 
and  F  <  F  ,  the  energy  gain  function  A  is  a  very  large  value 
of  600y  electron  volts  per'  ion  volume,  where  y  is  the 
mobility.  For  normal  values  of  y  this  magnitude  of  energy 
absorption  would  destroy  the  crystal  at  less  than  critical 
fields . 

Frohlich  has  also  proposed  a  model  for  amorphous 
materials  in  which  he  assumes  a  trap  distribution  as  shown 
in  Figure  2.  These  are  shallow  S  traps  just  below  the 
conduction  band  and  deep  D  traps  within  the  gap.  The 
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FIGURE  2 


Trap  level  representation  for  the  Frohlich  amorphous 
material  model. 
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density  of  electrons  in  the  traps  is  n  <  n  <<  nn  and  n  = 

C  S  il  c 

4  M  •» 

101  cm  nc  plus  ns  thus  determines  the  conduction 
electron  energy  distribution  and  the  electron  temperature  T. 
The  main  mechanism  of  energy  transfer  in  this  model  is 
emission  of  phonons  associated  with  the  change  of  electrons 
between  S  levels  and  the  conduction  band.  O'Dywer  has  ex¬ 
tended  this  theory  to  include  isolated  defect  levels  in 
9 

crystals. 

Criticism  of  the  amorphous  model  is  centered  on  the 
strong  possibility  of  emission  from  traps,  which  would  change 
the  electron  distribution.  Figure  2  is  also  a  very  ideal¬ 
ized  distribution  of  energy  levels  for  an  amorphous  solid. 

The  current  state  of  intrinsic  breakdown  theory  thus 
appears  to  be  that  the  actual  intrinsic  breakdown  mechanism 
or  mechanisms  may  be  more  complex  than  those  proposed  by 
present  breakdown  models.  Further  experimental  investiga¬ 
tion  to  reduce  experimental  difficulties  may  resolve  some 
of  the  complexity  or  provide  insight  into  modifications  of 
present  theory. 


thermal  Breakdown 

Thermal  breakdown  theory  also  is  based  on  an  energy 
balance  relation,  but  it  is  the  balance  between  heat  dis¬ 
sipated  by  the  sample  and  the  heat  generated  due  to  Joule 
heating,  dielectric  losses,  and  discharges  in  the  ambient. 
Hence,  it  is  the  lattice  temperature  and  not  an  electron 
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temperature  that  must  reach  a  critical  level  for  breakdown 
to  occur.  The  influence  of  the  applied  field  is  only  in¬ 
directly  felt  as  it  influences  the  heat  generating  mechanisms 
and  does  not  play  the  determining  role  evident  in  the  in¬ 
trinsic  theories  based  on  the  electron  temperature  in  the 
solid.  Because  of  the  relatively  weak  dependence  between 
field  and  temperature,  the  value  Tc  is  not  too  important  and 
the  actual  lattice  temperature  at  breakdown,  T^,  is  usually 
somewhat  greater. 

The  basic  relation  for  thermal  breakdown  is 

Cy  jjl  -  div (K  grad  T)  =  c(F,To)F2  (2) 

where  Cy  *  heat  capacity  of  the  material 

div(K  grad  T)  *  heat  conduction  of  a  volume  element 
2 

o(F,To)F  =  heat  generation  term 
There  is  no  charge  accumulation  so  the  current  is  contin¬ 
uous  . 

For  slow  field  applications  a  steady  state  is  assumed 

7 

and  the  Cv  dT/dt  term  can  be  ignored.  Calculations  employ¬ 
ing  this  model  show  that  breakdown  strength  is  inversely 
proportional  to  the  square  root  of  the  thickness.  This 
square  root  dependence  agrees  well  for  thin  samples  where 
a  uniform  temperature  can  be  realized  but  for  thicker 
samples  the  breakdown  strength  is  experimentally  observed 
to  be  inversely  proportional  to  the  thickness  itself.  The 
criteria  for  determining  if  a  sample  is  thick  or  thin 
depends  on  several  materials  constants,  namely  the  thermal 
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conductivity  as  well  as  the  pre-exponential,  oq,  and  the 
activation  energy,  Q,  for  the  conductivity- temperature 
dependence  of  the  solid,  i.e.,  o  *  cQ  exp(<j>/kT). 

When  the  field  pulse  is  fast  enough,  negligible  heat 
transfer  can  take  place.  Consequently,  the  heat  conduction 
term  can  be  ignored  and  the  electrodes  only  influence  the 
field  distribution  and  not  the  heat  flow,  as  in  the  previous 
case.  By  integration  it  is  possible  to  calculate  t  ,  the 

V 

time  for  breakdown  to  occur  after  reaching  T  : 

V* 


t 


c 


CydT 


°CFc’W 


(3) 


The  strong  dependence  of  t  on  F  is  easily  seen  in  the 
above  expression.  Similarly,  it  can  be  shown  that  F  ■ 

KT 

^ 2'  exp (<J>/2kT)  ,  thus,  the  critical  breakdown  field  is 
essentially  independent  of  Tc. 

Numerical  solutions  to  thermal  breakdown  models  must  be 
used  when  the  field  frequency  is  intermediate  and  either 
of  the  above  simplifications  can  not  be  made.  Each  case 
becomes  a  different  problem  with  different  boundary  con¬ 
ditions.  Numerical  techniques  can  also  be  extended  to 
cover  multi-dimensional  samples  or  different  waveforms  with 
some  success.  There  are  no  serious  criticisms  of  thermal 
breakdown  theory. ^ 
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Avalanche  Breakdown 


While  thermal  breakdown  theory  may  be  more  practical 
for  the  application  temperatures  of  many  ceramic  materials, 
if  the  material  geometry  approaches  that  of  a  thin  film, 
avalanche  theory  may  prove  more  useful.  Avalanche  break¬ 
down  theory  is  an  attempt  to  combine  features  of  intrinsic 
and  thermal  theories  since  an  electron  distribution  insta¬ 
bility  will  have  thermal  consequences.  The  electron  be¬ 
havior  is  described  with  an  intrinsic  theory  and  the  break¬ 
down  criteria  is  based  on  thermal  properties.  Avalanche 
theory  considers  the  gradual  or  sequential  buildup  of 
charge  rather  than  the  sudden  change  in  conductivity,  even 
though  the  charge  buildup  may  occur  in  a  very  short  time. 

Avalanche  theories  can  be  easily  classed  as  to  their 
initiation  mechanisms  which  are  either  field  emission  or 
ionization  collision.  Field  emission  assumes  that  the 
conduction  electron  density  increases  by  tunneling  from  the 
valence  band  to  traps  or  to  the  conduction  band.10  Consid¬ 
ering  that  the  probability  for  emission  is 

P  *  aF  exp(-bI2/F)  (4) 

shows  that  P  is  small  until  F  is  quite  high  (a  and  b  are 
constants) .  Using  this  development  gives  an  order  of  mag- 
nitude  calculation  of  F  =  10  volts/cm,  with  the  impulse 
thermal  criteria  T  *  Tc  being  the  critical  parameter. 
Criticism  of  initiation  by  field  emission  may  be  made  on 
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at  least  t h r . •<  points:  (1)  there  should  be  field  emission 
only  for  gaps  °f  '  ev  or  less,  i.e.,  for  semiconductors, 
or  there  (21  may  be  emission  from  traps,  and  (3)  thermal 
or  phonon  assisted  emission  is  not  included. 

The  simpler  single  electron  approach  to  an  ionization 

collision  theory  seems  to  be  preferred  to  the  more  complex 

many  electron  or  avalanche  multiplication  view.  The  single 

12 

electron  model  assumes  that  at  least  10  conduction  elec¬ 
trons  per  cm3  are  needed  to  disrupt  the  lattice.  'Vhen 
one  starting  electron- ionizing  collision  liberates  two 
electrons,  which  in  turn  liberate  four  and  so  on,  it  will 
take  about  40  such  collisions  to  achieve  breakdown.  This 
simple  approach  gives  a  critical  field  F  that  is  dependent 

v 

on  thickness:  the  sample  must  be  thi  k  enough  to  have  at 
least  40  mean  free  path  lengths.  However,  the  sample  can 
not  be  so  thick  that  the  conduction  electron  density 
becomes  so  large  that  electron-electron  collisions  occur, 
thus  limiting  the  electron  energy  to  that  below  the  ioni¬ 
zation  1  eve  1 . 

The  many  electron  avalanche  theory  gives  a  relation 
between  the  electron  energy  distribution  function  and  the 
ionization-recombination  rates.  There  are  two  views  of 

12 

this  approach  to  dielectric  breakdown.  One,  due  to  Heller, 
assumes  that  the  ionization  rate  is  so  large  that  recom¬ 
bination  can  not  keep  up.  It  is  assumed  in  the  other 
model  by  Franz1  and  Veelke,  that  recombination  can  be 
dismissed  since  ionization  increases  the  number  of  conduction 
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election  so  rapidly  that  no  recombination  mechanism  can 
quench  the  process.  The  actual  case  is  likely  to  be  a  com¬ 
bination  oi  the  two.  The  electron  source  is  ignored  because 
the  density  of  the  source  electrons  has  no  relation  to  F  . 

Criticism  on  several  points  makes  the  simpler  single 
electron  theory  appear  more  realistic  than  the  ones  in¬ 
volving  multiple  processes.  The  mu 1 tip le  electron  theory 
is  limited  to  short  path  lengths  because  breakdown  should 
occur  before  electron  interactions  become  important.  The 
imbalance  between  ionization  and  recombination  is  due  to 
rapid  saturation  of  traps,  which  is  always  unstable,  or  to 
some  emptying  mechanism,  which  is  unstable  only  if  traps 
empty  faster  than  they  fill.  Traps  might  fill  with  pre- 
breakdown  avalanches  having  slow  decay  times  but  phospho¬ 
rescence  studies  show  the  high  instability  of  traps,  which 
is  important  for  this  theory.  Also,  the  electron  energy 
distribution  occurs  only  from  electron- lattice  collision 

so  that  the  electron  density  is  lower  than  the  Frohlich- 

7 

Paranjape  collective  avalanche  proposal.  Finally,  the 
effects  of  the  cathode  on  such  a  short  path  of  ionization 
avalanche  is  ignored. 

O'Dywer  has  recently  proposed  a  space  charge  modified 
field  emission  model*'’  for  avalanche  breakdown.  Ilis 
criticism  of  the  single  electron  or  40  generation  model 
is  based  or  its  lack  of  a  continuity  of  current.  Assuming 
that  the  freed  electrons  and  holes  occurring  after  the 
fortieth  collision  form  a  parallel  plate  capacitor, 
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calculations  show  that  a  field  of  10*1  v/cm  would  be  needed 
to  maintain  this  charge  separation.  Since  such  a  large 
breakdown  strength  is  not  observed,  he  suggests  a  continuous 
current  model  involving  cold  cathode  field  emission  with 
collision  ionization  occurring  after  emission,  for  short 
breakdown  times.  Only  electrons  of  energy  kT  (i.e.,  room 
temperature)  are  assumed  to  be  accelerated  by  the  field  to 
the  ionization  level,  although  intermediate  electron 
energies  could  cause  ionization  or  some  electrons  with 
sufficient  energy  might  not.  Calculation  of  a  field  versus 
distance  relation*5  indicate  that  a  space  charge  of  low 
mobility  holes  is  influencing  the  field  distribution  and 
hence  breakdown  strengths,  as  seen  by  the  drop  in  field 
near  the  cathode,  shown  in  Figure  3.  Comparison  of  the 
theoretical  results  with  A^Oj  and  NaCl  thin  film  breakdown 
data  gives  good  agreement.  The  influence  of  space  charge 
on  the  temperature  dependence  at  higher  temperatures  is 
not  known. 

This  approach  combines  features  of  field  emission  and 
ionization  collision  which  were  mutually  excluded  in  previous 
proposals.  It  is  therefore  potentially  subject  to  the  same 
criticism  as  the  other  theories  although  it  does  account 
for  thickness  dependence  and  cathode  effects.  Several 
experimenters  who  were  inclined  to  interpret  their  results 

i  4  i  n 

in  terms  of  field  emission  or  ionization  collision  * 

before  this  theory  appeared  are  now  willing  to  consider  the 

18  19 

effect  of  space  charge  ’  or  at  least  of  the  possibility 
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of  both  mechanisms.  Others  *  also  feel  that  space 

charge  is  worth  considering  as  a  possible  mechanism.  In 

addition,  it  appears  Mott  feels  such  a  mechanism  is  probably 

responsible  for  switching  effects  in  some  amorphous  semi- 
23 

conductors . 


Experimental  Variables 

The  great  number  of  experimental  and  materials  variables 
that  can  influence  dielectric  breakdown  accounts  for  the 
wide  variance  between  the  results  of  different  investigators 
when  studying  the  same  material.  Comparison  of  breakdown 
strengths  between  different  materials  is  especially  diffi¬ 
cult  because  of  failures  to  report  materials  characteriza¬ 
tions  in  the  literature.  The  following  description  of  the 
importance  of  experimental  and  materials  factors  is  by  no 
means  exhaustive  but  should  indicate  the  large  amount  of 
work  that  is  still  to  be  done  in  breakdown  studies. 

Sample  geometry  influences  field  distribution  and  heat 
dissipation  considerably.  The  geometric  parameter  of 
major  study  has  been  the  thickness  of  the  breakdown  path 
on  the  breakdown  strength.  Figure  4  shows  a  general  result 
for  thin  films'^**4  and  indicates  that  the  breakdown 
strength  decreases  sharply  with  distance  at  some  critical 

value.  The  same  general  pattern  is  observed  for  thicker 
2  5 

specimens,  '  although  some  doubt  has  been  raised  with  regard 

2  6 

to  amorphous  materials. 
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Breakdown  Strength  (lO  V/cm  ) 


FIGURE  4 

Thickness  dependence  of  breakdown  strength  typical  of 
thin  films. 
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A  curve  showing  the  general  temperature  dependence  of 

many  materials  is  shown  in  Figure  5.  In  the  low  temperature 

region  the  breakdown  strength  increases  with  increasing 

temperature  and  then  declines  with  further  increases.  The 

location  of  the  maxima  depends  on  the  material  and  large 

deviations  from  this  pattern  have  been  reported,  including 

21 

temperature  independent  regions.  The  influence  of  thermal 
history  appears  to  have  been  generally  ignored  as  a  varia¬ 
ble,  although  the  affect  of  annealing  on  some  glasses  has 
been  reported  to  have  little  effect  on  breakdown  strength. 
The  breakdown  strength  of  measured  samples  generally 

increases  when  the  surrounding  environment,  such  as  special 

2 

oils,  has  a  high  dielectric  strength.  It  is  well  known 
that  the  ambient  atmosphere  can  alter  the  stoichiometry  of 
many  ceramic  materials  and  that  changes  in  stoichiometry 
influence  electrical  conductivity.  However,  a  detailed 
investigation  of  stoichiometry  effects  on  breakdown 
strength  is  not  known  to  the  authors.  It  might  be  pre¬ 
dicted  that  breakdown  strength  will  decrease  when  the 
stoichiometry  changes  such  that  the  conductivity  increases. 

The  applied  field  itself  is  an  important  parameter, 
in  that  not  only  its  magnitude,  but  its  frequency,  duration, 
waveform,  and  distribution  can  also  be  varied.  This  can 
be  very  important,  especially  for  intrinsic  type  breakdown, 
which  is  directly  field  dependent.  In  general,  it  has 

been  observed  that  dc  fields  produce  lower  breakdown 

27 

strengths  than  ac  fields  for  thin  specimens.  The  reverse 
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Breakdown  Strength 


Nominal  temperature  dependence  of  breakdown  strength 
for  many  materials. 


is  true  for  thick  specimens  due  to  dielectric  heating  and 
poor  heat  dissipation.  The  greatest  discrepancies  to  this 
generalization  occur  at  very  high  and  very  low  frequencies. 

28 

Many  researchers  have  used  squarewave  pulses,  either  singly 

or  in  series,  to  determine  intrinsic  breakdown  and  they 

observe  that  increased  breakdown  values  are  obtained  with 

shorter  pulses  and  faster  rise  times.  Prestressing,  or 

applying  a  field  less  than  breakdown  for  some  time,  has 

received  too  little  attention,  especially  since  the  influence 

of  space  charge  on  breakdown  strength  is  being  considered 

more  favorably.  Experiments  on  glass  indicate  a  very 

22 

definite  influence  of  prestressing  and  it  has  been  ob¬ 
served  for  some  time  that  there  is  a  threshold  field  such 
that  application  of  a  lesser  field  will  not  cause  breakdown 
in  any  time  interval.2 

The  electrodes  used  for  determining  breakdown  strength 

have  received  much  attention  from  investigators.  Vacuum 

depositions  of  aluminum  or  gold  appear  to  be  the  most 

reliable  electrode  metal"^  while  electrolytes  or  graphite 

17 

may  be  suitable  in  certain  cases.  Even  the  polarity  can 

make  a  difference  if  the  electrode  geometry  does  not  pro- 

19 

vide  a  uniform  field  intensity,  especially  in  light  of 
space  charge  considerations.  The  observation  of  micro- 

o  IQ 

projections  40-500  A  high  on  thin  film  electrodes  and 
evidence  of  current  asymmetry  due  to  such  projections"^ 
lends  further  support  to  theories  based  on  field  emission. 
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Electrode  configuration  has  been  a  subject  of  concern 
for  many  years.  However,  for  thick  samples  a  cup  and  plate 
electrode,  shown  in  Figure  6,  approximating  a  sphere  and 
plate  configuration,  appears  to  give  fairly  reliable 
results.  Generally,  the  values  of  t  and  r  are  about 
0.1mm.  The  ratio  d/t  must  be  varied  with  temperature  to 
obtain  reliable  results:  It  should  be  about  4/1  at  room 
temperature  but  increased  to  12/1  at  -195°C,  to  avoid  break¬ 
down  at  the  edges  of  the  sample.  For  amorphous  materials 

even  thicker  or  diffusely  deposited  electrode  edges  may  not 

32 

prevent  breakdown  over  undesired  paths.  The  success  with 
electrolytes  and  the  effect  of  projections  shows  the  im¬ 
portance  of  having  good  contact  between  the  electrode  and 
29 

the  sample.  Even  the  electrode  area  can  influence  results, 
as  shown  by  a  study  of  the  anisotropy  of  breakdown  in 
single  crystals. ^ 

The  most  important  factor  to  characterize  accurately, 
the  dielectric  material  being  investigated,  has  had  very 
little  systematic  investigation.  In  glasses,  it  is  gener¬ 
ally  found  that  the  effect  of  composition  on  breakdown 
strength  follows  the  change  in  resistivity  with  additions 
of  alkali  ions.  *  Less  is  known  about  compositional 
effects  in  ceramics,  either  as  impurities  or  mixtures. 

The  influence  can  be  great  as  is  shown  by  the  investigation 
of  La  additions  to  Ba1_x  Lax  ZrO^.34  For  x  =  0.06  the 
breakdown  strength  increased  175%,  apparently  clue  to  elec¬ 
tron-hole  recombinations  reducing  the  conductivity.  A 
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Typical  sample  and  electrode  configuration  for  dielectric 
breakdown  measurements. 
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similar  effect  was  observed  for  additions  of  Nb  in  place 
of  Zr. 

There  is  practically  no  specific  information  on  micro- 

structural  effects  on  breakdown  strength  known  to  the 

authors.  Aside  from  observing  that  extrusion  affects  the 

material  texture,  particle  alignment,  stress  distribution 

29 

and  hence  breakdown  strength  of  a  material,  the  only  de- 
tailed  paper  is  concerned  with  porosity.  The  investiga¬ 
tors  found  a  sharp  increase  in  the  statistical  distribution 
of  voids  in  the  direction  of  the  breakdown  path.  Experi¬ 
ments  with  Lucite  spheres  burned  out  of  PbZrOj  yielded 
fairly  good  agreement  with  the  theoretical  model  proposed 
for  porosity  effects. 

Investigations  concerned  with  the  influence  of  particle 
size,  glassy  phase,  and  phase  distributions  are  apparently 
absent  from  the  literature. 

Other  material  properties  that  should  be  examined  with 

respect  to  dielectric  breakdown  are  surface  condition  or 

grain  boundary  effects;  point  defects;  dislocations;  and 

dielectric  permittivity,  which  appears  to  be  related  to 

21 

breakdown,  in  BaTiOj.  There  is  evidence,  for  example, 

7  fL 

that  dislocations  increase  the  free  charge  in  dielectrics 
and  consequently  grinding,  polishing,  surface  treatments 
or  mechanical  history  can  influence  breakdown  behavior. 

Several  of  the  above-mentioned  experimental  parameters 
are  inter-related  and  others  might  also  be  mentioned,  but 
it  is  apparent  that  materials  studies  have  lagged  far 
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behind  theoretical  studies  of  dielectric  breakdown.  This 
is  unfortunate  because  the  theoretical  models  may  well  have 
overlooked  critical  materials  parameters  from  a  lack  of 
such  data. 


Conclusions 

The  conclusions  that  can  be  stated  are  similar  to  recom 
mendations  for  future  work.  Assumptions  made  in  the  theo¬ 
retical  approaches  often  limit  calculations  to  order  of 
magnitude  results,  although  including  the  effect  of  space 
charge  in  modifications  of  current  theories  may  produce 
improvements.  The  influence  of  prestressing  and  micro¬ 
structure  on  dielectric  breakdown  definitely  need  attention 
as  well  as  less  directly  related  phenomena,  such  as 
radiative  recombination.  Finally,  full  characterization 
of  the  materials  investigated  and  experimental  conditions 
must  be  reported  to  make  worthwhile  comparison  of  results 
practical . 
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B.  HETEROGENEOUS  SEMICONDUCTING  GLASSES  (H.  F.  Schaake) 

Introduction 

The  face  that  many  glasses  consist  of  several  phases 
in  metastable  equilibrium  has  been  recognized  for  many 

f' 

years.  The  separation  of  an  initially  homogeneous  glass 
into  two  or  more  phases  occurs  by  one  of  three  mechanisms: 
(1)  nucleation  and  growth  of  one  or  more  glassy  phases 
from  the  initially  homogeneous  phase;  (2)  spinodal  decom¬ 
position  of  the  initial  phase  into  two  new  phases;  and 
(3)  nucleation  and  growth  of  one  or  more  crystalline  phases. 
Each  of  these  processes  gives  rise  to  a  characteristic 
microstructure. 

That  the  existence  of  such  heterogeneities  affects  the 
electrical  properties  of  a  glass  has  been  demonstrated  for 
dielectric  glasses /^P^and  has  been  suggested  for  semi¬ 
conducting  glasses  It  appears,  however,  that  no  sys¬ 

tematic  study  of  the  effects  of  heterogeneities  on  the 
conductivity  of  a  semiconducting  glass  has  been  made.  This 
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situation  reflects  the  early  stage  of  development  of  the 
theory  of  homogeneous  amorphous  semiconductors,  and  the 
limited  structural  characterization  of  these  materials. 

An  example  of  a  composition  in  which  the  existence  of 
a  phase  separation  has  been  verified  is  the  701  y20s  * 

301  KPO^semiconducting  glass,  The  bulk  electrical 

properties  of  this  glass  in  both  quenched  and  heat  treated 
conditions  are  illustrated  in  Figure  1.  The  activation 
energy  for  DC  conduction  is  found  to  decrease,  while  the 
magnitude  of  both  AC  and  DC  conductivity  increases,  as  heat 
treatment  progresses.  Furthermore,  at  high  temperature, 
the  AC  conductivity  is  found  to  decrease  with  increasing 
frequency  in  the  heat  treated  glasses,  in  contrast  to  the 
increasing  conductivity  with  increasing  frequency  that  is 
found  in  the  quenched  glasses. 

Small  angle  scattering  x-ray  studies  have  confirmed 
the  occurrence  of  a  phase  separation  in  this  composition.^^ 
The  separation  is  present  in  the  quenched  glass  to  a  snail 
extent,  and  increases  markedly  with  small  amounts  of  heat 
treatment.  The  mean  particle  size  of  the  precipitated  phase 
is  found  to  increase  from  about  80  X  in  the  quenched  glass 
to  about  130  X  in  the  heat- treated  glasses.  At  the  same 
time,  the  volume  fraction  of  the  precipitate  increases,  but 
even  after  prolonged  treatment  the  volume  fraction  of 
crystallites  is  found  to  be  only  about  10\  (as  measured  by 
diffraction  in  a  Guinier  camera) .  Electron  microscopic 
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examination  of  glasses  in  this  system  have  also  confirmed 
the  occurrence  of  a  phase  separation  with  heat  treatment. ^ 
The  question  which  these  microstructural  studies  have 
raised  is  to  what  extent  can  the  observed  microstructural 
changes  affect  the  AC  and  DC  conductivity  of  the  bulk 
sample.  It  is  the  purpose  of  this  paper  to  review  the 
effect  of  a  heterogeneous  microstructure  on  the  AC  and  DC 
conductivity,  and  to  apply  these  results  on  data  obtained 
on  the  301  KPO^-  701  y2Ps\Sl*ss**« 


jfttrwtwt  HrtfU 

The  general  problem  of  the  electrical  behavior  of  a 
composite  in  which  random  nucleation  and  growth  leads  to 
intersection  of  the  new  phase  has  not  yet  been  treated. 

The  mixture  formulae  which  are  available  are  based  on 
specific  geometric  models,  none  of  which  are  applicable 
over  the  entire  range  of  volume  fractions  of  the  component 
phases  for  a  phase  separated  glass  A  semiquantitative 

understanding  of  the  electrical  properties  of  a  hetero¬ 
geneous  semiconductor  can  be  obtained,  however,  by  using 
an  idealized  model.  For  simplicity,  we  shall  consider 
the  precipitated  particles  to  be  non- intersecting  uniform 
spheres,  uniformly  distributed  in  the  matrix  phase. 

For  small  volume  fraction*  of  dispersed  particles,  the 
assumption  that  there  are  no  interactions  between  particles 
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is  applicable.  The  conductivity  of  such  a  composite  is 
given  bys^®^ 


Oj (2*2v2)  ♦  02(1+2V2-) 


0) 


where  o^is  the  conductivity  of  the  matrix  phase, 

is  the  conductivity  of  the  isolated  phase,  and 
Vjyi*  the  volume  fraction  of  the  isolated  phase. 

Assuming  that  the  conductivities  of  both  phases  can 
be  characterized  by  an  activation  energy,  the  apparent 
activation  energy  of  the  composite  may  be  calculated  from 
the  relation: 


dine'* 

— ST 


(2) 


Q* 


The  results  are: 

*  4«U-v2}  *  «*C1*2v2)1  ♦  9ay2. 
2C2"j2‘,';2x)  *  *  ^U*y2:2v22) 


(« 


with  a  ■ 

Three  distinct  regions  can  be  delineated  in  Eq.  (3): 


o  «  1  Q*  *  (4a) 

a  s  1  Q*  *  y2q2*Cl-V2)Qi,  (4b) 

o »  i  q*  =  91N  c«c) 


The  value  of  a  at  which  Eq.  (4a)  and  Eq.  (4c)  become  good 
approximations  to  Q*  is  dependent  on  the  volume  fraction, 


For  large  volume  fractions  of  the  isolated  phase,  the 
interactions  between  the  particles  cannot  be  ignored.  Hanai 
has  derived  an  equation  for  the  conductivity  valid  in  this 
limit  for  our  model : 


q ; 

°>1  *  q  2  ' 


(5) 


The  activation  energy  can  again  be  calculated  using  Eq.  (2): 

_  .♦  «i‘V-3  ♦  3?iTl/-3  «2Q2 

^ - -  C6) 


While  this  equation  is  slightly  more  complex  than  Eq.  (3), 
it  reduces  to  Eq.  (4)  in  the  same  limits.  These  results 
are  to  be  expected.  When  a  <<  1,  the  isolated  phase  is  an 
excellent  conductor,  and  the  conductivity  is  limited  by 
and  hence  determined  by  the  matrix  phase.  Finally,  when 
a  *  1,  both  phases  conduct  equally  well,  and  the  electric 
field  is  uniform  throughout  the  composite,  resulting  in 
the  simple  volume  fraction  dependence  of  the  activation 
energy. 

In  considering  the  AC  properties  of  the  DC  path  con¬ 
ductivity,  it  is  desirable  to  introduce  a  model  in  which 
the  isolated  particles  are  coated  with  a  thin  barrier  phase, 
as  is  shown  in  Figure  2.  Discussion  as  to  the  nature  of 
this  barrier  will  be  deferred  until  later.  The  sole  effect 
on  the  DC  conductivity  is  to  change  the  conductivity  of 
the  isolated  phase  to: 
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Fig.  2  Particle-barrier  model. 
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with  the  quantities  defined  in  Figure  2. 

The  effect  of  an  AC  field  on  the  conductivity  of  a 
composite  can  in  principle  be  computed  if,  in  addition  to 
the  real  conductivities  of  the  two  phases,  the  frequency 
behavior  of  the  dielectric  constant  is  known.  The  result¬ 
ing  equations,  however,  are  extremely  complicated,  even 
for  the  simplest  model.  The  reason  for  the  complication 
is  apparent  by  considering  the  distribution  of  the  electric 
field  in  the  composite.  If  a  »  1,  then  the  electric  field 
is  confined  almost  entirely  to  the  matrix  phase.  Hence  a 
relaxation  process  in  the  matrix  will  contribute  to  the 
AC  conductivity  much  more  greatly  than  will  the  identical 
process  in  the  isolated  phase.  Furthermore,  the  changing 
conductivity  with  frequency  will  result  in  a  change  in  the 
distribution  of  the  field  in  the  composite,  leading  to 
complex  frequency  dependent  interactions  between  the  two 
phases.  In  discussing  the  effects  of  frequency  dependent 
dielectric  constants,  we  shall  therefore  generally  confine 
our  attention  to  the  behavior  when  the  electric  field  is 
uniform  throughout  the  composite,  i.e.,  when  a  ■  1. 
Similarly,  when  discussing  the  effects  of  a  non-uniform 
electric  AC  field  in  the  composite,  we  shall  generally  treat 
the  dielectric  constant  as  being  frequency  independent. 


Mechanisms  and  Conductivity  in  Homogeneous  Glasses 


Before  proceeding  further  with  our  discussion  of 
heterogeneous  materials,  it  is  necessary  to  review  the  con¬ 
duction  mechanisms  which  have  been  proposed  in  the  homo¬ 
geneous  amorphous  materials  which  comprise  the  various 
phases . 

A  more  detailed  review  of  the  conduction  mechanisms 
in  homogeneous  semi-conducting  glasses  than  we  shall  present 
here  may  be  found  in  the  article  by  Davis  in  this  volume, 
in  several  papers  by  Mott,  and  in  various  papers  of 

the  SEAS  symposium.'^^" 

Three  distinct  mechanisms  have  been  recognized  as 
important  in  amorphous  materials  (Figure  3).  They  are: 

1.  Essentially  normal  band  conduction  for  carriers 
having  energies  above  some  critical  value,  E^. 

*2.  A  diffusion  mechanism,  in  which  the  mobility  of 
a  carrier  having  an  energy  below  Ec\but  above  a  lower  energy 
E^v  is  considerably  reduced  due  to  the  electronic  wave- 
function  being  highly  modulated  by  the  disorder  in  the 
lattice,  which  results  in  regions  in  the  diffusion  path 
where  the  probability  of  finding  the  carrier  is  extremely 
small.  The  behavior  of  the  carriers  in  the  diffusion  region 
is  still  band- like,  in  that  it  is  a  non- thermally  activated 
process.  There  may  be  a  substantial  interaction  between 
the  electron  in  these  states  and  the  lattice  displacement 
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Fig-  3  Conductivity  mechanisms  in  amorphous  semiconductors.  The  density  of  states, 
plotted  against  the  energy,  E.  Symbols  are  defined  in  the  text. 


which  is  polaronic  in  nature.  The  importance  of  the  inter¬ 
action  will  depend  on  the  various  material  parameters 
involved.  For  a  further  consideration,  see  Reference  3,  in 
which  this  diffusion  region  is  referred  to  as  a  polaronic 
band.  We  may  conclude,  however,  that  the  energies  and 
E^are  temperature  dependent,  increasing  with  increasing 
temperature. 

3.  Carriers  having  energies  below  p^will  be  localized 
in  the  lattice,  and  transport  will  be  by  a  thermally- 
activated  hopping  process. 

Assuming  the  carriers  to  be  electrons,  the  DC  conduc¬ 
tivity  of  a  homogeneous  material  may  be  written  as  follows: 

o  ■  /'b  en(E)^E)dE  ♦  en(E)ii^(E)dE 
EB\ 

♦  f  3  en(E)ySv(E)dE  (8) 

?c\ 

where  e  is  the  electronic  charge 

n(E)  is  the  density  of  carriers  with  energy  E 
l^(E)  is  the  mobility  of  the  carrier  in  the  hopping 
region 

p^E)  is  the  mobility  of  a  carrier  in  the  diffusion 
region 

lj£(E)  is  the  mobility  of  a  carrier  in  the  band  region 
and  E/^\is  the  enersy  at  t^ie  t0P  conduction  band. 
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When  an  AC  field  is  applied  to  a  homogeneous  amorphous 
semiconductor  containing  partially  filled  localized  states, 
localized  carriers  which  are  essentially  isolated  from  the 
DC  path  can  jump  to  nearby  localized  states  of  nearly  the 
same  energy,  giving  rise  to  a  relaxation  loss  peak.  The 
real  part  of  this  AC  conductivity,  which  is  additive  to 
the  DC  conductivity  will  be  given  by:'^^^-* 

R.  e^)  -  t  C(t)  (9) 


where  t  is  the  relaxation  time 

C(t)  is  a  proportionality  factor. 

The  complex  part  of  the  conductivity  is  related  through  a 
Kramers -Kronig  relation. 

The  relaxation  time  will  be  given  by  s'®'*’ 

t  ■  ro  ttnh  TFT  -Jr  >>  T  (10) 


or,  at  sufficiently  high  temperatures: 

r  -  C  exp t -(Ex.-  “  -  W^/kT] 

BjK  hwA 

j—  COth  jfo  »  1 

✓  °  • 

whore  EA\is  tho  activation  energy 

AE  is  the  difference  in  energy  between  initial  and 
final  state 

and  (j^is  the  optical  phonon  frequency. 


«  T;  (11) 
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The  presence  of  heterogeneities  in  the  semiconducting 
glass  gives  rise  to  the  possibility  of  several  additional 
conduction  mechanisms,  which  we  shall  now  discuss. 

The  classic  mechanism  in  a  heterogeneous  medium  is  the 
interfacial  polarization  or  Maxwell-Wagner  loss  mechanism, 
which  arises  when  the  conductivities  or  dielectric  constants 
of  the  two  phases  differ S!)*'  When  a  step  field  is  applied 
to  such  a  heterogeneous  medium,  the  electric  field  dis¬ 
tribution  in  the  various  phases  immediately  after  the 
application  of  the  field  (which  is  determined  by  the  dielec¬ 
tric  constants),  will  differ  from  the  field  distribution  at 
infinite  times  (which  is  determined  by  the  conductivities). 
The  transition  from  the  t  ■  Cr  regime  to  the  t  ♦  •  regime 
is  a  relaxation  process  dependent  on  the  geometry  of  the 
heterogeneous  medium,  as  well  as  the  electrical  constants. 

The  effect  of  the  inhomogeneities  on  an  AC  field  may 
be  calculated  for  the  low  concentration  model  of  Section  II 
by  replacing  the  real  conductivity  in  Eq.  (1)  with  the 
complex  admittance: 


o*  ♦  iue.\ 

,'i\  / 1 N 


(12) 


The  results  for  the  real  (total)  conductivity  when  the 
dielectric  constants  are  frequency  independent  are: 


£  »/  ♦ 


C*s'  * 


%  'V 

«  T 


1  ♦ 
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w  - T 


(is) 


84 


where : 


(2*v2)Cj  ♦  Cl-v2)e2 
CZ.v2)Oi  ♦  U-v2j«2. 


(M) 


■*\ 


Ry»2)t«2-«i)  •  (»«i**2»?i-«ii' 

^  S'2H  U>y2l^Mi-y2)?2]^  J 


e^is  given  by  Eq.  (1)  with  ej\substituted  for  o^,  and  6'S 
is  given  by  Eq.  (1) . 

In  order  for  the  losses  to  be  observable,  there  need 
only  be  an  observable  difference  between  e  and  .  It  will 
be  noticed  that  Eq.  (13)  is  the  equation  for  a  simple  Debye 
relaxation  process.  This  is  no  longer  the  case  if  we 
consider  the  high  concentration  model  of  Section  II,  or 
one  in  which  the  isolated  phase  is  no  longer  spherical. 

The  general  characteristics  of  the  loss  peak  remain  the 
same,  however,  with  other  powers  of  u,  in  addition  to 
appearing  in  the  equation  for  o(u). 

An  interesting  model  is  that  in  which  each  of  the 
isolated  particles  is  coated  with  a  barrier  phase.  Treating 
this  model  in  the  low  concentration  limit,  and  assuming 
iwc;  a ^  »  a^t  real  (see  Figure  2),  we  find  that 
Eq.  (13)  is  applicable  and  that 


We  shall  refer  to  these  results  later. 

In  addition  to  Maxwell -Wagner  losses,  several  other 
AC  phenomena  may  be  present  in  heterogeneous  semiconductors. 
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These  additional  mechanisms  depend  on  the  existence  of 
Schottky  barriers  or  other  types  of  barriers,  such  as  a 
compositional  change  at  the  interface  of  the  two  phases. 

The  Schottky  barriers  will  arise  when  the  Fermi  levels  in 
the  different  phases  in  isolation  are  different,  due  to 
differences  in  composition,  structure,  or  both.  In  order 
to  establish  a  constant  Fermi  level  when  the  two  are  in 
contact,  charge  may  be  transferred  from  the  phase  having 
the  higher  Fermi  level  to  that  having  the  lower  (Figure  4a). 
It  is  these  barriers  which  give  rise  to  the  barrier  phase 
of  Section  II. 

With  the  presence  of  these  barriers,  the  following 
additional  mechanisms  may  arise. 

1.  Barrier  Tunneling.  If  the  charge  transfer  is 
sufficient  to  completely  displace  the  conduction  states 
which  are  contributing  the  greatest  to  the  conduction  or, 
if  the  isolated  phase  is  separated  from  the  matrix  phase 
by  an  insulating  barrier,  then  barrier  tunneling  will 
result  (Figure  4).  If  conduction  is  by  carrier  hopping, 
then  the  hopping  time  through  the  barrier  will  be  given  by: 
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If  the  conduction  is  predominantly  by  a  band  diffusion 
mechanism,  then  band  tunneling  will  result.  A  successful 
quantitative  analytic  treatment  of  this  problem  in  crystal* 
line  semiconductors  has  not  been  made;  the  results  which 
are  available  depend  on  the  way  in  which  a  localized  wave- 
function  is  constructed  from  the  band  states,  and  on  the 
existence  of  localized  states  in  the  barrier. 

The  DC  conductivity  of  the  barrier  phase  may  be  approx¬ 
imately  calculated  from  the  relation: 

?bv“  /6  «xp(-2ad)  (17) 

where  oON is  the  conductivity  in  the  absence  of  the  barrier. 
The  results  under  an  applied  AC  field  are  more  interesting. 
Since  a  transition  through  the  barrier  takes  a  finite  amount 
of  time,  and  if  the  primary  DC  path  is  through  the  barrier 
(which  is  the  case  if  o2  >  Op  with  c^\defined  in  Eq.  (7)), 
then  as  the  period  of  the  applied  field  exceeds  the  transi¬ 
tion  time,  the  DC  path  conductivity  will  decrease.  Such 
a  circumstance  may  lead  to  a  decrease  in  the  total  conduc¬ 
tivity  with  increasing  frequency.  The  DC  path  conductivity 
should  then  be  approximately  given  by: 

o  ■  oj.(1-v2)n  ♦  [o’-o(l-v2X]  (»itx*Fy  (18) 

If,  on  the  other  hand,  the  isolated  particles  are 
essentially  isolated  from  the  DC  path  (^2v«  o  ^ ,  then  an 
additive  contribution  to  the  AC  conductivity  should  be 
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present,  in  accordance  with  Eqs.  (9)  and  (16),  due  to  the 
transition  of  carriers  through  the  barrier.  This  additive 
contribution  will  be  in  addition  to  the  Maxwell -Wagner  effect. 

2.  Thermal  Activation  Over  the  Barrier.  If  the  charge 
transfer  does  not  completely  displace  the  conduction  states 
which  are  contributing  the  greatest  to  the  conduction 
(Figure  4a) ,  then  thermally  activated  hopping  over  the 
barrier  may  result.  The  conductivity  through  the  barrier 
will  then  be  given  by: 

c  -  nejjg^  exp  (-W/kT)  (19) 

where  the  mobility  jij^is  essentially  unchanged  from  that 
of  the  constituent  material.  Since  the  mobility  is  essen¬ 
tially  unchanged,  such  a  process  will  make  no  additional 
contribution  to  the  AC  conductivity  --it  will  merely  affect 
the  DC  conductivity. 


Application  to  Hotoroaeneouo  Claopeo 

As  discussed  in  the  introduction,  glasses  with  the 
composition  30t  KPO^-  701  V^O^arc  found  to  be  heterogeneous 
after  a  short  heat  treatment  at  288*C.  These  glasses  (which, 
according  to  thermopower  measurements,  are  n-type  semicon¬ 
ductors)  are  therefore  potentially  useful  in  applying  the 
principles  of  the  previous  sections.  Our  task  will  be 
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complicated  by  the  lack  of  electrical  data  on  each  of  the 
constituent  phases. 

We  shall  begin  by  considering  the  AC  conductivity,  as 
this  gives  the  most  information  as  to  the  mechanism  opera- 
titf.  The  additive  AC  conductivity  for  a  quenched  and 

A 

heat  treated  glass  is  presented  in  Figures  S,  6  and  7.  To 
elucidate  the  conduction  mechanisms,  it  is  worthwhile  to 
note  that  the  AC  conductivity  above  SxlO^MHz  is  proportional 
to  both  the  quenched  and  heat  treated  glasses,  and 

that  the  activation  energy  for  the  AC  conductivity  is 
approximately  0.0S  ev  in  the  quenched  glass  at  all  frequen¬ 
cies,  and  in  the  heat  treated  glass  at  low  frequencies. 

This  is  suggestive  of  the  hopping  of  trapped  carriers  with 
x  appreciably  smaller  than  3xl0^/sec. ,  as  Eq.  (9)  reduces 
to: 


oA^«)  -  Z  C(t) 


\2/^ 

W*T 


(20) 


when  the  product  ut  <<  1.  The  temperature  dependence  of 
Eq.  (20)  is  contained  in  the  product  C(t)t.  The  temperature 
dependence  in  C(t)  is  due  to  the  product  of  the  Fermi 
factor  f (E)  and  a  distribution  function  g(E,  AE)  which  tells 
how  many  pairs  of  states  at  energy  E  have  a  difference  in 
energy  AE,  and  a  factor  1/T  which  arises  from  the  Einstein 
equation  for  the  diffusion  coefficient,  while  the  temperature 
dependence  of  t  is  given  by  Eq.  (10)  or  Eq.  (11).  The 
present  data  is  insufficient  to  separate  these  temperature 


90 


FREQUENCY,  Hx 


Fig.  5  AC  conductivity  in  quenched  and  heat  treated  glasses 
(70%  V205  -  30%  KP03)  at  300#K. 
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Fig.  6  Temperature  dependence  of  AC  conductivity  in  quenched 
glass. 
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dependences.  Owen,  however,  has  measured  the  AC  conduc¬ 
tivity  of  similar  glasses  and  has  found  t  to  be  of  the 
order  of  lCT^lVsec.  and  essentially  independent  of  tempera¬ 
ture.^1^  It  therefore  appears  that  all  of  the  temperature 
dependence  is  in  the  factor  C(t). 

If  the  high  frequency  portion  of  the  AC  conductivity 
is  due  to  thermally  activated  hopping  of  trapped  charges, 
these  results  force  us  to  the  following  conclusions: 

1.  That  the  interaction  between  the  localized  elec¬ 
tron  and  the  lattice  displacement  is  very  small.  This  in 
turn  means  that  localized  wavefunction  has  a  substantial 
radius  in  comparison  with  an  interatomic  distance. 

2.  That  localized  states  occur  either  isolated,  or 
in  pairs,  with  a  well  defined  distance  and  small  energy 
difference  between  the  pairs.  This  latter  situation  is 

a  result  of  the  fact  that  if  the  localized  stat.'s  are  ran¬ 
domly  distributed,  then  the  AC  conductivity  should  be 

\2  'fin 

approximately  proportional  to  w  rather  than  to 

3.  That  the  hopping  is  adiabatic  (see  Appendix). 

The  fact  that  there  is  a  substantial  activation  energy 
for  DC  conduction  (0.2  to  0.6  ev,  Figure  8),  in  conjunction 
with  the  above  conclusions,  indicates  that  the  DC  conduction 


is  determined  by  a  mechanism  other  than  thermally  activated 
hopping.  The  width  of  the  energy  gap,  however,  is  of  the 
order  of  2.5  ev  for  these  materials /^Vindicating  that 


they  are  extrinsic,  rather  than  intrinsic  semiconductors 
(with  V4^  acting  as  donors) . 
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Fig.  8  Temperature  dependence  of  DC  conductivity  in  quenched  and 
heat  treated  glass.  The  dashed  line  is  explained  in  text. 


We  are  therefore  led  to  postulate  two  models  for  the 
conduction  mechanisms  in  the  glass: 

1.  (Figure  9a)  That  due  to  the  disorder,  a  large 
number  of  localized  states  are  split  off  from  the  conduction 
band,  with  a  large  concentration  located  in  the  vicinity 

of  the  impurity  states,  ~0.2  ev  below  the  conduction  band, 
and  with  barriers  of  from  0  to  0.4  ev  present,  depending 
on  the  heat  treatment.  The  high  frequency  AC  conductivity 
is  then  due  to  impurity  hopping,  with  the  conditions  pre¬ 
viously  discussed  being  true.  The  DC  conduction  is  barrier- 
regulated  by  a  band  mechanism. 

2.  (Figure  9b)  That  the  impurity  states  lie  just 
underneath  the  conduction  band,  but  that  there  are  barriers 
raised  in  the  conduction  band  which  determine  the  DC  acti¬ 
vation  energy.  The  DC  conduction  is  therefore  a  barrier- 
regulated  band  mechanism,  while  the  AC  conductivity  results 
from  a  migration  of  carriers  at  the  bottom  of  the  band 
(which  are  localized  states  when  viewed  over  the  entire 
glass)  under  the  influence  of  the  AC  field.  The  height  of 
the  barriers  must  decrease  with  heat  treatment. 

Turning  now  to  the  low  frequency  AC  conductivity,  we 
see  that  in  the  quenched  glass  there  is  a  region  where  the 
conductivity  is  proportional  to  w,  while  in  the  heat 
treated  glass  there  appears  to  be  a  well  defined  relaxation 
process  with  x  *  SxlO’/^sec.  There  are  two  explanations 
for  the  behavior  in  the  quenched  glass.  (1)  The  hopping  of 
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le  electronic  energy  structures  in  70%  V2O,.  -  30%  KPO, 
the  Fermi  level. 


localized  carriers  between  randomly  distributed  centers. 

This  would  be  in  addition  to  the  pair  hopping  previously 
discussed;  a  much  smaller  number  of  carriers  would  be 
participating  in  this  process  than  in  the  pair  hopping 
process.  (2)  Maxwell -Wagner  losses  due  to  a  non-spherical 
second  phase.  There  is  some  x-ray  evidence  that  the  iso¬ 
lated  phase  in  the  quenched  glass  is  non-spherical.^ 

In  the  heat  treated  glass ,  however,  only  the  Maxwell- 
Wagner  loss  mechanism  explains  the  data.  The  relaxation 
time  for  this  loss  is  found  to  be  a  function  of  temperature, 
with  an  activation  energy  of  0.2  ev  (Figure  10) .  This  is 
also  the  activation  energy  for  DC  conduction  in  the  composite. 

Using  the  low  concentration  model,  and  the  assumptions 
that  a  >>  1,  and  that  the  dielectric  constants  of  the  two 
phases  are  temperature  independent  in  the  region  of  inter¬ 
est,  the  relaxation  time  should  be  proportional  to  a. 

In  the  low  concentration  model  when  the  particles  are 
coated  with  a  capacitive  barrier  phase  (Eq.  (15)),  the 
relaxation  time  is  proportional  to  A  physical  consid¬ 

eration  of  this  latter  model  in  the  limits  of  high  concen¬ 
tration  and  with  the  barrier  phase  having  a  small  real 
conductivity,  indicates  that  the  relaxation  time  should  be 
proportional  to  the  reciprocal  of  the  composite  conductivity, 
for  the  barrier  forms  a  series  RC  circuit,  with  the  resis¬ 
tance  being  determined  essentially  by  the  entire  composite. 
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Turning  now  to  the  DC  path  conductivity  in  the  heat 
treated  glass,  it  will  be  noticed  that  at  high  temperatures 
and  frequencies,  the  conductivity  decreases  with  increasing 
frequencies  of  lO^and  4xlO^MHz  (Figure  8) .  The  conduc¬ 
tivity  of  the  composite  at  4  MHz  less  the  additive  AC 
conductivity  (which  is  evident  at  low  temperatures  and  is 
assumed  to  be  present  also  at  high  temperatures)  is  plotted 
in  Figure  10.  It  will  be  noticed  that  this  forms  an  exten¬ 
sion  of  the  low  temperature  DC  conductivity.  This  data  is 
consistent  with  a  model  in  which  the  isolated  phase  is 
coated  with  a  tunneling  barrier,  for  which  the  transit  time 
through  the  barrier  is  -10"?  seconds.  The  activation  energy 
of  the  matrix  is  0.05  ev,  while  that  of  the  isolated  phase 
is  of  the  order  of  0.2  ev.  At  330*K,  the  conductivity  of 
the  isolated  phase  (particle  plus  barrier)  becomes  equal 
to  that  of  the  matrix  phase,  marking  the  onset  of  the 
decreasing  conductivity  with  frequency.  With  this  inter¬ 
pretation,  we  would  expect  the  total  conductivity  to  level 
off  at  frequencies  above  4.5  MHz,  and  then  to  again  increase 
with  frequency,  due  to  the  domination  of  the  carrier  hopping 
contribution  to  the  total  conductivity. 

Attempts  to  elucidate  the  nature  of  the  barriers  are 
of  questionable  value,  due  to  the  limited  structural  and 
electrical  data,  and  an  incomplete  knowledge  of  the  elec¬ 
tronic  structure  of  the  component  phases.  Schottky  barriers 
may  be  justified  on  the  basis  of  the  difference  in  the 
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Fig.  1 0  Loss  peaks  in  heat  treated  glass. 


composition  of  the  matrix  and  isolated  phases.  On  the 
other  hand,  insulating  barriers  could  exist  as  a  result  of 
a  vanadium  deficiency  at  the  interphase  boundary,  which 
results  from  a  minimization  of  the  interfacial  energy. 
Further  structural  and  electrical  studies  will  be  required. 

The  discussion  applied  to  the  data  of  Figures  5,  6, 

7,  and  8  does  not  appear  to  apply  to  all  data  on  heterogen¬ 
eous  glasses,  however.  For  example,  the  decrease  in  AC 
conductivity  with  frequency  illustrated  for  the  glasses 
heat  treated  0.5  hours  and  30.8  hours  in  Figure  1  is  not 
the  same  as  that  just  discussed,  in  that  the  DC  conduction 
process  does  not  show  a  break  in  activation  energy.  An 
explanation  for  this  difference  may  lie  in  differences 
between  the  microstructures  of  these  glasses.  A  liquid- 
liquid  separation  and  a  liquid-crystal  separation  may  be 
competing  with  each  other.'CS^  In  that  these  glasses  were 
cast  from  different  temperatures,  their  microstructures 
before  heat  treatment  may  be  different,  in  which  case  the 
kinetics  and  microstructural  development  during  heat  treat¬ 
ment  will  be  different.  Such  a  possibility  is  of  consid¬ 
erable  interest  in  that  the  activation  energies  for  DC 
conduction  differ  in  the  two  quenched  glasses  (0.4  ev  for 
that  in  Figure  1,  0.6  ev  for  that  in  Figure  8). 
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Appendix 


In  the  derivation  of  the  transition  times  (Eqs.  (10) 
and  (11)) ,  it  was  assumed  that  the  value  of  the  overlap 
integral,  J,  between  the  initial  and  final  states: 

J  *  <m|V|n> 

is  sufficiently  small,  so  that  the  transition  takes  place 
over  several  periods  of  the  phonon  frequency.  This  assump¬ 
tion  leads  to  a  dependence  of  the  transition  time  for  a 
one  phonon  process: 

_  „  1 
W 

If,  however,  J  is  sufficiently  large,  the  conditions  for 
the  approximations  made  in  Ref.  3  are  no  longer  valid. 

This  problem  (adiabatic  hopping)  has  been  investigated  for 
the  multiphonon  process  by  Holstein,'H5^who  finds^the 
pre-exponential  factor  in  Eq.  (11)  becomes  approximately 
constant,  independent  of  the  overlap  integral.  Physically, 
this  means  that  the  time  it  takes  the  carrier  to  tunnel 
from  one  site  to  another  is  small  in  comparison  to  the 
time  it  takes  for  a  phonon  to  be  absorbed  or  emitted;  hence 
the  phonon  coupling  process  becomes  the  rate  determining 
mechanism.  We  expect  the  same  analysis  to  be  true  for  a 
one  phonon  process  (Eq.  (10)). 
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The  existence  of  adiabatic  hopping  means  that  the 
transition  time  between  sites  having  a  separation  greater 
than  some  critical  value  r^will  be  dependent  on  the 
separation,  while  the  transition  time  for  hops  between 
sites  having  a  separation  less  than  r^will  be  approxi¬ 
mately  independent  of  the  separation. 

The  AC  conductivity  will  therefore  consist  of  two 
components;  one  proportional  to  u  due  to  non-adiabatic 
hopping,  and  one  proportional  to  ue  to  adiabatic 
hopping.  If  r.^is  sufficiently  great,  such  that  there 
are  a  large  number  of  site  pairs  contributing  to  adiabatic 
hopping,  then  the  iib&'dependence  will  dominate  at  high 
frequency.  This  is  the  behavior  found  in  the  quenched 
glass  (Figure  S) . 
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C.  POWDER  SAMPLE  PREPARATION  ERRORS  IN  X-RAY  SPECTROCHEMICAL  ANALYSIS 
(D.  E.  Clark  and  L.  L.  Hench) 

Introduction 

In  recent  years  the  use  of  x-ray  spectrochemical  analysis 
has  become  very  popular  in  the  fields  of  mining,  glass,  and 
ceramics  for  determining  the  relative  composition  of  ores,  raw 
materials,  and  finished  products.  In  each  of  these  cases  analy¬ 
ses  are  made  on  particulate  samples.  Errors  in  x-ray  spectro¬ 
chemical  analysis  are  due  to  fluctuations  in  equipment,  mis-use 
of  equipment,  or  problems  in  sample  preparation.  This  paper 
is  concerned  primarily  with  the  error  resulting  from  sample 
preparation. 

Neglecting  the  ever-present  problem  of  obtaining  a  repre¬ 
sentative  10-20  gram  sample  from  a  10  ton  stockpile,  sample 
preparations  generally  stands  out  as  the  most  difficult  and 
time-consuming  step  in  the  entire  analytical  process.  The 
x-ray  spectrochemical  method  may  be  applied  to  a  sample  in  the 
liquid,  solid,  glassy  or  powdered  state.  Consequently,  there 
are  many  methods  of  sample  preparation.  We  are  concerned  here¬ 
in  with  only  the  problems  associated  with  the  general  method 
of  pressing  dry  powders  into  small  pellets  and  exposing  them 
to  radiation. 

All  matter  absorbs  x-rays  to  a  greater  or  less  extent 
depending  on  (1)  the  wavelength  of  the  x-rays,  and  (2)  the  mass 
absorption  coefficient,  u/p,  of  the  materials  for  that  partic¬ 
ular  wavelength.  For  completely  solid  materials  with  no  voids 
the  absorption  equation  relates  the  transmitted  intensity  of 
the  x-ray  beam,  It,  to  the  initial  intensity  before  striking 
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the  sample  (IQ) , 

It  -  IQ  exp(-(y/p)  pt) 

where  "y"  is  the  linear  absorption  coefficient,  "p"  is  the 
density  of  the  material,  "tM  is  the  sample  thickness,  and  "u/p" 
is  the  mass  absorption  coefficient.^  For  the  heavier  elements 
u/p  is  such  that  the  incoming  x-rays  penetrate  only  a  few  ang¬ 
stroms.  However,  in  the  lighter  elements,  such  as  magnesium, 
u/p  is  small  and  thus  x-rays  can  penetrate  to  a  considerable 
depth  in  the  sample.  As  much  as  801  porosity  can  occur  in 
pressed  powder  samples  at  low  forming  pressures  and,  consequently, 
the  mass  absorption  coefficient  (y/p)Q  for  the  pure  compound  is 
changed  to  (y/p^  for  the  compound  plus  (y/p)2  for  pores  between 
the  particles.  Of  course,  the  coefficient  of  absorption  for 
pores  is  much  less  than  for  a  solid  material.  Thus,  one  would 
expect  the  x-rays  to  penetrate  a  pressed  powder  sample  to  a  greater 
depth  than  the  pure  crystalline  material  with  no  porosity.  The 
accuracy  of  x-ray  spectrochemical  analyses  is  thus  potentially 
highly  susceptible  to  variations  in  density  and  thickness  of 
pressed  samples.  This  potential  source  of  analytical  error  is 
examined  in  this  paper. 

The  objectives  of  the  research  are:  (1)  To  determine  the 
characteristic  x-ray  intensity  obtained  from  pressed  powders  as 
a  function  of  forming  pressure  and  sample  density.  (2)  To 
determine  the  x-ray  intensity  of  pressed  powder  samples  as  a 
function  of  sample  thickness. 
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Experimental  Procedure 

Four  chemical  compounds,  MgCO^ ,  Ti02,  Cr203,  Th02,  were 
studied.  All  four  compounds  were  reagent  grade  and  were  200 
mesh  or  finer.  Samples  were  formed  by  pressing  loose  powders 
in  a  single  action  cylindrical  Carver  die  with  a  pressing  sur¬ 
face  area  of  one  square  inch.  Samples  were  prepared  with 
forming  pressures  from  2,000  psi  to  24,000  psi  at  increments 
of  2,000  psi.  A  second  series  of  low  density  samples  was  pre¬ 
pared  by  compacting  loose  powders  directly  in  the  x-ray  sample 
holder.  The  sample  holder  was  filled  with  powder  and  exposed 
to  x-radiation.  The  powder  was  pressed  at  50  psi  by  hand  and 
irradiated  again.  This  process  was  performed  in  increments  of 
50  psi  up  to  --200  psi.  A  sheet  of  1/4  mil  mylar  was  placed 
over  the  end  of  the  sample  holder  to  prevent  the  loose  powders 
from  falling  into  the  x-ray  unit.  The  same  sheet  of  mylar  was 
used  for  the  pressed  pellets  as  well. 

A  Norelco  Vacuum  X-ray  Spectrometer  was  used  for  the  experi¬ 
ments.  The  samples  were  exposed  to  Tungsten  La  radiation  and 
the  time  necessary  to  count  a  fixed  number  of  characteristic 
x-rays  was  recorded.  The  fixed  count  divided  by  the  time  neces¬ 
sary  to  produce  the  count  is  reported  as  the  measured  intensity. 
The  counting  times,  sample  thicknesses  and  experimental  condi¬ 
tions  employed  are  summarized  in  Table  I. 

The  intensities  of  the  characteristic  MgCOj  and  Ti02  radi¬ 
ation  were  measured  with  a  modified  flow-gas  proportional 
counter.  90%  argon- 10%  methane  gas  flowed  through  the  counter 
at  the  rate  of  0.5  cfh.  This  type  of  detector  is  very  effective 


108 


O 

N  <t  N  N 
©  I  •  •  •  • 

fvi|  vo  vo  04  t" 

©  to  to  04 


44 

c 

a 

o 

u 

o 

o 

o 

O 

o 

o 

o 

o 

T3 

o 

o 

o 

o 

in 

a> 

<«*■ 

00 

«*■ 

* 

X 

VO 

VO 

3 

•H 

•n 

(JU 

* 

o 

a 

(/> 

G 

<D 

o 

in 

E 

Vi 

o 

3 

U 

o 

3 

•H 

1 

» 

1 

iH 

U 

E 

1 

1 

I 

(4 

1 

1 

1 

0) 

> 

o 

45 

o 

«H 

<44 

O 

G 

G 

o 

o 

44 

•H 

•H 

45 

J 

4-* 

4-> 

to 

O 

45 

45 

(4 

(4 

•H  • 

«H 

<w 

<44 

iH 

iH 

3  u 

U« 

u 

u 

iH 

iH 

45  0> 

•i4 

•H 

*3 

</> 

U) 

to 

4-* 

4-* 

a>  iH 

«4 

• 

• 

C 

G 

45  O 

O 

o 

o 

•H 

•H 

44  45 

U 

O 

in 

in 

O  4) 

4->  .-4 


iH 

v>  E 

<4 

•3  t4 

44 

H 

H 

G  in 

t/J 

a. 

o 

Q 

bu 

O 

>s 

Q 

o 

Q 

•H 

0-  >S 

U 

< 

w 

W 

m  3 

u 

3  M 
M  i 

Vi  X 

•3 

O 

G 

U  3 

3 

4S 

o 

to 

to 

i-4  44 

a 

o 

N 

O 

04 

00 

B 

u 

o 

OJ 

o 

*  G 

o 

40 

•H 

Vi 

45 

CM  *rt 

u 

X 

H 

CJ 

H 

« 

109 


for  elements  with  low  atomic  number  (Z  <  20} .  The  characteris¬ 
tic  intensities  of  C^Oj  and  ThC^  were  measured  with  a  scintil¬ 
lation  detector. 

Experimental  Results 

The  increase  in  bulk  density  of  the  samples  versus  forming 
pressure  is  shown  in  Figure  1.  The  density  change  is  repre¬ 
sented  as  the  measured  fraction  of  theoretical  (or  single  crys¬ 
tal)  density  of  the  compound.  Initially,  the  loose  powders 
contain  about  60-90%  porosity.  Cold  pressing  at  24,000  psi 
produces  samples  with  a  range  of  40-60%  porosity.  Figure  1 
shows  that  the  powders  have  an  initial  large  increase  in  density 
with  increasing  pressure.  The  higher  density  compounds  appear 
to  reach  their  maximum  density  at  lower  forming  pressures  than 
do  compounds  containing  lighter  elements. 

Magnesium  carbonate  has  the  most  interesting  pressing 
characteristics  of  the  materials  studied.  Its  pressure-density 
curve  has  the  smallest  initial  slope  and  reaches  its  maximum 
at  much  higher  pressures  than  the  other  materials.  Pressing 
above  10,000  psi  causes  an  interesting  phenomenon  to  occur  in 
thick  samples  of  magnesium  carbonate.  A  cone  shaped  portion  of 
the  pellet  can  actually  be  separated  from  the  rest  of  the  sam¬ 
ple.  This  is  caused  by  shearing  due  to  internal  density 
gradients,  as  discussed  by  Train. ^ 

The  x-ray  intensity-density  curves  for  all  four  compounds 
are  similar  as  shown  in  Figure  2.  There  are  three  distinct 
regions  that  occur  in  the  curves.  Region  I  corresponds  to  a 
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range  of  densities  in  which  the  change  in  intensity  with  density 
is  very  small.  The  second  region  (II)  in  the  density-intensity 
curve  shows  that  a  small  change  in  density  results  in  a  large 
change  in  intensity.  This  effect  is  highly  pronounced  in 
magnesium  carbonate.  At  high  densities  another  regime  of 
nearly  independent  x-ray  intensities  appear  (Region  III).  In 
the  case  of  Region  I  is  preceded  by  another  inflection 

in  x-ray  intensity,  denoted  as  Region  IA. 

The  range  in  density  over  which  the  three  regions  occur 
and  the  incremental  change  in  x-ray  intensity  (AI)  accompanying 
them  is  summarized  in  Table  II.  An  average  measured  intensity 
(I  ),  independent  of  density  fluctuations,  is  also  given  in 
Table  II.  The  ratio  AI / Ift  represents  a  potential  analytical 
error  if  samples  are  pressed  and  analyzed  without  control  of 
the  pressed  density.  A  summary  of  the  density  dependent  errors 
for  each  of  the  regions  of  the  compounds  studied  is  also  given 
in  Table  II . 

The  magnesium  carbonate  samples  showed  another  very  inter¬ 
esting  phenomenon.  Figure  3  shows  that  the  x-ray  intensity 
from  samples  pressed  at  2,000  psi,  with  a  density  range  of 
0.73  g/cc  to  0.75  g/cc,  exhibited  a  significant  dependence  on 
the  thickness  of  the  pellet.  The  intensity  increases  from  65 
cps  to  73  cps  as  the  thickness  is  increased  from  0.4  cm  to 
1.51  cm.  The  same  effect  was  observed  for  samples  pressed  at 
6,000  psi  with  a  density  range  of  0.93  g/cc  to  0.95  g/cc.  The 
intensity  increased  from  65  cps  to  76  cps  as  the  thickness 
increased  from  0.2  cm  to  0.75  cm.  Samples  irradiated  from  the 
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TABLE  II 


Compound 

Region 

pth* 

MgC°3 

I 

8-22 

II 

22-36 

III 

36-43 

Ti02 

I 

14-36 

II 

36-41 

III 

41-54 

Cr2°3 

IA 

15-20 

I 

20-44 

II 

44-49 

III 

49-58.5 

Th02 

I 

19-36 

II 

36-42 

III 

42-58.2 

*1  of  theoretical  density 


AI  Average  I  %  Error 


12 

cps 

12  cps 

23 

55 

80 

69 

2 

105 

2 

50 

cps 

275 

18 

60 

330 

18 

10 

360 

3 

500 

cps 

2850 

18 

10 

3100 

3 

400 

3450 

12 

3800 

0 

10 

cps 

1850 

5 

1100 

2500 

44 

20 

3150 

6 
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surface  directly  exposed  to  the  plunger  in  the  pressing  opera¬ 
tion  give  correspondingly  higher  intensities  than  those  irradi¬ 
ated  from  the  bottom  of  the  sample,  away  from  the  plunger. 

The  other  three  compounds  investigated  showed  no  dependence 
of  intensity  on  thickness  or  variation  of  intensity  when  irradi¬ 
ated  on  the  top  and  bottom. 

Discussion  of  Results 

Magnesium  carbonate  is  considered  separately  from  the 
others  due  to  its  unusual  thickness  dependence.  The  explanation 
for  its  behavior  is  based  on  the  fact  that  magnesium  is  a  less 
efficient  x-ray  absorber  in  comparison  with  heavier  elements. 
Therefore,  the  x-ray  beam  penetrates  to  an  appreciable  depth 
in  the  sample  before  being  totally  absorbed.  Figure  4a  illus¬ 
trates  the  point  that  when  the  incident  radiation  penetrates 
the  porous  powder  sample  of  density  Pj  to  a  depth  of  Xj,  the 
measured  intensity  is  Ij.  If  the  internal  density  of  the 
sample  is  uniform,  the  characteristic  intensity  can  be  expected 
to  remain  constant  at  1^  even  if  the  bulk  density  of  the  sample 
changes.  To  see  this  behavior,  compare  Figure  4a  with  Figure 
4b.  Although  the  incident  radiation  penetrates  to  a  greater 
depth  in  a  low  density  sample  (X^ ,  the  amount  of  characteristic 
radiation  that  gets  to  the  counter  from  this  depth  is  the  same 
as  that  from  a  smaller  depth  in  a  denser  sample  (X 2 )  .  Thus, 
intensity  is  not  a  function  of  sample  density  if  the  density 
is  uniform. 
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However,  assume  that  some  forming  pressure  produces  the 
density  configuration  shown  in  Figure  4c.  If  the  incident 
radiation  penetrates  to  the  surface  of  the  high  density  region, 
a  greater  amount  of  characteristic  radiation  will  be  produced 
from  this  area  than  the  area  above  it.  The  more  characteristic 
radiation  excited  from  a  given  depth  within  the  sample  the 
greater  the  probability  that  a  given  fraction  will  reach  the 
detector.  The  measured  intensity  will  be  increased  to  ^  at 
the  pressure  that  creates  the  inhomogeneous  regions  of  density. 
Thus,  intensity  can  be  a  function  of  forming  pressure  if 
density  gradients  exist  within  the  sample. 

David  Train  et  al.  ^  has  shown  that  when  powders  of 
magnesium  carbonate  are  pressed  density  gradients  are  introduced 
into  the  sample.  The  density  gradient  distribution  depends  on 
the  pressing  pressure  as  well  as  the  thickness  of  the  sample. 

The  initial  region  of  the  density- intensity  curves  in 
which  di/dp  is  small  can  be  explained  by  the  fact  that  although 
the  bulk  density  is  increasing,  the  local  density  which  is 
exposed  to  incident  radiation,  is  remaining  fairly  constant. 
Therefore,  the  measured  inten.sity  does  not  change. 

The  second  region  in  the  density-intensity  curve  corres¬ 
ponds  to  a  large  di/dp.  Train’s  results  show  that  within  this 
range  of  densities  there  is  a  large  change  in  local  density 
with  only  a  small  change  in  bulk  density.  Therefore,  the  inten¬ 
sity  changes  by  a  large  factor  with  just  a  small  change  in  bulk 
density.  In  other  words,  the  small  change  in  bulk  density  is 
sufficient  to  bring  a  high  density  region  close  enough  to  the 
top  of  the  sample  to  be  exposed  to  the  incident  radiation. 
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The  third  region  of  the  density- intensity  curve  is  a 
plateau.  Train's  results  show  that  within  this  region  there  is 


not  much  change  in  local  density  with  bulk  density.  Thus, 
dl /dp  8  0  (Figure  2) . 

The  x-ray  results  show  that  there  is  also  a  dependence  of 
intensity  on  sample  thickness  for  magnesium  carbonate  powders 
(Figure  3).  Train's  data  are  evidence  that  samples  of  differ¬ 
ent  thicknesses  pressed  to  the  same  density  have  different 
density  gradients.  An  adaptation  of  their  results  is  shown 
in  Figure  5.  The  figure  illustrates  that  the  higher  the  form¬ 
ing  pressure  the  greater  the  changes  in  local  density  with  a 
change  in  thickness.  It  was  seen  that  a  sample  pressed  at 
2,000  psi  and  1.50  cm  thick  gives  a  lower  intensity  than  a 
sample  pressed  at  6,000  psi  and  1.50  cm  thick.  A  sample  pressed 
at  6,000  psi  gives  a  greater  change  in  intensity  if  the  sample 
thickness  is  increased  to  1.75  cm  than  does  a  sample  pressed 
at  2,000  psi  with  a  corresponding  increase  in  thickness.  Conse¬ 
quently,  the  observed  thickness  and  pressure  dependence  of 
di/dp  is  consistent  with  Train's  density  gradient  results. 

The  regions  of  the  intensity-density  curves  for  the  other 
three  compounds  can  be  explained  in  a  similar  manner.  The  curve 
for  titanium  dioxide  resembles  that  of  magnesium  carbonate 
because  titanium  is  a  light  element,  also.  The  curves  for 
chromium  oxide  and  thorium  oxide  show  a  larger  di/dp  in  the 
second  region  than  do  the  lighter  compounds,  because  the  local 
density  in  the  irradiated  volume  changes  more  rapidly  for  the 
heavier  compounds.  The  fact  that  there  is  no  intensity 
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Fig.  5  Density  distribution  in  HgCO-,  sample 
pressed  at  10,000  psi.  ^ 
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dependence  on  thickness  for  titanium  dioxide,  chromium  oxide 
and  thorium  oxide  is  due  to  the  larger  mass  absorption  coeffi¬ 
cients  of  these  compounds  and  that  changes  in  local  density 
gradients  as  a  function  of  thickness  are  not  as  large  as  for 
magnesium  carbonate.  Very  large  length  to  diameter  ratios 
must  be  used  before  thickness  variations  become  important  for 
compounds  containing  heavier  elements. 

Conclusions 

The  results  of  this  paper  indicate  that  in  order  for  an 
accurate  x-ray  spectrochemical  analysis  of  pressed  powder  speci 
mens  to  be  determined,  the  following  rules  should  be  observed: 

1)  The  density  of  the  sample  should  be  carefully 
controlled  for  all  materials  if  the  pressed 
density  lies  in  the  region  where  dl/'dp  is  large. 

The  best  results  are  obtained  when  samples  are 
pressed  in  the  density  range  where  di/dp  -  0. 

This  is  usually  Region  III. 

2)  The  thickness  of  the  sample  is  very  important 
for  compounds  containing  the  lighter  elements. 

All  samples  should  be  the  same  thickness  for 
reproducible  results  to  be  obtained. 

3)  The  sample  thickness  is  not  an  important  facte 
for  compounds  containing  the  heavier  elements 
(Z  >  20)  . 


121 


References 


1.  L.  S.  Birks,  X-ray  Spectrochemical  Analysis.  Interscience  Publishers,  Inc. 
New  York. 

2.  D.  Train,  "Transmission  of  Forces  Through  a  Powder  Mass  During  the  Process 
of  Pelleting,"  Trans.  Instn.  Chem.  Engrs..  Vol.  35,  pp.  258-266,  1957. 


122 


IV.  Mpflanrpmpnt  Tprhnigups  (R.  W.  Gould,  L.  L.  Hench,  and  J.  J.  Hren) 

A.  STRUCTURE  OF  SURFACE  DEFECTS  (D.  L.  Stoltz  and  J.  J.  Hren) 

Introduction 

In  the  past  decade  the  properties  of  solid  surfaces  in  electronic, 
mechanical  and  chemical  phenomena  have  been  gradually  recognized  as  being 
more  and  more  important.  The  influence  of  different  types  of  defects  and 
their  distribution  in  the  bulk  structure  of  materials  has  been  studied  in 
depth  and  utilized  to  explain  various  observed  bulk  properties.  The  advent 
of  ultra-high  vacuum  capabilities  and  sophisticated  observational  methods, 
such  as  transmission  electron  microscopy,  field-ion  and  electron-emission 
microscopy  and  low  energy  electron  diffraction,  now  creates  the  possibility 
to  extend  defect  structural  analysis  to  surfaces  and  to  correlate  this 
analysis  with  the  physical  properties  of  solid  surfaces. 

Ideally,  surfaces  of  crystalline  materials  represent  the  abrupt  ter¬ 
mination  of  the  periodic  array  of  the  crystalline  structure.  Because  of 
the  inherent  complexity  of  real  surfaces,  however,  several  viewpoints 

i* 

should  be  considered  to  gain  a  complete  picture.  For  example,  the  atomic 
arrangement,  electronic  configuration,  and  nature  of  surface  defects  need 
to  be  correlated  with  various  surface-dependent  physical  properties.  Such 
properties  are  directly  dependent  upon  the  nature  of  the  surface.  Such 
properties  include:  electrical  surface  conductance,  electronic  work  func¬ 
tion,  surface  tension,  and  interaction  of  the  surface  with  gas,  liquid  and  solid 
phases.  Pertinent  Phenomena  for  solid-gas  interfaces  include  physical  adsorp- 
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tion,  chemisorption,  and  oxidation  or  film  formation. 
Examples  of  solid-liquid  interfacial  phenomena  are  elec¬ 
trical  double  layers,  electrode  equilibrium  and  kinetics, 
and  corrosion.  Solid-solid  interfaces  are  associated  with 
grain  boundaries,  p-n  junctions,  rectifying  contacts,  and 
dissimilar  solids  in  physical  contact. 

To  analyze  surface  defects  one  must  utilize  techniques 
involving  both  "clean"  and  "real"  surfaces.  Studies  of 
"clean"  surfaces,  which  ideally  represent  an  atomically  flat 
surface  void  of  any  contaminating  atoms,  provide  a  necessary 
means  for  developing  fundamental  concepts  and  the  required 
theoretical  models  to  explain  various  observed  intrinsic 
surface  structures.  Experimental  information  for  "real" 
surfaces,  on  the  other  hand,  will  provide  the  basis  for 
understanding  the  kinetics  and  energetics  involved  in  dif¬ 
ferent  observed  phenomena  and  physical  properties. 

The  importance  of  bulk  defects  is  already  known  for 
semiconductor  devices.  Point  defects,  such  as  impurity 
atoms,  strongly  influence  the  electrical  properties  of  the 
devices,  particularly  the  minority  carrier  lifetime.  In 
addition,  dislocations  have  been  shown  to  increase  the 
diffusion  rate  of  both  acceptor  and  donor  atoms  as  compared 
to  the  perfect  material.  Queisser  e  t  a  1 .  ,  ^  for  example, 
have  attributed  this  rate  increase  to  the  formation  of 
dopant  atmospheres  at  dislocations,  enhanced  diffusion  by 
vacancy  mechanisms,  and  lowered  activation  energies.  In 
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turn,  this  enhanced  mass  transport  disrupts  the  sharp  im¬ 
purity  concentration  gradients  within  the  devices  by  the 
build-up  of  impurity  atoms  into  dislocation  atmospheres 
and,  sometimes,  in  the  formation  of  precipitates.  These 
structural  changes  then  create  short  circuits  between  the 
emitter  and  collector  regions  and  cause  reduced  device 
lifetime.  A  similar  structural  change,  involving  precipi¬ 
tates  of  heavy  metals  in  p-n  junctions,  results  in  raising 
the  reverse  currents  (called  softening)  to  such  a  degree 
that  drastic  failures  occur.  Several  workers'-  1  have  men¬ 
tioned  that  these  precipitates  nucleate  at  dislocations 
and  second-order  twin  boundaries  and  report  that  diodes 
with  high  dislocation  contents  and  second-order  boundaries 
were  soft,  while  diodes  with  lower  line  defect  densities 

f  31 

showed  desirable  reverse  current-voltage  characteristics.  1 
It  seems  apparent  that  defect-free  materials  would  be 
desirable  for  semiconductor  devices.  However,  in  1960 
empirical  tests  of  transistors  and  diodes  made  from  dislo¬ 
cation-free  silicon  indicated  no  superiority  compared  with 
devices  made  from  materials  with  moderate  dislocation 
densities.  1  It  was  interesting  that  subsequent  investi¬ 
gations  revealed  that  the  reason  for  this  (at  first  blush) 
surprising  lack  of  superiority  of  the  perfect  material  was 
simply  that  the  material  was  not  free  from  defects  follow¬ 
ing  the  manufacturing  processes. 

Briefly,  one  can  state  various  possible  and  recognized 
means  of  introducing  defects  in  solid  state  electronic 


125 


devices . 


These  include: 


1.  variation  in  dopant  concentrations  in  single 
crystals  which  produce  spiral-like  patterns  of 
impurity  atoms.  Such  spirals  of  impurity  atoms, 
according  to  Goetzberger , ^  result  in  diodes 
exhibiting  a  striated  surface  with  low 
breakdown  voltages. 

2.  introduction  of  dislocations  during  slicing  of 
single  crystals  and  mishandling  prior  to  heat 
treatment. 

3.  nucleation  and  growth  of  precipitates  of  im¬ 
purity  atoms  during  heat  treatment. 

4.  introduction  of  interstitial  point  defects  at 
the  surface  of  the  single  crystal  by  diffusion 
of  impurity  atoms. 

5.  coalescence  of  point  defects  and  rearrangement 
of  dislocation  arrays  during  application  of 
high  electric  fields. 


Methods  of  Observation 


From  the  preceding  one  recognizes  that  the  bulk  and 
surface  defects  in  materials  play  a  major  role  in  determining 
the  materials  performance  in  various  applications.  It  seems 
appropriate,  therefore,  to  briefly  describe  various  means  of 
observing  surface  and  bulk  defects. 
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Several  optical  methods  exist  for  revealing  the  pres¬ 
ence  of  dislocations  on  the  surface  and  interior  of  the 
crystal.  Etching  techniques  expose  dislocations  at  the  sur¬ 
face  by  forming  pits  at  their  points  of  emergence  with  the 
surface.  Initiation  of  these  pits  appears  to  be  related  to 

impurity  atmospheres  and  local  lattice  strain  associated 
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with  dislocations  intersecting  the  surface.  1  This  technique 
is  particularly  valuable  for  determining  low  dislocation 
densities  and  observing  general  dislocation  arrangements. 
Another  method  of  optical  observation  of  dislocations  is 
the  decoration  technique.  Wherein  one  intentionally 
causes  the  separation  of  a  second  phase  to  occur  preferen¬ 
tially  on  the  dislocations.  This  technique  permits  surface 
and  volume  arrays  to  be  observed  optically.  Figure  1  repre¬ 
sents  a  dislocation  network  of  mixed  type  in  KC1,  decorated 
with  particles  of  Ag.  Dash^has  also  observed  copper 
decorated  dislocations  in  silicon  by  using  infra-red  radia¬ 
tion.  However,  the  decoration  technique  is  clearly  restric¬ 
tive,  since  the  thermal  diffusion  treatments  required  to 
decorate  the  dislocations  can  modify  the  dislocation  arrange¬ 
ment  and  in  any  event  leave  it  permanently  altered.  The 
general  phenomenon  of  photoelasticity  can  be  also  used  to 
reveal  dislocations,  since  the  lattice  strain  around  a  dis¬ 
location  causes  an  isotropic  crystal  to  become  birefringent . 
Kear  and  Pratt, v  '  for  example,  used  a  polarizing  microscope 
to  study  quenching  strains  in  rock  salt,  so  that  the  colored 
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Fig.  1.  Network  of  mixed  dislocations  in  KC1,  decorated 
with  silver  particles.  (Amelinckx,  S.,  Acta 
Met.,  6,  34  (1958).) 
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birefringent  bands  were  sharply  visible  and  could  be  related 
to  the  different  signs  of  the  dislocations.  Indenbohm^-1^ 
utilized  all  three  of  the  above  optical  techniques  on  the 
same  sample  and  obtained  a  perfect  correlation. 

Due  to  the  opaqueness  and  high  dislocation  content  of 
many  materials,  however,  other  observational  techniques  have 
been  developed  that  utilize  various  higher  energy  incident 
beams.  One  such  method  is  the  transmission  electron  micro¬ 
scopy.  Hirsch^11)  and  Bollman,^1^  independently,  first 
obtained  high  magnification  images  of  line  defects  in  thin 
foil  sections  of  aluminum  and  stainless  steel  respectively. 
These  dislocations  are  made  visible  in  the  electron  micro¬ 
scope  by  means  of  diffraction  contrast  arising  from  the 
strain  field  associated  with  the  dislocations.  Thus,  the 
trace  of  the  dislocation  line  itself  is  observed.  This 
well  known  technique  has  proven  itself  very  successful  in 
observing  and  studying  dislocations,  antiphase  boundaries, 
point  defect  aggregates,  dislocation  loops,  stacking  faults, 
dislocation  interactions  and  arrays,  second  phase  particles, 
and  Guinier-Preston  zones.  Overall,  this  method  has  prob¬ 
ably  contributed  more  to  our  understanding  of  bulk  crystalline 
defects  than  any  other,  but  we  will  not  comment  further  on 
it  here  because  of  the  already  extensive  literature  extant. 

A  monochromatic  x-ray  beam  will  also  pass 
through  a  crystal  and  yield  contrast  effects  very  similar 
to  transmission  electron  microscopy.  This  observational 
method  is  particularly  suited  to  materials  with  low  defect 
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density  and  has  been  used  increasingly  in  recent  years.  Two 
such  methods  are  the  Berg-Barrett ^ and  Lang^**^ 
techniques.  Of  these  the  Lang  method  appears  most  promising 
for  producing  topographs  from  crystals  with  optically  flat 
"real"  surfaces.  These  projection  topographs  represent 
dark  field  images  and  the  diffracting  conditions  are  readily 
controlled  and  defined.  One  application  of  this  technique 
is  the  determination  of  Burgers  vectors  from  paired  projec¬ 
tion  topographs . Figure  2  is  an  example  and  shows  an 
enlarged  section  of  a  (224)  topograph  of  a  silicon  wafer 
after  boron  diffusion.  Notice  the  dark  dislocations  (thin 
bands)  and  precipitates  (black  spots) .  Stereographic  pairs 
of  such  topographs  (using  (hkl)  and  (hkl)  reflections) 
create  a  three  dimensional  representation  of  any  dislocation 
structure. 

All  of  the  observational  techniques  just  described 
have  very  limited  utility  in  determining  surface  defect 
structures.  The  remaining  discussion  will  therefore  be 
primarily  concerned  with  the  use  of  the  field-ion  microscope 
which  is  probably  most  applicable  for  such  studies.  The 
method  was  initially  developed  by  Muller as  a  modi¬ 
fied  version  of  the  electron-emission  microscope .  A 

schematic  diagram  of  the  instrument  is  shown  in  Figure  3. 

The  specimen  is  generally  a  small  diameter  wire  which  has 
been  electrochemically  polished  to  a  very  sharp  point 
('1000  A) ,  and  field  evaporated  to  produce  a  nearly  atomically 
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Fig.  2.  224  topograph  of  silicon  after  boron  diffusion. 

(Schwuttke,  G.  H.,  AFCRL-66- 192 ,  p.  6-6  (1966).) 


© 
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Fig.  3.  Schematic  diagram  of  field- ion  microscope. 
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smooth,  hemispherical  surface.  An  imaging  gas,  usually 
helium,  is  introduced  into  the  chamber  and  its  atoms  are 
ionized  at  the  vicinity  of  the  discrete  atomic  surface. 

These  positive  ions  are  repelled  by  an  accelerating  electri¬ 
cal  field  of  approximately  450  MV/cm  to  the  phosphor  coated 
screen  at  ground  potential. 

The  image  formed  by  the  field-ion  microscope  repre¬ 
sents  an  enlarged  projection  of  the  detailed  field  distribu¬ 
tion  of  the  surface  and  depends  sensitively  upon  local  field 
strengths.  The  atomic  structure  of  the  surface  greatly 
influences  this  field  distribution,  i.e.,  the  protruding 
atoms  create  high  local  fields  which  in  turn  produce  bright 
image  points  as  shown  in  Figure  4.  Note  the  numerous  small 
circular  facets  present  and  that  the  atoms  at  the  edges  of 
these  facets  appear  brightest  because  their  protrusion  leads 
to  a  high  localized  field  strength.  The  symmetry  of  this 
image  allows  the  various  crystal  faces  to  be  identified 
crystallographically.  Figure  5  illustrates  that  the  path 
traveled  by  the  ions  is  approximately  normal  to  the  local 
tip  surface.  The  geometrical  magnification  is  normally 
about  10**,  as  determined  by  the  ratio  of  the  distance  from 
tip  to  screen  to  the  radius  of  the  tip. 

It  should  be  added  that  field- ion  microscopy  not  only 
allows  individual  atomic  observation  of  a  surface,  but  by 
use  of  field  evaporation,  individual  atomic  layers  can  be 
removed  to  provide  an  atom  by  atom  reconstruction  of  a 
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Fig.  4.  Field-ion  micrograph  of  platinum. 
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Fig.  5.  Schematic  illustration  of  path  of  ionized  gas 
from  specimen  tip  to  phosphorus  screen. 
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volume.  Also,  by  using  field  evaporation  in  ultra-high 
vacuums  (better  than  10  ^  Torr) ,  "clean"  surfaces  can  be 
formed  and  retained  for  extended  periods. 


Interpretation  of  Surface  Defects 
by  Field-Ion  Microscopy 

Despite  the  various  experimental  difficulties  in 
observing  defects  in  the  small  hemispherical  surfaces  of  the 
field-ion  specimens,  a  number  of  researchers  (of  Ref.  (21)) 
have  obtained  images  of  dislocations,  grain  boundaries, 
stacking  faults  and  point  defects  and  analyzed  these  in 
depth.  We  will  describe  several  of  these  in  detail. 

One  of  the  most  commonly  appearing  defects  detected 
in  field-ion  images  is  the  spiral  structure.  Most  often  in 
discussions  of  crystal  growth,  these  spirals  have  been  des¬ 
cribed  as  resulting  from  the  intersection  of  a  screw  dis¬ 
location  with  the  surface.  In  fact,  however,  if  one 
scrutinizes  the  basic  literature  on  the  geometry  of  dislo¬ 
cations,  it  becomes  clear  that  any  dislocation  can  give 
rise  to  an  unending  step  or  spiral  ramp  when  it  intersects 
a  surface.  For  example,  Cottrell^^)  pointed  out: 

"Suppose  that  we  have  a  family  of  parallel 
planes  and  project  through  them  a  dislocation 
line  with  a  suitable  Burgers  vector  which  also 
projects  through  them.  Then  the  structure  is 
changed  by  the  dislocation  from  a  family  of 
parallel  planes  into  a  single  spiral  surface 
or  "helocoid."  This  property  is  true  of  edge, 
screw  and  mixed  dislocations,  but  can  be  seen 
most  easily  in  the  case  of  screw  dislocations 
which  project  vertically  through  the  family  of 
planes . " 
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Papers  by  George ^23^  and  Ranganathan  et  al .  ,  ^2^  that 

explained  spiral-like  structures  at  grain  boundaries  and 

streaks  in  field-ion  images  respectively,  were  close  in 
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recognizing  the  correct  interpretation.  Pashley^  1  pre¬ 
sented  a  quantitative  partial  interpretation  of  dislocation 
images  that  included  a  mathematical  description  of  the  pitch 
of  the  helix  structure  about  any  dislocation  line  for  which 
fi  •  S  f  0,  where  fi  is  a  unit  vector  normal  to  the  (hkl) 
planes  and  S  is  the  Burgers  vector.  He  stated:  "It  is 
important  to  realize  that  the  effect  of  dislocations  on  the 
field- ion  micrograph  depends  only  on  the  component  of  the 
Burgers  vector  perpendicular  to  the  (hkl)  planes  concerned. 
Thus,  a  dislocation  showing  a  spiral  structure  is  not 
necessarily  a  pure  screw  dislocation..."  Ranganathan^26-^ 
independently  developed  these  ideas  while  attempting  to 
explain  grain  boundary  structures  in  terms  of  dislocation 
networks.  His  work  is  more  detailed  than  Pashley's,  but  is 
quantitatively  the  same.  Ranganthan's  terminology  deduced 
the  nature  of  a  spiral  from  the  dot  product  gj^j  •  ft  in 
analogy  to  electron  microscope  interpretations,  where  the 
reciprocal  lattice  vector  is  equivalent  to 

ghki  *  h  a‘  +  k  S*  H  c*  (1) 

By  definition  ghkl  is  perpendicular  to  (hkl)  and  has  magni¬ 
tude  irrespective  of  crystal  structure.  The  Burgers 

vector,  ft,  is  a  real  lattice  vector  and  proceeds  from  one 
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lattice  point  to  another,  if  it  is  a  total  dislocation.  It 
need  not  be  (and  normally  is  not)  coincident  with  one  of 
the  basis  vectors  of  the  lattice,  but  it  may  always  be  written 
in  terms  of  a  constant  K  and  the  real  basis  vectors  a,  £,  c 
as : 

S  =  K(u  a  +  v  S  +  w  c)  (2) 

The  dot  product  is  then  the  following: 

ghki  -  6  =  Ch  a*  *  k  S*  +  1  S*)  •  (u  a  *  v  «  *  w  c)K 
=  K(h  u  +  k  v  +  1  w)  (3) 


where  the  quantity  in  brackets  is  clearly  either  zero  or 
an  integer.  The  constant  K,  therefore,  determines  whether 
the  dot  product  is  an  integer  or  not.  If  t  is  a  total  dis¬ 
location,  K  is  either  zero  or  an  integer  since  S  goes  from 
one  lattice  point  to  another.  On  the  other  hand,  if  ft  is  a 
partial  dislocation  (superlattice  dislocations  may  be  in¬ 
cluded  in  this  definition) ,  K  may  be  such  as  to  give  a 
fraction  for  the  dot  product.  The  gj^  *  S  criterion  can 
be  illustrated  by  considering  the  real  lattice  vectors  of 
a  cubic  lattice.  The  reciprocal  lattice  vector  may  be 

expressed  as: 


t  *  fihkl 
hkl  dhkl 


(4) 


where  is  unit  normal  vector  to  (hkl)  and  d 

planar  spacing.  Using  the  following  identities 


hkl 


is  inter- 
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and 


fihkl 

dhkl 


Thkll 

'.2  i — 7J7T7J 

(h  *  k  +  1  )  ' 


a. 


(h2  +  k2  +  12)1/2 


(5) 

(6) 


where  aQ  is  the  unit  cell  dimension,  equation  (4)  becomes 


Ihki  *  [hkl]/a0 


(7) 


Consider  a  Burgers  vector  in  the  face  centered  cubic  struc¬ 
ture,  S  B  aQ/2  <110>,  then  by  substituting  into  equation 
(4) ,  one  obtains 

Ihkl  •  *  *  7  th  ♦  k]  (8) 

or  the  equivalent  expression  depending  on  the  exact  Burgers 
vector  value.  Since  the  possible  (hkl)  planes  in  face 
centered  cubic  structures  have  all  odd  or  all  even  indices, 
the  last  equation  is  always  zero  or  an  integer.  The  quan¬ 
tity  *  b  may  be,  thus,  thought  of  as  being  the  projec¬ 

tion  of  the  Burgers  vector,  5,  on  the  unit  vector  ft  in 
units  of  d^j. 

An  alternate  way  to  state  this,  due  to  Smith  et_al. ,(27) 
is  as  follows:  "Characterize  a  surface  facet  of  a  field- 
ion  specimen  by  its  normal  fthkl.  Then  the  Burgers  vector 
of  a  dislocation  is  given  by  the  closure  failure  in  the 
Burgers  circuit. ^28^  But  this  is  independent  of  the  circuit 
selected,  as  long  as  it  is  always  taken  in  good  crystal." 
Consider  the  schematic  description  of  this  statement  in 
Figure  6.  One  notes  that  when  the  Burgers  vector  lies  in 
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Fig.  6.  Pure  edge  dislocation's  intersection  with 
rational  surfaces,  illustrating  creation 
of  spiral  surface. 
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the  plane  of  the  surface,  it  is  defined  by  a  Burgers  circuit 
lying  in  that  surface  also.  If  the  Burgers  circuit  is 


taken  for  this  identical  edge  dislocation  on  an  arbitrary 
plane  surface,  it  will  be  converted  into  a  spiral  surface 
to  retain  its  original  direction  and  magnitude.  Utilizing 
the  *  &  expression  to  define  the  pitch  of  the  disloca¬ 

tions  spiral,  where  t  is  equal  to  the  lattice  vector  and 
the  represents  a  rational  plane,  this  dot  product  must 

be  an  integral  number  of  interplanar  spacings. 

The  conclusion  to  be  emphasized  is  that,  except  for 
the  relatively  rare  case  where  the  Burgers  vector  of  a 
dislocation  intersecting  the  surface  lies  in  the  plane  of 
the  surface,  all  dislocations  will  produce  an  unending 


surface  ramp  which  extends  over  large  regions  of  the  sur¬ 


face.  Thus,  if  the  surface  is  truly  atomically  flat,  the 
step  would  proceed  from  the  emergence  point  of  the  disloca¬ 
tion  line  to  the  edge  of  the  piece,  as  shown  in  Figure  6. 

If  the  surface  has  a  long  range  curvature,  such  as  the 
hemispherical  field  emitter  surface,  the  step  adopts  the 
character  of  a  spiral  when  viewed  from  above.  Figure  7 
illustrates  that  the  sense  of  the  spiral  can  be  established 
by  the  sense  of  the  Burgers  vector  with  respect  to  the 
(hkl)  planes  of  interest,  i.e.,  when  the  Burgers  vector 
points  inward  the  spiral  has  one  sense  and  when  it  points 
outward  the  spiral  has  the  opposite  sense. 

Since  the  character  of  the  dislocation  is  not  required 
to  predict  the  nature  of  the  expected  spiral  structure, 
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Fig.  7.  Schematic  drawing  indicating  relationship 
between  sense  of  spiral  formed  on  surface 
of  hemisphere  when  intersected  by  a  disloca¬ 
tion  with  the  direction  of  the  normal  com¬ 
ponent  of  the  Burgers  vector. 
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computer  models  can  be  <*"vised,  based  only  upon  the  Burgers 
vector  and  the  plane  intersected  by  the  dislocation  to  test 
the  predictions  of  the  g^j  *  &  criterion.  Utilizing 
arbitrary  dislocations  and  available  mathematical  descrip¬ 
tions  of  the  atomic  displacements  around  the  dislocation, 
the  geometry  of  the  unending  ramp  may  be  accurately  simu¬ 
lated  in  a  digital  computer.  Sanwald  and  co-authors^®) 
first  selected  a  pure  screw  dislocation,  defined  by  a 
Burgers  vector,  S  ■  ^  {110],  intersecting  a  family  of  the 
{420}  type  planes.  The  (204)  pole  is  reproduced  in  Figure 
8  and  illustrates  the  single  spiral  from  the  g^j  •  S 
criterion.  Figure  9  corresponds  to  an  experimental  field- 
ion  image  of  a  single  spiral  and  verifies  the  geometrical 
foundations  of  the  spiral  hypothesis. 

A  more  complicated  spiral  structure  was  developed  by 
assuming  two  closely  spaced  dislocations  with  aQ/2  <110> 
Burgers  vectors  intersecting  the  image  near  the  (113)  planes 
of  a  FCC  lattice.  By  selecting  for  this  computer  simula¬ 
tion  the  Burgers  vectors  ■  aQ/2  [011}  and  S2  *  a0/2  [011} 
and  assuming  them  to  be  separated  by  10  A  in  the  specimen, 
the  resulting  image  shown  in  Figure  10  is  obtained.  The 
experimental  image  of  iridium  near  the  (113)  plane,  found 
in  Figure  11,  shows  good  qualitative  agreement  with  the 
computer  simulated  image. 

When  the  dislocation  intersection  does  not  emerge 
from  the  center  of  the  pole  (Figure  8) ,  one  obtains  images 
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(204)  POLE 


Fig.  8.  Computer  simulated  image  of  pure  screw 

dislocation  emerging  at  center  of  the  (204) 
pole  to  produce  single  leaved  spiral. 
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Fig.  9.  Experimental  field-ion  micrograph  of  an  Ir 
specimen  containing  a  single  spiral. 
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Fig.  10. 


Computer  simulated  image  of  triple  spiral 
caused  by  pair  of  noninteracting  dislocations. 
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Fig.  11.  Triple  spiral  on  (113)  plane  of  an  Ir  specimen. 
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similar  to  Figure  12.  Figure  12a  represents  a  pure  screw 
dislocation  causing  a  single  spiral  on  the  (204)  plane. 

Note  that  the  only  difference  between  Figures  8  and  12a  is 
the  point  of  emergence  of  the  pure  screw  dislocation.  When 
analyzed  in  terms  of  the  expected  single  spiral  structure , 
the  result  is  obviously  correct,  yet  if  such  a  structure 
was  observed  experimentally  it  could  be  easily  misinterpreted 
as  showing  an  extra  half-plane  and  be  mistaken  for  an  edge 
dislocation. 

Figure  12b,  on  the  other  hand,  illustrates  the  simu- 
lated  image  of  a  pure  edge  dislocation  in  the  same  location 
of  the  crystal  facet  as  the  pure  screw  dislocation  in 
Figure  12a.  There  is  little  difference  in  the  appearance 
of  the  two  images.  The  same  single  spiral  structure  is 
predicted  in  both  cases  and  only  the  diameters  of  the  rings 
appear  to  be  different  (in  the  computed  images).  This 
latter  effect  is  due  simply  to  the  sensitivity  of  the  com¬ 
puter  model  to  the  exact  tip  radius, which  was  slightly 
different  in  the  neighborhood  of  the  pole  in  each  case. 
Figures  12a  and  12b  point  up  the  need  for  careful  interpre¬ 
tation  of  experimental  field-ion  images  to  accurately 
characterize  the  spiral  structures. 

So  far  attention  has  been  directed  primarily  to  the 
interpretation  of  surface  steps  caused  by  total  dislocations. 
Straightforward  application  of  the  *  ft  criterion  to 

partial  dislocations  of  the  type  aQ/6  <112>  and  stacking 
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Fig.  12.  Pure  screw  and  edge  dislocations  causing 
single  spiral  on  (204)  plane  edges. 

(*  -  point  of  emergence  of  dislocation) 


faults  will  be  illustrated  in  the  FCC  structure.  Although 
each  partial  dislocation  will  not  yield  a  step  with  an  inte¬ 
gral  number  of  d-spacings,  a  pair  of  partial  dislocations 
together  yields  a  whole  dislocation  with  a  Burgers  vector 
from  one  lattice  point  to  another.  (This  argument  also 
applies  for  super  dislocations  in  ordered  alloys.)  Conse¬ 
quently,  the  long-range  effect  must  give  a  step  with  an 
integral  number  of  d-spacings,  as  illustrated  previously 
for  total  dislocations.  The  stacking  fault,  between  the 
pair  of  partial  dislocations,  will  contribute  a  fractional 
step.  Figure  13  shows  a  pair  of  Shockley  partial  disloca¬ 
tions,  aQ/6  <1 1 2> ,  intersecting  a  (220)  plane.  The  total 
5  is  j  [101].  The  displacement  of  the  rings  along  a  great 
circle  joining  the  two  dislocations  corresponds  to  the 
intersection  of  the  stacking  fault  plane  with  the  surface. 
Since  the  two  Shockley  partial  dislocations  were  separated 
by  a  small  distance,  this  disturbance  appears  along  a 
straight  line  joining  the  two  emergence  points  of  the  dis¬ 
locations.  Away  from  each  of  the  dislocations  one  notices 
a  single  spiral  which  is  expected  from  the  total  dislocation 
from  which  the  pair  has  dissociated.  Experimental  confirm¬ 
ation  of  this  analysis  comes  from  a  recent  study  by  Newman 
et  al.,^31’3^  which  revealed  a  Frank  loop  and  its  associ¬ 
ated  stacking  fault.  The  stacking  fault  can  be  seen  as  a 
streak  connecting  the  two  spirals  of  opposite  sense  in 
Figure  14.  This  configuration  was  simulated  as  a  Frank 


150 


Stocking  Fault  on  (220)  Pole 


Fig.  13.  Pair  of  Shockley  partial  dislocations  and 

associated  st|cking  fault  intersecting  (220) 
plane;  total  b  =  a/2  [101].  (*  =  point  of 

emergence  of  dislocation) 
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Fig.  14.  Computer  simulation  of  Frank  loop  with 
associated  stacking  fault. 
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loop  on  the  computer  and  the  result  is  shown  in  Figure  15. 
The  agreement  is  clearly  very  good. 

Since  all  dislocations  which  have  a  component  of  their 
Burgers  vectors  normal  to  the  surface  will  produce  steps  if 
they  intersect  it,  it  is  instructive  to  calculate  the  amount 
of  new  surface  created  thereby.  This  will  also  serve  to 
illustrate  the  potential  importance  of  surface  defects  on 
physical  and  chemical  properties.  Using  an  extension  of 
the  g  *S  criterion,  it  is  possible  from  pure  geometry  to 
write  the  expression: 


*, r  ■  *0  -W  pbt  l£l*  c» 

where  A  =  surface  area  created  by  the  dislocation 
®hkl  5  reciprocal  lattice  vector  of  (hkl) 

.  =  Burgers  vector  of  i**1  dislocation 

dj^  =  interplanar  distance  (hkl) 

5  density  of  dislocations  with  Burgers  vector  b^ 

| L |  =  parameter  which  corresponds  to  average  step 
length  from  point  of  emergence  of  any  dislo¬ 
cation  to  edge  of  planar  area* 

An  order  of  magnitude  for  A  may  be  obtained  from  Equation 
(9)  by  assuming  the  following  values: 


Pbi  »  10®  cm/cm3 

«hki  "  I111! 

•  a/2  [110] 
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Fig.  15.  Extreme  blow-up  of  region  containing  a  Frank 
loop  with  associated  stacking  fault.  Note 
the  streak  connecting  the  two  spirals  of 
opposite  sense. 
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*  o 

dhki  *  2  x  10~°  cm 
|£|  *  1/2  cm 

2 

The  resulting  value  of  As,  created  by  dislocations,  is  1  cm  , 
the  same  as  the  nominal  initial  surface  area  of  the  sample  1 
This  amount  of  new  surface  area  must  markedly  affect  the 
properties  of  the  material's  surface.  For  instance,  in 
Metallic  samples  these  surface  steps  could  develop  electric 
double  layers,  whereas  in  ionic  materials  localized  dipole 
regions  could  be  formed.  For  covalently  bonded  materials, 

such  as  silicon  and  germanium,  these  steps  could  create  a 

« 

significant  increase  in  the  number  of  dangling  bonds.  It 
follows  that  surface  step  structures  from  dislocation 
intersections  must  markedly  affect  the  electronic,  chemical 
and  mechanical  properties  of  real  surfaces  in  general. 

Point  defects  also  may  be  observed  and  analyzed  by 
field-ion  microscopy.  Use  of  controlled  field  evaporation 
permits  identification  of  single  lattice  sites  on  the 
specimen's  surface.  Since  field  evaporation  removes 

atoms  almost  exclusively  from  edge  sites,  the  appearance 
of  a  void  in  the  interior  of  a  fully  resolved  atomic  surface 
corresponds  to  a  vacancy.  Therefore,  the  concentration 
and  distribution  of  vacancies  can  be  determined  from  the 
field-ion  micrographs.  Point  defects  may  affect  many 
other  properties  of  the  surface  as  well.  For  example, 
vacancies  and  aggregates  of  vacancies  on  the 
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surfaces  of  ultra-high  vacuum -cleaved  rock  salts  contribute 
significantly  to  the  nucleation  and  growth  of  single  crystal 
thin  films.  Rhodin  et  al . have  reported  recently  that 
surface  vacancies  created  by  electron  irradiation  provide 
nucleation  sites  favoring  the  formation  of  oriented  nuclei 
for  an  epitaxial  layer  of  gold,  silver,  or  copper  on  KC1. 
Another  recent  study  of  epitaxial  deposition  of  platinum 
on  iridium  was  made  by  Graham  et  al.  using  the  field- 
ion  microscope.  Their  results  indicate  that  the  incident 
platinum  atoms  preferred  to  be  trapped  selectively  in  the 
most  energetically  stable  positions  on  the  surfaces,  i.e., 
those  positions  representing  vacancies,  not  depressions, 
in  the  iridium  surface. 

In  similar  ways,  defects,  other  than  those  described, 
can  be  quantitatively  studied  via  fie Id- ion  microscopy  and 
the  image  interpretation  aided  by  computer  simulated 
models.  These  include,  for  example,  grain  boundaries  and 
anti-phase  boundaries  in  ordered  alloys.  Some  studies 
have  been  made,  but  many  more  remain  to  be  carried  out. 

Sunmary  and  Conclusion! 

Although  the  observational  techniques  and  identifica¬ 
tion  procedures  of  various  surface  defects  have  been  dis¬ 
cussed,  the  actual  role  surface  defects  play  in  affecting 
electronic  and  chemical  properties  still  needs  to  be 
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evaluated  in  detail.  A  more  complete  understanding  of  the 
complex  nature  of  surface  defects  requires  an  integrated 
investigation  involving  not  only  the  characterization  of 
the  surface  and  its  defects  but  also  the  correlated  deter¬ 
mination  of  particularly  important  physical  properties  as 
well.  The  experimental  techniques  would  certainly  include: 
field  emission,  low-energy  electron  diffraction,  Auger 
spectroscopy,  photo-electric  emission  and  optical  proper¬ 
ties.  Such  studies,  for  instance,  would  involve  determina¬ 
tions  of  surface  conductivity  and  electronic  work  functions 
for  similar  specimens  with  varying  defect  content  and  type. 
By  controlling  the  environment  and  surface  structure 
measurements  could  be  made  on  "clean’'  and  "real"  surfaces 
of  the  same  bulk  material.  Another  possibility  would 
include  evaluation  of  a  material's  chemical  reactivity  as 
a  function  of  its  surface  area  and  structure,  as  in  cataly¬ 
sis,  oxidation,  or  adsorption.  Similarly,  further  studies 
involving  the  nucleation  and  growth  of  thin  films  (or 
epitaxial  layers)  on  various  substrates  must  be  conducted. 
Various  means  of  computer  aided  calculation  of  both  the 
geometric  and  electronic  properties  of  surfaces  exist  and 
should  be  employed  with  all  of  the  experiments  suggested. 
Suitable  tools  are  now  at  hand  and,  with  perhaps  slight 
modifications,  we  should  employ  them  actively  to  fulfill 
a  long  vacant  information  gap. 
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B.  THE  EVOLUTION  OF  PARTICLE  SIZE  DISTRIBUTION  FROM  DATA  ON  THE  LOW 
TEMPERATURE  AGING  OF  ALUMINUM-ZINC  ALLOYS  (R.  W.  Gould) 

Abstract 

In  this  paper,  a  method  is  used  by  which  the  particle  size 
distribution  can  be  determined  as  a  function  of  aging  time  from 
small-angle  x-ray  data  on  the  low  temperature  aging  (-4S°C)  of 
three  aluminum- zinc  alloys.  A  series  of  particle  size  distribu¬ 
tion  curves  and  particle  diameter  growth  paths  results  from  this 
analysis,  and  the  evolution  of  the  particle  size  distribution  can 
be  followed. 

Introduction 


During  the  past  fifteen  years,  small-angle  x-ray  scattering 
and  transmission  electron  microscopy  have  been  used  to  study  the 
growth  of  Guinier-Preston  zones  in  metallic  solid  solutions.  The 
aluminum- zinc  system  has  been  of  particular  interest  as  it  contains 
spherical,  zinc-rich  Guinier-Preston  zones  which  give  rise  to  high 
mass  contrast  in  the  transmission  electron  microscope  and  which 
produce  easily  observable  small-angle  scattering  effects.  When 
small-angle  scattering  was  first  applied  to  this  problem,  one  of 
three  size  parameters  was  generally  used  to  characterize  the  growth 
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of  the  zones,  Guinier  radius,  Porod  or  surface  area  "radius,"  and 
a  radius  determined  from  the  position  of  the  small-angle  scatter¬ 
ing  intensity  maximum.  These  single  parameters  were  difficult  to 
correlate  with  the  transmission  electron  microscopy  measurements 
which  became  available  in  the  mid-1950s  as  shown  clearly  by  Baur 
and  Gerold.  Harkness,  Gould  and  Hren,  in  a  recent  publication 
have  shown  that  the  Guinier  radius  alone  does  not  give  a  good  mea¬ 
sure  of  the  mean  zone  radius.  But  when  the  Guinier  radius  and  the 
Porod  radius  are  both  available,  the  combination  of  these  two  mea¬ 
surements  enables  the  particle  size  distribution  to  be  determined 
if  the  particles  are  spherical  and  a  size  distribution  function 
can  be  assumed. 

For  the  present  study,  the  system  a' uminum- zinc  was  chosen 
and  data  was  taken  from  a  publication  by  Bonfiglioli  and  Guinier3 
in  which  the  aluminum- zinc  alloys  were  aged  at  -45°C  in  order  to 
study  the  very  early  stages  of  zone  growth.  Three  alloys  were  ex¬ 
amined  by  Bonfiglioli  and  Guinier,  hereafter  referred  to  as  B  and 
G,  i.e.,  10,  20  and  30  wt.  I  zinc  in  solid  solution  in  aluminum. 

Analysis  of  Data 

B  and  G  studied  the  structure  of  the  Guinier-Preston  zones  in 
the  first  stages  of  their  formation,  thus  requiring  the  low  temp¬ 
erature  of  aging.  During  the  course  of  their  investigation,  they 
measured,  as  a  function  of  aging  time,  several  parameters  charac¬ 
teristic  of  G.P.  zone  formation  which  are  listed  in  Table  I.  Table 
I  shows  the  time  rate  of  change  of  the  Guinier  radius  (Rq)  ^  and 
the  Porod  radius  (Rp)  as  well  as  the  integrated  intensity  Q(0) . 
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TABLE  I 


3 

Zonal  State  Parameters  of  Three  Al-Zn  Alloys  Aged  at  -45°C. 


Alloy 

Aging 

Time 

Minutes 

Volume  Fraction 
of  Zones 

Q(0)  =  2tt  r  sj(s)ds 

XlV 

rg 

0 

(A) 

Rp 

o 

CA) 

Al-10  w/o 

Zn 

10' 

0.0318 

2.07 

6.34 

5.40 

30 

2.53 

7.34 

5.44 

90 

3.19 

9.48 

6.33 

270 

3.35 

10.7 

7.28 

810 

3.58 

11.5 

7.30 

2430 

3.69 

11.5 

7.26 

Al-20  w/o 

Zn 

0 

0.0889 

3.10 

5.0 

3.10 

10 

6.50 

13.6 

7.58 

30 

6.90 

14.9 

12.2 

90 

8.35 

15.0 

11.8 

270 

8.35 

15.3 

12.1 

810 

8.72 

15.4 

13.9 

2430 

8.57 

16.4 

15.2 

Al-30  w/o 

Zn 

10 

0.135 

10.65 

14.6 

9.62 

30 

11.88 

15.2 

11.7 

90 

11.75 

16.1 

14.1 

270 

12.15 

16.5 

13.9 

810 

12.10 

16.3 

13.9 

2430 

11.90 

16.1 

14.2 
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Particle  size  distribution  curves  for  Al-10%  Zn  after  various 
aging  times  at  -45°C. 


Fig. 


a  2  N 

»iO|x*ujo/S3NOZ 


Particle  size  distribution  curves 
times  at  -45*C. 
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0.30 


20%  Zn  after  various  aging 


A!- 10%  Zn  T*»  — 45*C 
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ZONE  RADIUS  A 


In  order  to  measure  the  Porod  radius  it  is  considered  good  experi- 

3  3  3 

mental  procedure  to  show  that  a  plot  of  s  J(s)  versus  s  results 

in  a  straight  line  of  slope  zero  over  a  rather  large  angular  range 
which  includes  the  Porod  region  of  the  small-angle  scattering  curve. 
Such  a  straight  line  plot  insures  that  the  Porod  law  is  valid,  i.e., 
that  there  is  a  sharp  interface  between  the  zone  and  the  surround¬ 
ing  matrix.  B  and  G  in  their  Figure  2  show  that  this  situation 
does  exist  for  all  of  the  alloys  with  the  single  exception  of  the 
aluminum  30  wt.  \  zinc  immediately  after  the  quench.  Porod  radii 
taken  from  the  Al-301  Zn  in  this  condition,  therefore,  have  not 
been  included  in  the  present  analysis. 

The  volume  fraction  of  segregated  zones  and  the  composition 
of  the  zones  during  the  course  of  aging  can  be  monitored  by  meas¬ 
uring  the  change  in  integrated  small-angle  scattering  intensity. 

This  measurement  was  also  performed  by  B  and  G  and  they  show  that 
the  total  integrated  intensity  Q(0)  does  indeed  change  for  all  of 
the  alloys  up  to  100  minutes  or  so  after  quenching.  They  explain, 
however,  that  this  increase  in  Q(0)  represents  an  increase  in  the 
percent  of  zinc  in  the  zones  and  not  an  increase  in  the  volume 
fraction  of  the  zones.  If  their  explanation  is  accepted,  then  the 
calculated  volume  fraction  (6)  (see  Table  I)  can  be  used  from  the 


4  Gerold  (Physica  Status  Solidi  1,  37  1961)  discusses  and  derives 
the,  Guinier  radius,  the  Porod  radius  and  the  integrated  intens¬ 
ity.  The  Guinier  particle  radius  is  obtained  from  the  slope  of 
a  plot  of  log  I  versus  whereas  the  Porod  radius  is  obtained 
from  the  outer  or  tail  portion  of  the  SAXS  curve.  The  Porod 
Tadius  is  a  surface  area  sensitive  measurement. 
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early  stages  of  decomposition  to  obtain  the  number  of  zones  per 
unit  volume  in  each  size  class  "i"  once  the  particle  size  distri¬ 
bution  curve  has  been  measured.  The  analysis  of  particle  size 
distribution  follows  that  given  by  Harkness  and  Gould  in  a  recent 
publication. ^ 

Utilizing  the  Guinier  and  Porod  radii  in  Table  I,  particle 
size  distributions  were  calculated  at  each  aging  time  for  each 
alloy.  The  results  of  such  calculations  are  found  in  Figures  1-3. 
Visual  observations  of  these  figures  indicates  changes  in  the  posi¬ 
tion  of  the  size  distribution  maximum,  which  does  not  coincide  with 
either  the  Guinier  or  Porod  radius.  In  Figures  4-6,  the  total  num¬ 
ber  of  particles  per  unit  volume,  the  mean  zone  radius  and  the 
change  in  the  breadth  (FWHM)  of  the  particle  size  distribution  are 
plotted  as  a  function  of  aging  time.  In  Figures  7  and  8,  the 
growth  paths  of  a  particular  size  particle  can  be  followed  as  a 
function  of  aging  time.  These  curves  will  be  analyzed  in  the  dis¬ 
cussion  section  which  follows. 


Determination  of  NV(R) 


The  log  normal  size  distribution*  as  given  below: 


=  V(2it)  R  Inc 


exp 


-»  f 


lny  -  In  R 
lna 


is  completely  characterized  when  y^  and  a  are  known.  It  has  been 


lny  is  the  geometric  mean  and  ln*a  is  the  geometric  variance  of 
the  log  normal  distribution 
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^  Total  nuaber  of  zones/cm  as 
a  function  of  aging  time  at 
-45*C. 
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Fig.  5  Mean  zone  radius  as  a  function  of  aging  tine  at  -45*C. 
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Fig.  7  The  number  of  particles  per  unit  volume  greater  than  size  i(N>>)  versus  size  i 
as  a  function  of  aging  time  for  Al-20  x/o  Zn  aged  at  -45*C.  v 
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Particle  growth  curves 
for  G.P.  zones  in  Al-20  w/o 
Zn  alloy  aged  at  -45*C. 
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2 

shown  that  y  and  a  may  be  written  in  terms  of  and  Rp.  F(R) 
then  gives  the  form  of  the  experimentally  determined  size  distibu- 
tion  curve  and  one  needs  only  the  volume  fraction  B  of  particles 
to  put  the  curves  on  an  absolute  scale,  i.e.,  Ny(R)  [number  of 
particles  per  cm^  of  radius  R] .  The  volume  fraction  (B)  of  G.P. 
zones  is  given  by  B  and  G“*  for  each  alloy  (see  Table  I)  and  is 
assumed  constant  throughout  aging.  Now 

B  »  4/3*  Z  Ny(R.)  R* 

and 

F(R)  -  k  NyCR). 

Thus , 

B  -  4/3*  k  (FCRj)  R*  ♦  F(R2)  Rlj  ♦  ...  F(Rn)  R*}. 

A  program  has  been  written^  in  PL  1  language  which  calculates 
k  and  provides  the  number  of  particles  in  each  size  class  R.  as 
well  as  the  total  number  of  particles.  Input  consists  of  B,  Rq, 
Rp,  the  initial  value  of  R^ ,  AR  and  the  final  value  of  Rj. 


Discussion 


The  zone  radius  distribution  curves  for  all  three  alloys  are 
shown  in  Figures  1-3.  The  evolution  of  the  size  distribution  at 
-45°C  can  easily  be  follovred.  Note  that  the  mean  zone  radius  of 
the  Al-10  w/o  Zn  alloy  is  constant  over  the  entire  period  while  the 
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curves  for  Al-20  w/o  Zn  and  Al-30  w/o  Zn  alloy  clearly  indicate 
particle  coarsening  (Figure  5)  and  the  total  number  of  particles 
is  decreasing  as  aging  occurs  (Figure  4) .  Thus  the  decrease  in 
the  total  number  of  particles  with  aging  time  indicates  that  smal¬ 
ler  zones  are  dissolving  in  favor  of  the  growth  of  larger  zones. 

Note  that  the  Al-10  w/o  Zn  alloy  is  undergoing  a  redistribu¬ 
tion  process,  although  this  fact  is  not  reflected  in  the  change 
of  mean  zone  radius  with  aging  time  as  shown  in  Figure  5.  The 
small  maximum  of  the  Al-10  w/o  Zn  curve  in  Figure  4  at  approximately 
30  minutes  aging  time  would  appear  to  indicate  that  this  alloy  may 
not  have  fully  decomposed  during  the  quench.  The  change  in  the 
breadth  of  the  zone  size  distribution  curve  for  each  alloy  as 
measured  at  half  maximum  intensity  appears  to  pass  through  a  weak 
maximum  (Figure  6)  indicating  that  the  quenched  state  and  the  well 
aged  state  are  more  monodisperse  than  intermediate  states  of  growth. 
The  zero  aging  times  were  not  included  in  this  figure  as  the  valid¬ 
ity  of  the  Porod  radius  may  be  questioned  for  this  initial  state. 
However,  the  Al-20  w/o  Zn  alloy  appears  to  have  a  sharp  zone -matrix 
interface  even  immediately  after  the  quench.  The  FWHM  breadth  of 

O 

this  curve  can  be  measured  from  Figure  2  to  be  approximately  l.S  A, 
indicating  a  very  monodisperse  initial  zone  state. 

If  the  zone  radius  distribution  curves  are  available,  one  is 
able  to  calculate  the  growth  paths  of  various  sized  zones.  This 
method  has  been  illustrated  by  Harkness  and  Gould^  for  Al-Ag  and 
Al-Zn  alloys  and  the  reader  is  referred  to  that  paper  for  details. 
Ny>  ^  curves  can  be  obtained  from  the  particle  size  distribution 

4  N  >  curves  plot  the  number  of  zones  larger  than  a  particular  zone 
size  versus  that  size. 
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data  if  the  volume  fraction  of  particles  is  known.  A  typical  Ny> 
curve  for  the  Al-20  w/o  An  alloy  is  reproduced  in  Figure  7.  The 
resulting  growth  path  curves  are  shown  in  Figure  8.  As  anticipated, 
the  critical  zone  diameter,  D  critical,  tends  to  increase  with  in¬ 
creasing  aging  time. 

Su— ary 

A  method  for  obtaining  particle  size  distribution  curves  from 
SAXS  data  has  been  applied  to  G.P.  zones  in  the  Al-Zn  system  aged 
at  -45°C.  A  series  of  G.P.  zone  diametei  distribution  curves  have 
been  derived  using  Guinier  and  Porod  radii  measured  by  Bonfiglioli 
and  Guinier.  These  curves  were  then  analyzed  and  the  low  tempera¬ 
ture  G.P.  zone  growth  behavior  can  be  described  as  a  typical  part¬ 
icle  coarsening  process  with  competitive  growth.  It  appears  that 
the  zonal  state  immediately  after  the  quench  is  characterized  by 
a  more  monodisperse  zone  size  distribution  than  at  later  aging  times. 
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V.  Magnetic  Materials  (J.  K.  Watson) 

A.  ANALOG  APPROXIMATIONS  OF  TRANSIENT  FLUX  DIFFUSION  (J.  K.  Watson  and  C.  T.  Leis  ) 
Introduction 

The  trapped  flux  interaction  between  magnetic  films  and  nearby  elec** 

trical  conductors  is  a  well-known  phenomenon  experimentally  ^  and  has 

35 

been  used  as  the  basic  mechanism  of  NDRO  memories  '  .  However,  a  theoretical 

4  5 

understanding  of  this  process  has  been  difficult.  Previous  analyses  '  are 
formidable  enough  computationally  and  conceptually  to  limit  their  usefulness. 

The  new  approach^  presented  here  has  a  surprisingly  simple  form  and  leads  to 
results  that  agree  well  with  more  complicated  methods^.  Our  approximations 
and  assumptions  are  indicated  in  detail,  since  this  analog  method  is  thought 
to  be  applicable  to  other  magnetic  diffusion  problems. 

The  underlying  mechanism  stems  from  Maxwell's  equations  which  require 
that  changes  of  flux  density  B  within  a  non-magnetic  conductor  must  satisfy 
the  diffusion  equation 

V2  B  -  y  o  3B/3t  (1) 

o 

where  is  the  permeability  of  space  and  a  is  the  conductivity  of  the 
material.  In  the  geometry  of  interest  here,  B  is  a  complicated  spatial  func¬ 
tion  which  is  not  easily  represented  in  closed  form,  even  in  the  static  case 
where  it  defines  the  stray  field  from  the  film. 


*C.  T.  Leis  is  now  with  the  Hewlett-Packard  Co.,  Cupertino,  California. 
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When  a  film  is  switched,  or  remagnetized  to  another  orientation,  the 
stray  field  must  also  be  reoriented.  The  transient  interaction  exists  while 
the  old  field  collapses  and  the  new  field  builds  up.  The  major  objective  of 
this  paper  is  to  provide  a  basis  for  estimating  the  magnitude  and  time  duration 
of  the  transient,  taking  into  account  the  geometry  of  the  film, insulators, 
and  conductors. 

The  general  approach  which  follows  is  based  on  the  analogy  to  (1)  of 
an  RC  transmission  line,  justified  by  assumptions  described  below.  Subse¬ 
quent  use  of  the  analog  permits  geometrical  factors  to  be  represented  as 
electric  circuit  components,  with  certain  conceptual  advantages  that  will  be 
evident. 

System  Model 

In  this  section  we  consider  those  assumptions  which  lead  to  the  repre¬ 
sentation  of  the  3-dimensional  nonuniform  field  problem  by  a  conceptually 
simpler  one-dimensional,  uniform  field  approximation. 

For  convenience  of  developing  the  concepts,  the  geometry  of  Fig.  1(a)  is 
assumed  in  which  the  magnetic  film  is  halfway  between  two  conducting  plates 
of  a  stripline. 


The  total  stray  flux,  $,  emanating  from  the  film  has  a  free  component,  4^, 
and  a  trapped  component,  4>t ,  where  the  latter  penetrates  the  conductors  in  the 
static  situation.  Immediately  after  an  instantaneous  reversal  of  film  magneti¬ 
sation,  Fig.  1(b),  still  retains  its  original  orientation  within  the  con¬ 
ductors.  However  the  total  flux  ^  io  temporarily  squeezed  between  the  strip¬ 
lines  in  the  new  orientation,  along  with  the  return  from  the  trapped  flux, 
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Fig.  1.  Geometrical  notation,  and  flux 
patterns  around  thin  film.  (a)  Before 
switching.  (b)  After  switching. 


because  of  the  relatively  slow  diffusion. 

We  consider  the  flux  density  in  the  vicinity  of  the  z-axis,  just  above 
and  below  the  film  which  is  in  the  x-y  plane.  Here  the  flux  is  parallel  to 
the  plane  of  the  film  and  conductors.  The  flux  density  along  the  z-axis  has 
only  an  x-component,  which  is,  however,  a  function  of  z  and  time.  Hence  the 
problem  of  interest  is  reduced  to  one  spatial  dimension.  For  a  circular  film 
of  radius  r,  thickness  6,  and  saturation  magnetization  M,  the  static  flux 
density  along  the  z-axis  is  given  by 


B(z) 


dr2M 


.  ,  2  .  2.3/2 

4(r  +  z  ) 


(2) 


Interestingly,  it  is  not  necessary  to  find  a  diffusion  process  which  approaches 
(2)  as  a  final  value.  Instead,  we  shall  use  a  surprisingly  simple,  rectangular 
approximation  to  (2),  based  on  the  key  concept  that  the  total  stray  flux  and 
the  trapped  fraction  of  flux  are  the  pertinent  quantities  of  the  diffusion 
process.  Using  the  geometry  of  Fig.  1,  with  insulators  near  the  film  for 
| z |  <  h  and  conductors  for  |z|  >  h,  the  free  flux  is  defined  in  terms  of  its 

average  value  B  as  . 

a  "  22  -1/2 

-  2h  B  -  /  B(z)  dz  «  2h  Q(r  +  h*  ) 

1  a  -h 

where  Q  ■  6M/^.  Choosing  the  rectangular  approximation  of  Fig.  2,  of  value 
B#  inside  the  conductor  and  out,  gives  rise  to  the  concept  of  "effective 
penetration  depth"  d. 


Thus  the  trapped  flux  is 

♦-2  /  B(z)  dz  *  2d  B 

, -  h  a 

•y  2  2 

where  d  *  ir  +h  -h.  In  the  transmission  line  analogy  to  be  developed 
subsequently,  d  will  be  a  useful  measure  for  defining  the  effective  termi- 
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nation  impedance  for  the  lines.  The  important  role  of  the  conductor  sepa¬ 
ration  2h  can  be  discerned  from  t lie  equation  for  parameter  d. 

The  original  problem  has  now  been  reduced  to  a  much  simplified  form. 
The  assumed  flux  patterns  are  illustrated  in  Fig.  3. 


Before  switching,  the  flux  density  along  the  z-axls  Is  zero  except  for 

a  constant  value  of  -B  from  -(h  +  d)  <  z  <  (h  +  d). 

a 

Just  after  a  fast  rotation  of  M  the  entire  flux  from  the  film  plus  the 
closing  path  of  the  trapped  flux  is  squeezed  between  the  plates.  The  new 
value  of  flux  density  in  the  region  is  then 

Bs  "  h(f  +  “  Ba<h  +  2d)/h*  U) 

As  time  progresses,  flux  from  the  space  between  the  plates  diffuses  into  and 
through  the  conductors  until  equilibrium  is  reached  in  the  opposite  direction. 
At  all  times  the  total  flux 


4  -  /  B(z)  dz  «  2Ba(h  +  d) 


must  stay  constant. 

The  model  for  the  physical  system  is  now  complete.  Of  particular  interest 
are  the  two  components  of  flux  between  the  plates  during  the  diffusion  process. 
They  represent  the  mechanism  of  the  interaction  torque  between  film  and  strip- 
line. 
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z — ► 


-h-d 


-h 


h 

-+ 


h+d 


Fig.  3.  The  assumed  s-axis  flux  density 
distributions.  (a)  Before  switching. 

(b)  Immediately  after  "fast"  switching. 

(c)  After  diffusion. 
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Electric  Analog 


The  assumptions  cited  above  permit  the  diffusion  equation  to  be  re¬ 
duced  to  one  dimension  along  the  z-axis,  thereby  providing  access  to  solu¬ 
tions  which  are  available^.  For  our  purposes  it  is  instructive  to  consider 
the  analogous  RC  transmission  line  equation 

D  3^B/3z^  *  3B/3t  analogous  to  j^V/3z^  «  RC  3V/3t 

where  diffusivity  D  ■  1 /pqo.  In  this  comparison  flux  is  analogous  to  elec¬ 
tric  charge  as  the  diffusing  quantity.  Outlined  in  the  appendix  and  illus¬ 
trated  in  Fig.  4,  the  linear  dimensions  along  the  z-axis  are  analogous  to 
circuit  elements  for  y-dimensions  normalized  to  unity. 


In  this  figure,  the  central  capacitor  CQ  corresponds  to  the  space  between 
the  plates  so  its  value  is  the  z-dimension  2h.  As  shown,  each  stripline 
plate  is  represented  by  RC  transmission  line  of  length  a^  and  each  is  termi¬ 
nated  by  a  capacitor  Cj  ■  d  -  a^  which  represents  the  amount  of  free  space 
beyond  the  conductor  available  for  flux  distribution.  If  the  conductor  thick¬ 
ness  is  greater  than  d,  the  transmission  line  will  be  an  open-circuited  line 
of  length  d. 

We  will  be  concerned  with  finding  V(t),  the  voltage  across  C^,  which  is 
equivalent  to  finding  the  flux  density  between  the  conductors  as  a  function  of 
time.  As  suggested  by  Fig.  3,  the  initial  conditions  after  an  instantaneous 
film  reversal  would  place  a  positive  voltage  B#  on  CQ,  and  a  voltage  -B^  on 
all  the  other  capacitors.  At  time  t  ■  0  the  system  will  be  released  and  cur¬ 
rent  will  flow  until  equilibrium  is  reached. 
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Conventional  circuit  analysis  techniques,  including  lumped  approximations 
to  distributed  lines,  can  now  be  brought  to  bear  on  the  problem.  Results  from 
two  special  eases  are  considered,  to  illustrate  the  principle  of  finding  the 
Input  impedance  which  is  connected  to  C^.  The  feasibility  will  be  evident  of 
extending  the  technique  to  other  geometry. 

Case  1.  a  «  d 

When  the  conductors  are  very  thin,  the  capacitance  of  the  RC  line  can 
be  neglected  with  little  loss  of  accuracy  with  the  result  shown  in  Fig.  5. 

Each  resistor  has  the  value  a/D,  from  the  appendix,  and  the  capacitor  has  the 
value  d  -  a  *»  d.  Thus  the  input  impedance  is 

Z(s)  ■  a/D  +  1/ds 


Case  2.  d  -  a 

When  the  conductors  are  very  thick,  the  electric  circuit  consists  of  an 

open  circuited  RC  line  of  length  d  connected  on  each  side  of  C  .  The  input 

o  * 

Impedance  for  such  a  line  is 

Z(s)  -  —  coth  (d 

yDs 

as  noted  in  the  appendix. 

It  is  sometimes  useful  to  consider  a  lumped  approximation  to  the  RC  line. 

for  a  circuit  model  in  the  same  form  as  Fig.  5.  Using  the  two  lowest-order 

8 

terms  of  the  series  expansion  about  zero  of  coth,  which  is  equivalent  to  a 
low-frequency  model^  for  Z,  the  result  is 


Z(s)  * 


d_ 

3D 


+ 


i_ 

ds 
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Pig.  5".  Complete  analog  circuit  incorporating 
approximations  for  very  thick  or  very  thin 
strip-line  plates. 
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System  Solutions 


We  now  indicate  transient  solutions  to  Llie  analogous  circuits  which  were 
developed  In  the  previous  section,  using  the  approximation  of  Fig.  5.  We 
consider  three  cases:  both  stripline  conductors  are  very  thin,  both  conductors 
arc  very  thick,  and  one  is  thick  and  the  other  is  thin.  In  each  case  the  fact 
that  and  have  identical  initial  conditions  stems  from  the  assumption 
that  the.  film  is  located  midway  between  the  conductors.  Using  conventional 
circuit  analysis  techniques,  the  transform  of  the  voltage  is 


B  +  B  7, 

•  ■  /  \  8 fl  r  1 

V(s)  - - - —  [— 


cl +  k 


B 

]  -  -a- 
J  S 


where  is  the  parallel  combination  of and 
the  same  thickness,  the  final  result  is 


If  both  conductors  have 


B(t)  *=  B  +  [B  -  B  ]  c"t/r 

8  b  8 

after  replacing  the  electrical  symbols  by  their  physical  counterparts.  In 
that  equation  B  is  given  by  (3)  as  before.  The  value  of  time-constant  \ 
ranges  from  t  =  ah'/D  for  very  thin  conductors  to  approximately  t  <=  dh'/3U 
for  very  thick  conductors,  where  p^o  =  1/D  as  before,  h'  ■  dh/d  +  h  *=  h/1  +  X, 
where  X  is  the  ratio  of  free  to  trapped  flux. 

A  particularly  interesting  geometry  consists  of  one  thin  conductor  above 
the  film  and  a  thick  conducting  ground  plane  below.  The  circuit  of  Fig.  5  is 
then  solved  for  V(t)  when  R^  <<  R^.  Two  important  time  constants  result 

^  ■  R^K^  with  given  by  Cq  in  series  with  C^, 

and  ^  e  ^2^2  S*vcn  ^y  ^2  *n  scr*cs  +  C^. 

The  physical  interpretation  is  that  the  flux  diffuses  rapidly  through  the 
thin  plate  until  a  temporary  equilibrium  is  reached.  Then  the  slow  diffusion 
into  the  ground  plane  takes  place.  To  find  the  intermediate  value  of  flux 
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density  B^  it  is  assumed  that  the  diffusion  into  the  thick  conductor  is  negli 
gible  during  the  fast  mode.  Using  a  sketch  of  the  above  statement,  and  the 
principle  that  total  flux  (4)  must  remain  constant,  it  follows  that 


B 


i 


B  U.t- 2A). 
a(l  +  2A) 


The  complete  estimate  of  the  time  response  of  flux  density  between 


conductors  is 

-t/r.  -t/t- 

B(t)  =■  B  +  [B  -  B . )  e  1  +  (B,  -  B  ]  e 
a  s  l  la 


the 


(5) 


where  the  two  terms  on  the  right  are  due  to  the  influence  of  eddy  currents. 


The  form  of  the  eddy  current  component  is  the  same  ns  an  empirical  result  . 

2  ah 

In  physical  units  the  two  dominant  time  constants  are  ■  DO~+  2T) 


and 


ciCl  +  2A) 

t2  *=  -u&rrx)  w,icrc  a  *s  t*lc  thickness  of  the  top  conductor.  Fig.  A 


6D(1  d  A) 

shows  a  sketch  of  (?),  which  we  now  examine  in  gre  ter  detail. 


Comparison  of  Results 

In  an  attempt  to  establish  a  measure  of  the  validity  of  (5)  and  its 

associated  time  constants,  a  comparison  was  made  with  some  results  of  Jutzi^. 

Using  a  complicated  numerical  procedure,  he  obtained  solutions  for  the  magnetic 

field  intensity  H  =  B/p^  between  the  conductors  as  a  function  of  time  after 

an  instantaneous  change  in  the  magnetization  vector  from  a  value  of  zero  to 
2 

1  W/m  .  His  procedure  differs  from  the  above  analysis  in  having  no  trapped 
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Fig.  6  Estimation  of  the  time  response  of  flux 
density  between  the  conductors  (one  thin  plate, 
the  other  thick) . 
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flux  in  tlto  initial  condition,  thus  the  initial  (li  )  and  the  intermediate  (II. ) 

S  X 


values  of  flux  density  need  to  be  adjusted.  The  time  constants  remain  un¬ 
affected  since  they  depend  solely  on  the  geometry.  The  physical  dimensions 
in  the  comparison  are:  <5  «  .07  y,  r  B  200  p,  h  c  2.5  p,  and  a  ■  2.5  p.  The 

ratio  of  h/r  implies  X  '  0  and  d  s  200  p.  The  resulting  time  constants  are 
Tj  *  0.915  ns  and  T2  B  487  ns.  After  omitting  the  trapped  flux,  the  field 
intensities  are  computed  to  be  I!  B  B  /p  ■  70  A/m;  H,  *  140  A/nt;  II  ■  5600  A/m. 

8  8  O  1  8 

The  fast  portion  and  the  slow  portion  of  the  transient  are  compared 
separately  with  Jutzi's  calculated  results  in  Figs.  ^  and  8,  respectively,  and 
the  agreement  is  good.  It  is  concluded  that  (5)  is  an  adequate  representation 
of  the  physical  situation. 


Conclusion 

This  paper  has  sought  to  establish  an  approximation  method  for  flux 
diffusion,  then  to  apply  it  to  the  computation  of  parameters  which  are  of 
considerable  practical  interest.  These  objectives  motivated  our  attempt  to 
carefully  define  the  assumptions  and  approximations  of  the  system  model.  The 
approach  is  attractive  because  of  the  simplicity  of  the  results,  the  ease  with 
wnich  the  geometry  is  considered,  and  certain  conceptual  advantages  which  are 
gained  from  the  analog  representation  of  the  problem.  The  circuit  analog  is 
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H(t) 


Fig.  7  Fast  decay  mode  comparison.  Jutzi's 
computation ,  A,  and  the  estimation,  B,  with 
time  constant  ij  =  .915  ns. 
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particularly  helpful  to  the  electrical  engineer. 

10  , 

The  model  developed  here  provides  a  basis  for  the  observation  that 
eddy-current  torques  have  not  been  a  significant  factor  in  film  flux-reversal 
studies  where  the  stripline  separation  is  large. 

The  question  of  accuracy  has  not  yet  been  fully  answered  due  to  the 
difficulty  of  measurements,  but  these  results  compare  well  with  computations 
of  others^. 

APPENDIX 

ELECTRICAL  ANALOGY  FOR  DIFFUSION 
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The  equation  for  an  RC  transmission  line  is  well-known  to  be  — j  =  RC 
which  has  an  evident  analogy  to  D  — r  *=  —  on  the  basis  of  the  similar 

„  A  01 
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equations.  The  purpose  here  is  to  further  interpret  the  parameters  of  resis¬ 
tance  and  capacitance  per  unit  length  in  terms  of  line  geometry. 

First,  it  is  argued  on  the  basis  of  the  following  tabulation  that  capaci¬ 
tance  is  a  spatially  related  parameter  which  is  independent  of  conductivity, 
and  which  serves  to  contain  a  diffusing  substance.  Analogous  quantities  are 
vertical. 

(capacitance)  y  (voltage)  «  (charge) 

(area)  Y  (flux  density)  *  (quantity  of  flux 

On  that  basis  we  take  capacitance  per  length  as  analogous  to  area  per  length, 

where  the  total  area  is  yz,  normal  to  the  x-direction  of  flux.  Taking  length  y 

as  unity  leads  to  unit  dimension  for  C,  hence  R  -  1/D.  Applying  these  defini- 
io  ttvory;  ( - 

tions^  y  *  ^  RCs  and  Zq  =  ^  R/C $,  and  to  the  equation  for  input  impedance  of 

an  open-circuit  line  of  length  d,  Z,  =  Z  coth  Yd»  leads  to  the  equation 

in  o 

of  case  2, 
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VI.  Discussion 


The  research  of  this  fifth  semiannual  period  has  yielded  results 
of  technical  significance  and  it  has  laid  the  foundations  for  future 
study  and  growth. 

The  study  of  the  photomagnetoelectric  effect  in  silicon  doped  with 
gold  and  phosphorus  has  yielded  a  basic  characterization  of  such  material 
and  has  suggested  several  potential  applications.  Among  these  is  the 
possibility  of  a  laser  beam  demodulator  utilizing  silicon  overcompensated 
with  gold,  with  applications  to  the  field  of  communications.  Research  ex¬ 
ploring  this  possibility  is  now  proceeding. 

From  the  study  of  the  density  of  quantum  states  in  a  semiconductor  has 
come  a  basic  understanding  of  the  dependence  on  impurity  concentration  of 
this  density.  The  impurity  concentration  dependence  can  manifest  itself 
importantly  in  the  behavior  of  the  highest-performance  transistors  of  the 
present  day,  and  thus  the  study  has  yielded  information  relevant  to  the 
many  fields  of  application  that  utilize  transistors:  communications, 
telemetry,  detection,  navigation,  control,  bio-medicine,  etc. 

The  electric  field  that  a  ceramic  can  withstand  may  determine  its 
applicability  in  various  uses  in  the  electronic  technology:  as  an  insulator, 
as  a  protective  surface  in  various  devices,  as  an  isolator  in  integrated 
circuits,  and  as  a  material  for  encapsulation.  The  study  reported  here 
suggests  the  relative  importance  of  the  many  variables  pertinent  to  the 
design  of  ceramic  materials  for  use  in  these  applications. 

From  the  investigation  of  the  relationship  borne  between  heterogeneities 
present  in  semiconducting  glass  and  the  consequent  electrical  properties 
has  come  an  understanding  that  will  contribute  ultimately  to  the  evaluation 
of  the  potential  merit  of  semiconducting  glass  as  a  device  material  for 
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various  military  applications,  particularly  for  those  involving  system 
operation  in  high-radiation  environments. 

The  trapped-flux  interaction  between  magnetic  films  and  nearby  elec¬ 
trical  conductors  is  a  well-known  phenomenon  experimentally  and  has  consti¬ 
tuted  the  basic  mechanism  of  NDRO  memories.  Its  study  therefore  relates 
directly  to  the  design  of  digital  computers  and  thus  to  many  fields  of 
relevance  to  the  D.O.D.  Because  of  its  simplicity  relative  to  methods  of 
analysis  previously  proposed,  the  new  approach  reported  here  to  the  analysis 
cf  this  phenomenon  promises  a  substantial  advantage  to  the  designer  of 
magnetic  memories. 

The  plan  of  the  research  remains  unchanged.  Emphasis  will  continue 
to  focus  on  the  relationship  between  the  structural  properties  of  materials, 
on  the  one  hand,  and  the  resultant  electrical,  quantum,  chemical,  and  thermal 
properties,  on  the  other.  Device  applications  made  possible  by  the  under¬ 
standing  gained  from  this  materials  research  will  receive  increasing  attention. 
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IS.  abstract  This  report  describes  technical  findings  in  four  main  subject 
areas:  semiconductors  and  semiconductor  devices,  insulating  and  semi¬ 
conducting  glasses,  measurement  techniques,  and  magnetic  materials. 

For  silicon  doped  with  gold  and  phosphorus,  a  study  °f  the  photo¬ 
magnetoelectric  effect  yields  a  basic  theoretical  relationship  between 
the  excess  hole  and  excess  electron  concentrations  and  an  experimental 
determination  of  the  effective  lifetime  for  temperatures  between  21  and 
84°K.  A  pragmatic  approach  is  proposed  and  empirically  tested  con¬ 
cerning  tne  dependence  that  the  density  of  quantum  states  in  a  semicon¬ 
ductor  shows  on  the  impurity  concentration.  For  a  variety  of  ceramic 
materials,  the  electric-field  breakdown  shown  is  related  to  the  various 
environmental  and  material  variables.  Based  on  available  data  concern¬ 
ing  material  structure,  a  study  is  made  to  characterize  the  relationshij 
borne  in  semiconducting  glasses  between  heterogeneity  in  structure  and 
resultant  electrical  properties.  Errors  in  x-ray  spectrochemical  anal¬ 
ysis  of  powdered  materials  due  to  pressing  variations  are  quantitatively 

examined.  The  several  means  for  analysis  of  the  defect  structure  of 
solid  surfaces  are  discussed  together  with  the  relation  that  this  de¬ 
fect  structure  bears  to  the  surface  electronic  and  chemical  properties. 

Aat- a  r»n  t-Vwa  1  nu-fomnprflfiirp  ao-fno  of  aluminum— zinc 
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abd  evaluated  for  characterizing  the  trapped-flux  interaction  between 
magnetic  films  and  nearby  electrical  conductors. 
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